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■ ABSTRACT 

A new class of giant resonances in nuclei, namely double giant resonances, is discussed. 
^ I They are giant resonances built on top of other giant resonances. Investigation on their 

I properties, together with similar studies on low- lying two-phonon states, should give an 

' answer on how far the harmonic picture of boson-type excitations holds in the finite 

fermion systems like atomic nuclei. 
^ ■ The main attention in this review is paid to double giant dipole resonances (DGDR) 

' which are observed in relativistic heavy ion collisions with very large cross sections. A 

■ great experimental and theoretical effort is underway to understand the reaction mecha- 

00 ' ■ ■ 

■ nism which leads to the excitation of these states in nuclei, as well as the better microscopic 

0^ ' understanding of their properties. The Coulomb mechanism of the excitation of single and 

double giant resonances in heavy ion collision at different projectile energies is discussed in 
details. A contribution of the nuclear excitation to the total cross section of the reaction 
is also considered. The Coulomb excitation of double resonances is described within both, 
^ ■ the second-order perturbation theory approach and in coupled-channels calculation. The 

properties of single and double resonances are considered within the phenomenologic har- 
monic vibrator model and microscopic quasiparticle-RPA approach. For the last we use 
the Quasiparticle-Phonon Model (QPM) the basic ideas and formalism of which are pre- 
sented. The QPM predictions of the DGDR properties (energy centroids, widths, strength 
distributions, anharmonicities and excitation cross sections) are compared to predictions 
of harmonic vibrator model, results of other microscopic calculations and experimental 
data available. 
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1 Introduction 



The phenomenon of a Giant Resonance (GR) in a nucleus is now known for more than 60 
years. The first article on this subject was published in 1937 by Bothe and Gentner who 
observed an unexpectedly large absorption of 17.6 MeV photons (from the ^Li(p, 7) reaction) 
in some targets. They noticed that the cross section for ®^Cu was surprisingly high and they 
suggested that this might be due to a resonance phenomenon. These observations were later 
confirmed by Baldwin and Klaiber (1947) with photons from a betatron. In 1948 Goldhaber 
and Teller p[ interpreted these resonances (named by isovector giant dipole resonances) with a 
hydrodynamical model in which rigid proton and neutron fluids vibrate against each other, the 
restoring force resulting from the surface energy. Steinwendel and Jensen later developed 
the model, considering compressible neutron and proton fluids vibrating in opposite phase 
in a common fixed sphere, the restoring force resulting from the volume symmetry energy. 
The standard microscopic basis for the description of giant resonances is the Random Phase 
Approximation (RPA) in which giant resonances appear as coherent superpositions of one- 
particle one-hole (Iplh) excitations in closed shell nuclei or two quasi-particle excitations in 
open shell nuclei (for a review of these techniques, see, e.g., Ref. @]). 

The isoscalar quadrupole resonances were discovered in inelastic electron scattering by Pit- 
than and Walcher (1971) and in proton scattering by Lewis and Bertrand [^. Giant monopole 
resonances were found later and their properties are closely related to the compression mod- 
ulus of nuclear matter. Following these, other resonances of higher multipolarities and giant 
magnetic resonances were investigated. Typical probes for giant resonance studies are (a) 7's 
and electrons for the excitation of GDR (isovector giant dipole resonance), (b) a-particles and 
electrons for the excitation of isoscalar GMR (giant monopole resonance) and GQR (giant 
quadrupole resonance), and (c) {p,n), or ('^He, t), for Gamow- Teller resonances, respectively. 

Relativistic Coulomb Excitation (RCE) is a well established tool to unravel interesting 
aspects of nuclear structure 0. Examples are the studies of multiphonon resonances in the SIS 
accelerator at the GSI facility, in Darmstadt, Germany |^ |, p. Important properties of nuclei 
far from stability |T0| , |Tl| have also been studied with this method. 

Inelastic scattering studies with heavy ion beams have opened new possibilities in the field 
(for a review the experimental developments, see Ref. 0, |^). A striking feature was observed 
when either the beam energy was increased, or heavier projectiles were used, or both |T^]. 
This is displayed in Fig. 0, where the excitation of the GDR in ^°^Pb was observed in the 
inelastic scattering of ^^O at 22A MeV and 84A MeV, respectively, and ^"^Ar at 95A MeV 



p!3| , 14|. What one clearly sees is that the "bump" corresponding to the GDR at 13.5 MeV is 



appreciably enhanced. This feature is solely due to one agent: the electromagnetic interaction 
between the nuclei. This interaction is more effective at higher energies, and for increasing 
charge of the projectile. 

Baur and Bertulani showed in Ref. |T^ that the excitation probabilities of the GDR in 
heavy ion collisions approach unity at grazing impact parameters. They also showed that, if 
double GDR resonance (i.e. a GDR excited on a GDR state) exists then the cross sections 
for their excitation in heavy ion collisions at relativistic energies are of order of a few hundred 
of milibarns. These calculations motivated experimentalists at the GSI ^ and elsewhere 



p!3| , |Tj] to look for the signatures of the DGDR in the laboratory. This has by now become a 



very active field in nuclear physics with a great theoretical and experimental interest & H. 



In the first part of tliis review, we study tfie reactions meclianism in RCE of giant reso- 
nances in several collisions between heavy ions. In section 2 we start with the description of the 
semiclassical theory for relativistic Coulomb excitation, then we consider the effects of recoil 
and later we describe the fully quantum mechanical approach. The role of nuclear excitation in 
relativistic heavy ion collisions is also discussed here. We demonstrate that the experimental 
data on the excitation and decay of single giant resonances are well described by these for- 
malisms. In section 3 the process of the excitation of multi-phonon resonances in relativistic 
heavy ion collisions is considered within the second-order perturbation theory and in coupled- 
channels calculations. Giant resonances are treated in this section within the phenomenologic 
harmonic vibrator model. Some general arguments for the width of multi-phonon resonances 
are discussed here as well as an influence of giant resonances width on the total cross section 
of their excitation. A good part of this report (sections 4 and 5) is dedicated to a review of the 
microscopic properties on the giant resonances in the Quasiparticle-Phonon Model. In section 
4 we present the main ideas and formalism of this model. The particle-hole modes of nuclear 
excitation are projected into the space of quasi-bosons, phonons, and matrix elements of in- 
teraction between one- and multi-phonon configurations are calculated on microscopic footing 
within this approach. In section 5 we use this model as a basis for a detailed investigation 
of the interplay between excitation mechanisms and the nuclear structure in the excitation of 
the DGDR. Different aspects related to the physical properties of the DGDR in heavy nuclei 
(energy centroids, widths, strength distributions, anharmonicities and excitation cross sections) 
as predicted by microscopic studies are discussed in this section and compared to experimental 
data. 



2 Heavy ion excitation of giant resonances 
2.1 Coulomb excitation at relativistic energies 

In relativistic heavy ion collisions, the wavelength associated to the projectile-target separa- 
tion is much smaller than the characteristic lengths of system. It is, therefore, a reasonable 
approximation to treat r as a classical variable r(t), given at each instant by the trajectory 
followed by the relative motion. At high energies, is also a good approximation to replace this 
trajectory by a straight line. The intrinsic dynamics can then be handled as a quantum me- 
chanics problem with a time dependent Hamiltonian. This treatment is discussed in full details 
by Alder and Winther in Refs. |jl6|, O, |18|. We will describe next the formalism developed by 



Canto et al. [0, which explicitly gives the time-dependence of the multipole fields, useful for 
a coupled-channels calculation. 

The intrinsic state \ip{t) > satisfies the Schrodinger equation 

[h + Vim)] m) >= ^^^^^ • (1) 

Above, h is the intrinsic Hamiltonian and V is the channel-coupling interaction. 

Expanding the wave function in the set {|m >; m = 0, A^} of eigenstates of h, where A^ is 
the number of excited states included in the Coupled-Channels problem, we obtain 

N 

IHi) >= am{t) \m> exp ( - iE^mt/h) , (2) 

m=0 



where Em is the energy of the state \m >. Taking scalar product with each of the states < n\ 
we get the set of coupled equations 



N 



in<Xn{t)=Y. < n|F|m > e*(^"-^'")*/^ a„(t) n = OtoA^. (3) 



m=0 



It should be remarked that the amplitudes depend also on the impact parameter h specifying 
the classical trajectory followed by the system. For the sake of keeping the notation simple, we 
do not indicate this dependence explicitly. We write, therefore, anif) instead of an{b,t). Since 
the interaction V vanishes ast—* ±00, the amplitudes have as initial condition an{t —00) = 
6{n, 0) and they tend to constant values as t — 00. Therefore, the excitation probability of an 
intrinsic state \n > in a collision with impact parameter b is given as 

Pn{b) = Moo)\^ . (4) 

The total cross section for excitation of the state \n > can be approximated by the classical 
expression 

an = 2n I Pnib) bdb . (5) 

Since we are interested in the excitation of specific nuclear states, with good angular mo- 
mentum and parity quantum numbers, it is appropriate to develop the time-dependent couphng 



interaction V{t) into multipoles. In Ref. |jT8[, a multipole expansion of the electromagnetic ex- 
citation amplitudes in relativistic heavy ion collisions was carried out. This work used first 
order perturbation theory and the semiclassical approximation. The time-dependence of the 
multipole interactions was not explicitly given. This was accomplished in Ref. |T^, which we 
describe next. 

We consider a nucleus 2 which is at rest and a relativistic nucleus 1 which moves along the 
z-axis and is excited from the initial state |/jMj > to the state \IfMf > by the electromagnetic 
field of nucleus 1. The nuclear states are specified by the spin quantum numbers Jj, // and 
by the corresponding magnetic quantum numbers Mj, Mf, respectively. We assume that the 
relativistic nucleus 1 moves along a straight-line trajectory with impact parameter b, which is 
therefore also the distance of the closest approach between the center of mass of the two nuclei 
at the time t = 0. We shall consider the situation where b is larger than the sum of the two 
nuclear radii, such that the charge distributions of the two nuclei do not strongly overlap at 
any time. The electromagnetic field of the nucleus 1 in the reference frame of nucleus 2 is given 
by the usual Lorentz transformation [^] of the scalar potential 0(r) = Zie/\r\, i.e., 

<P{r',t) = 7 0[b'-b, -fiz'-vt)_, 

A{r',t) = ^ 7 0[b'-b, 7(^'-t;t)]. (6) 

Here b (impact parameter) and b' are the components of the radius- vectors r and r' transverse 
to V. 

The time-dependent matrix element for electromagnetic excitation is of the form 

Vf.it) =< IfMf\[pir') - ^ ■ J(r')] </)(r', t)|/,M, > , (7) 

where p (J) is the nuclear transition density (current). 

A Taylor-series expansion of the Lienard-Wiechert potential around r' = yields 



0(r', t)=70[r(t)]+7V0[r(t)]-r' + 



(8) 



where r = (b, ^vt) , and the following simplifying notation is used: 

= V>(r', t) 



r'=0 



Thus, 



Vf.it) ^< IfMf\[p{r') - - 



J(r')l f70(r) + 7r' • V(/)(r) 



Using the continuity equation 

V ■ 3 — —i cu p , 
where cu — {Ef — Ei)/h, and integrating by parts, 

-V V 



In spherical coordinates 



70(r + 7r' • V0(r) 



\IiMi > 



where 



and are the spherical unit vectors 

e± = ^-^(ex ±ey). 

We will use the relations 



eo = ez 



r V ^ 



V At: 




v(ryro) 



and 

where L = —iV x r. 
Then, one can write 



V X L(r'^ll^) = t{k + l)V(r'=ll^) 



J • - ^) [70 + 7r' ■ V0 



-7J 



X 



c 



The last term in the above equation can be rewritten as 

c 



2c2 ^ 

'V -r 



eo(r' • V0) + (r' • eo)V0 
eo(r' • V0) - (r' • eo)V0 



(9) 

(10) 
(11) 

(12) 

(13) 
(14) 



(15) 
(16) 



(17) 



(18) 



The first term in this equation is symmetric under parity inversion, and contributes to the 
electric quadrupole {E2) excitation amplitudes, since 



2r-2 



eo(r' • V0) + (r' • eo)V(/) = — J • V z'{r' ■ V0) 



(19) 



The second term in Eq. (|T8D is antisymmetric in J and r', and leads to magnetic dipole (Ml) 
excitations. Indeed, using Eqs. (p!3|-[T^), one finds 



— J 



eo(r' ■ V<P) - (r' ■ eo)V0 



— J 



(20) 



Thus, only the first two terms on the right-hand-side of Eq. ([T7|) contribute to the electric 
dipole (-E1) excitations. Inserting them into Eq. ([12D , we get 



where 



and 




M{E\, -f^) = ^ J d\j{r).V{rY,^) = J p(r) r Y^^{ 



(21) 



(22) 



/3± = -a^ = -{V(t)-e 



9b 



/3o 



The derivatives of the potential (j) are explicitly given by 



r'=0 



r'=0 



— z 7 ft 



-X 6^ 

Zie 



Zie 



[62 +72^2^2] 3/2 



[62 + 72^2^2] 3/2 • 

Using the results above, we get for the electric dipole potential 



/27r 



+ V2 



7 



tSi{t) - S2{t) 



where r = 7^/6, and 



62 [1 + r2]3/2 



and £-2ij) 



A</,(E1,0) , 



h [l + r2]i/2 



(23) 



(24) 



(25) 



(26) 



are the transverse and longitudinal electric fields generated the relativistic nucleus with charge 
Zie, respectively. From the definition 



2c 



(27) 



and Eq. 



we find 



'27r V 



^/f H^) = M/^ 77 ^i(^) A</,,(M1,1) +A^^,(M1,-1) 



(28) 



The current J in Eq. (^) is made up of the usual convective part and a magnetization part, 
proportional to the intrinsic (Dirac and anomalous) magnetic moment of the nucleons. 



To obtain the electric quadrupole {E2) potential we use the third term in the Taylor expan- 
sion of Eq. (P). Using the continuity equation, a part of this term will contribute to E3 and 
M2 excitations, which we neglect. We then find that 



r) 



7 3^3(t) Mf^{E2,2) + Mf^{E2,-2) 



+ 7 



+ 76 7^ 



7C^ 

2t' - l) £,{t) 



6 rSsir) - t — S,{r) Mf,{E2, -1) +Mfi{E2, 1) 



t^rS,iT)]Mf,iE2,0)\, 

7c 



where S^^t) is the quadrupole electric field of nucleus 1, given by 

Zie 



53 [1 + ^2]5/2 • 



(29) 



(30) 



The fields Si{T) peak around r = 0, and decrease fastly within an interval At ~ 1. This 
corresponds to a collisional time At ~ b/'-fv. This means that, numerically one needs to 
integrate the Coupled- Channels equations (Eq. (|^)) only in a time interval within a range 
n X At around r = 0, with n equal to a small integer number. 



2.1.1 First-order perturbation theory 

In most cases, the first-order perturbation theory is a good approximation to calculate the 
amplitudes for relativistic Coulomb excitation. It amounts to using ak = Sko on the right hand 
side of Eq. (|). The time integrals can be evaluated analytically for the VEiit) perturbations, 
given by Eqs. (p3), (PB), and (pQl). The result is 



"1st 



3 hvh 



.V2 



7 



(31) 



where Ki {K2) is the modified Bessel function of first (second) degree, and ^ = ub/'jv. For the 
E2 and Ml multipolarities, we obtain respectively. 



AE2) 
'Ist 



2i , — 



^'"^ e\K2{0 Mf,{E2,2) + Mf,{E2,-2) 



30 'jhvb'^ 

+ 27 (2-^) K,iO [Mf,{E2,-l)+Mf^{E2,l 



and 



(Ml) 
'Ist 



IStt Zie 



^K,{0 Mf,{Ml,l)-Mf,{Ml,-l] 



(32) 



(33) 



3 hcb 

These expressions are the same as those obtained from the formulae deduced in Ref. [0. 
We note that the multipole decomposition developed by those authors is accomplished by a 
different approach, i.e., using recurrence relations for the Gegenbauer polynomials, after the 
integral on time is performed. Therefore, the above results present a good check for the time- 
dependence of the multipole fields deduced here. 



The formulas above have been derived under the assumption of the long-wavelength approx- 
imation. When this approximation is not valid the matrix elements given by Eqs. (p2| , p7D are 



to be replaced by the non-approximated matrix-elements for electromagnetic excitations \lt 
i.e., 



M{M\ii) = -i 



(2A + 1)!! 

/s;^+ic(A + 1) 

(2A + 1)!! 



J(r) ■ V X L 



jx{K-r)Yx^,{r) 



JW-L 



]x{Kr)Yx^,{Y) 



(34) 



(35) 



K^C{X + 1) 

for electric and magnetic excitations (k = uj/c), respectively. However, the other factors do not 
change (see, e.g.. Ill 



2.1.2 Excitation probabilities and virtual photon numbers 

The square modulus of Eqs. (|3ll , ^|33D gives the probability of exciting the target nucleus from 



the initial state | ItMi > to the final state | IfMf > in a collision with impact parameter b. If 
the orientation of the initial state is not specified, the probability for exciting the nuclear state 
of energy Ef and spin If is 

^ Mi,Mf 

Integration of ( pGf ) over all energy transfers e = huj, and summation over all possible final states 
of the projectile nucleus (making use of the Wigner-Eckart theorem and the orthogonality of 
the properties of the Clebsch-Gordan coefficients) leads to the Coulomb excitation probability 
in a collision with impact parameter b: 

Pc = lll ^-/(^) (37) 
/ 

where pf{s) is the density of final states of the target with energy Ef = Ei + e. 
Inserting (|3l]j3|,|3|) into (g^) one finds 

Pc{b,e) = Y: P.xib,e) = T.J^ ^-a(&,^) <(^) (38) 



ttA ttA 



where 



^-^'^'^ = ApA + ^)nf T.Pfie)^''-' B{.\h - //) (39) 

are the photonuclear absorption cross sections for a given multipolarity vrA. The total photonu- 
clear cross section is a sum of all these multipolarities, 

^7 = E<(^)- (40) 



The functions nT^\{e) are called the virtual photon numbers, and are given by 

^j. jji y I 7^ 



nEi(h,e) = ?^%{y[Kl + ^Kl\ (41) 



Z{a c 4 J 4 



nE2{h, e) = ^ (-)^ I [Ki + iK,K, + eK',\ + - v'/c^Ki ) (42) 



and 



nMiib,e) 



7r2 62 



(43) 



where all K^'s are functions of ^(6) = 

Since all nuclear excitation dynamics is contained in the photoabsorption cross section, 
the virtual photon numbers |42| , |43|) do not depend on the nuclear structure. They are 
kinematical factors, depending on the orbital motion. They may be interpreted as the number 
of equivalent (virtual) photons that hit the target per unit area. These expressions show that 
Coulomb excitation probabilities are exactly directly proportional to the photonuclear cross 
sections, although the exchanged photons are off-shell. This arises from the condition that the 
reaction is peripheral and the nuclear charge distributions of each nuclei do not overlap during 
the collision. This result can be proved from first principles, and has been shown in some 



textbooks (see, eg. |2l[|). 

The usefulness of Coulomb excitation, even in first order processes, is displayed in Eq. (|38|). 
The field of a real photon contains all multipolarities with the same weight and the photonu- 
clear cross section (|40|) is a mixture of the contributions from all multipolarities, although only 
a few contribute in most processes. In the case of Coulomb excitation the total cross section is 
weighted by kinematical factors which are different for each projectile or bombarding energy. 
This allows one to disentangle the multipolarities when several ones are involved in the excita- 
tion process, except for the very high bombarding energies 7^1 for which all virtual photon 
numbers can be shown to be the same 



2.1.3 Cross sections and total virtual photon numbers 

The cross section is obtained by the impact parameter integral of the excitation probabilities. 
Eq. (|38D shows that we only need to integrate the number of virtual photons over impact 
parameter. One has to introduce a minimum impact parameter in the integration. Impact 
parameters smaller than are dominated by nuclear fragmentation processes. One finds 



da. 



(44) 



where the total virtual photon numbers Nt^x^s) = 2tt J db b n{b,e) are given analytically by 



NEi{e) 



2Zi a ,c,o 

71 V 



2c2 



Kl) 



(45) 



NeAs] 



2Zfa ,c,4 



TT V 



2(1 - "^)Kl + - "^fK,K, + ^-^{K^ - Kl 



(46) 



and 



NmAs) 



71 6^ 



(47) 



where all K^j^^s are now functions of ^(6) 



■yv 



2.2 Coulomb excitation at intermediate energies 



2.2.1 Classical trajectory: recoil and retardation corrections 

The semiclassical theory of Coulomb excitation in low energy collisions accounts for the Ruther- 



ford bending of the trajectory, but relativistic retardation effects are neglected |17|. On the 



other hand, in the theory of relativistic Coulomb excitation |T8[ recoil effects on the trajectory 
are neglected (one assumes straight-line motion) but retardation is handled correctly. In fact, 
the onset of retardation brings new important effects such as the steady increase of the exci- 
tation cross sections with bombarding energy. In a heavy ion collision around lOOA MeV the 
Lorentz factor 7 is about 1.1. Since this factor enters the excitation cross sections in many 
ways, like in the adiabaticity parameter 

m = ^ , (48) 

one expects that some sizable modifications in the theory of relativistic Coulomb excitation 
should occur ||2^. Recoil corrections are not negligible either, and the relativistic calculations 
based on the straight-line parametrization should not be completely appropriate to describe 
the excitation probabilities and cross sections. The Coulomb recoil in a single collision is of the 
order of 

ao = ^ , (49) 

which is half- distance of closest approach in a head-on collision, with mo equal to the reduced 
mass of the colliding nuclei. Although this recoil is small for intermediate energy collisions, the 
excitation probabilities are quite sensitive to it. This is important for example in the excitation 
of giant resonances because the adiabaticity parameter is of the order of one (see, Eq. (^8])). 
When ^{b) <^ 1, the excitation probabilities depends on b approximately like l/&^, while when 
^(6) becomes greater than one they decrease approximately as e~'^'^^^''^ /b"^. Therefore, when 
^ ~ 1 a slight change of b may vary appreciably the excitation probabilities. 

In the semiclassical theory of Coulomb excitation the nuclei are assumed to follow classi- 
cal trajectories and the excitation probabilities are calculated in time-dependent perturbation 
theory. At low energies one assumes Rutherford trajectories for the relative motion while at 
relativistic energies one assumes straight-line motion. In intermediate energy collisions, where 
one wants to account for recoil and retardation simultaneously, one should solve the general 
classical problem of the motion of two relativistic charged particles. But, even if radiation 



is neglected, this problem can only be solved if one particle has infinite mass [|24[. This ap- 
proximation should be sufficient if we take, e.g., the collision ^^O + ^osp]^ ^^j, system. An 
improved solution may be obtained by use of the reduced mass, as we show next, in a formalism 
developed by Aleixo and Bertulani [ P5| . 

In the classical one-body problem, one starts with the relativistic Lagrangian 

£ = -m^c' {1 - i (r^ ^ ^2^2)|'/' _ ^ ^^Q^ 

where f and <p are the radial and the angular velocity of the particle, respectively (see Fig. 
H). Using the Euler-Lagrange equations one finds three kinds of solutions, depending on the 
sign of the charges and the angular momentum in the collision. In the case of our interest, the 
appropriate solution relating the coUisional angle (p and the distance r between the nuclei is 



m 



where 



- = A\e cos{W(p) - 1] 
r 



W 



A 



^ cLo ^ 
ZiZ2e^E 



1/2 



cLn 



ZiZ2e^E 



,mocZiZ2e^ 



1/2 



(51) 

(52) 
(53) 
(54) 



E is the total bombarding energy in MeV, irio is the mass of the particle and Lq its angular 
momentum. In terms of the Lorentz factor 7 and of the impact parameter b, E = 'yrrioc'^ 
and Lq = 'jrriovb. The above solution is valid if Lq > ZiZ2e^/c. In heavy ion collisions at 
intermediate energies one has Lq ^ ^1^26^/0 for impact parameters that do not lead to strong 
interactions. It is also easy to show that, from the magnitudes of the parameters involved in 
heavy ion collisions at intermediate energies, the trajectory (51) can be very well described by 
approximating 

76^' ' " V «g 



W = 1 



A 



+ 1 



(55) 



where Oo is half the distance of closest approach in a head on collision (if the nuclei were 
pointlike and if non-relativistic kinematics were used), and e is the eccentricity parameter. In 
the approximation ( p5D e is related to the deflection angle ■(9 by e = (0-0/7) cot-i?. 

The time dependence for a particle moving along the trajectory ( |^ may be directly ob- 
tained by solving the equation of angular momentum conservation. Introducing the paramet- 
rization 

^(x) = — [e coshx + 1] (56) 



we find 



t 



Oo 

■yv 



[x + e sinh x] ■ 



(57) 



Using the scattering plane perpendicular to the Z-axis, one finds that the corresponding com- 
ponents of r may be written as 

X = a [cosh x + e] ! (58) 

(59) 
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1 sinhx , 







(60) 



where a = ao/7. This parametrization is of the same form as commonly used in the non- 
relativistic case [0, except that a^, substituted by ao/7 = 

In the limit of straight-line motion e ~ 6/a ^ 1, and the equations above reduce to the 
simple parametrization 

y = vt , X = b , and z = . (61) 

As we quoted before, the classical solution for the relative motion of two relativistic charges 
interacting electromagnetically can only be solved analytically if one of the particles has infinite 
mass. Non-relativistically the two-body problem is solvable by introduction of center of mass 
and relative motion coordinates. Then, the result is equivalent to that of a particle with 
reduced mass tUo = mprriT/imp + mr) under the action of the same potential. The particle 



with reduced mass iiio is lighter than those with mass mp and m^, and this accounts for the 
simultaneous recoil of them. An exact relativistic solution should reproduce this behavior as 
the relative motion energy is lowered. We shall use the reduced mass definition of m^ as usual 
in the parametrization of the classical trajectory of Coulomb excitation in intermediate energy 
collisions, as outlined above. In a ^^O + ^osp^ collision this is not a too serious approximation. 
For heavier systems like U+U it would be the simplest way to overcome this difficulty. But, as 
energy increases, this approximation is again unimportant since the trajectories will be straight- 
lines parametrized by an impact parameter b. A more exact result can be obtained numerically 
using the Darwin Lagrangian to determine the classical trajectory in collisions at intermediate 



energies But, the parametrization of the classical trajectory as given by Eqs. ( pq , p^ ,|60|) 
with a reduced mass particle, besides reproducing both the non-relativistic and the relativistic 
energies, gives a reasonable solution to the kind of collisions we want to study. 

2.2.2 Excitation amplitudes 

Including retardation, the amplitude for Coulomb excitation of a target from the initial state 
I i > to the final state | / > is given in first order time- dependent perturbation theory by 

1 /" 1 
a/i = y {pfiir) r) + - J fi{r) ■ A{iu, r)} £r (62) 

where p/j (J/i) is the nuclear transition density (current) and 

^{u, r) = Zie / e^-* -— dt (63) 

J — oo I r — r (^t j I 

I Jwt ^ 



A{uj, r) = ^ / V{t) e™* dt (64) 

c J —oo I r — r (^tj I 

are the retarded potentials generated by a projectile with charge Z2 following a Coulomb tra- 
jectory, and K = u/c. When the magnitude of the amplitudes (|62D is small compared to unity, 
the use of first order perturbation theory is justified. 
We now use the expansion 

piK,\r—r'\ 

^— — = 47r z/s:^j,(«:r<)r,;(r<)/i,(Kr>)n;.(i->) , (65) 

where j\ {h\) denotes the spherical Bessel (Hankel) functions (of first kind), r> (r<) refers to 
whichever of r and r' has the larger (smaller) magnitude. Assuming that the projectile does 
not penetrate the target, we use r> (r<) for the projectile (target) coordinates. At collision 
energies above the Coulomb barrier this assumption only applies for impact parameters larger 
than a certain minimum, below which the nuclei penetrate each other. 

Using the continuity equation ([Tl|) for the nuclear transition current (we changed the no- 
tation: p = pfi, J = jfi), we can show that the expansion ( [65| ) can be expressed in terms of 



spherical tensors (see, e.g., Ref. ||2T|, Vol. II) and Eq. ( |62D becomes 



(66) 



where A4(7iX,n) are the matrix elements for electromagnetic transitions, as defined in (|34 



The orbital integrals S{'7iX,fi) are given by 
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7moC A(2A- 1)!! ° 
where Lq is the angular momentum of relative motion, which is constant: 

Lo = 'yarrioV cot — 

with -d equal to the (center-of-mass) scattering angle. 
In non-relativistic collisions 

, ujr' V uur' v 
Kr = = <-<l 

c c V c 



(67) 

(68) 

(69) 



(70) 



because when the relative distance r' obeys the relations ur'/v > 1 the interaction becomes 
adiabatic. Then one uses the limiting form of hx for small values of its argument ||2^ to show 
that 

/oo 
r'-'-\t)Yx,{e'it),(l>'it)}e'-'dt (71) 
-oo 



and 
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V{r'-'-\t)Yx, [e\t),<j>'{t)]} e'-'dt 



(72) 



which are the usual orbital integrals in the non-relativistic Coulomb excitation theory with 
hyperbolic trajectories (see Eqs. (II. A. 43) of Ref. [p!7[] ). 



In the intermediate energy case the relation (^) is partially relaxed (of course, for relativistic 
energies, v ~ c, it is not valid) and one has to keep the more complex forms ( |67| , |68| ) for the 
orbital integrals. 

Using the Z-axis perpendicular to the trajectory plane, the recursion relations for the spher- 
ical Hankel functions and the identity 

dx dr 
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r = aev sinh x ■ 



(73) 



we can rewrite the orbital integrals, in terms of the parametrization (|59| , |72|) , as 
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(74) 
(75) 



and 

uja ujao 

V = — = — , 

V 7f 

and with all h\s as functions of {y/c)ri (ecoshx + !)• 
For convenience, we define 

I{EX,^^) = — S{E,,^^) 



(76) 



c 



(77) 



and we translate the path of integration by an amount in/ 2 to avoid strong oscillations of the 
integral. We find, 
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where all /i\'s are now functions of 



z = -r] {ie sinh x + 1) . 



(79) 



In the case of magnetic excitations, one may explore the fact that Lq is perpendicular to 
the scattering plane to show that 
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The magnetic orbital integrals become 
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we obtain, using the parametrization (|58|j59|j60|) , and translating the integral path by i7r/2. 
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(83) 



Generally, the most important magnetic excitation has Ml multipolarity. The orbital inte- 
grals (^,B1) can only be solved numerically. 



2.2.3 Cross sections and equivalent photon numbers 



In the high-energy hmit the classical trajectory reduces to a straight-line. One can show that 
using the approximation e = 6/a ^ 1 the orbital integrals ( fTHD and ( pSf ) can be expressed 
in terms of simple analytical functions. However it is instructive and useful to deduce the 
excitation amplitudes from the first principles again. 

The square modulus of Eq. (|66|) gives the probability of exciting the target nucleus from 
the initial state | JjMj > to the final state | IfMf > in a collision with cm. scattering angle d. 
If the orientation of the initial state is not specified, the cross section for exciting the nuclear 
state of spin If is 

ah^ 1 



da. 



4 2Ii + 



- J2 dn, (84) 



where a'^e'^dQ/A is the elastic (Rutherford) cross section. Using the Wigner-Eckart theorem and 
the orthogonality properties of the Clebsch-Gordan coefficients, one can show that 

~ln~-~h^''' ^ (2A + 1)3 1^^^^'^)! ' ^s^) 

where tt = E or M stands for the electric or magnetic multipolarity, and the reduced transition 
probability is given by 

B(n\;I,-^Ij) = — ^ J2 \<hM,\MinX, ^^)\IJMf> 

= ^^|</.||-M(^A)||/^>f . (86) 

Integration of (^) over all energy transfers e = hu, and summation over all possible final 
states of the projectile nucleus leads to 



dQ ^ J dVt 

where Pf{e) is the density of final states of the target with energy Ej = Ei + e. Inserting (^) 
into (|87|) one finds 

d(Tc daT,x ^ f de dn.„x ^ 

ttA ttA 

where o"^^ are the photonuclear absorption cross sections for a given multipolarity vrA. The 
virtual photon numbers^ riT^xi^), are given by 

dn^x Zfa A[(2A + l)!!f 

"dfT = ^ (A + 1)(2A + 1)3 ^ I I • 

In terms of the orbital integrals I {EX, /x), given by ([79|) , and using the Eq. (^9l), we find for 
the electric multipolarities 
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In the case of magnetic excitations we find 
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Only for the £^1 multipolarity the integrals can be performed analytically and we get closed 
expression 

^ = I? e-C {ij H- IK[,(K)?] , (92) 

where e = 1/ sin(6'/2), a = 1/137, ( = uao/^yv, ao = ZiZ2e'^ /2ELab, Kic^ is the modified Bessel 
function with imaginary index, K[(- is the derivative with respect to its argument. 

Since the impact parameter is related to the scattering angle by 6 = acot'i9/2, we can also 
write 

which are interpreted as the number of equivalent photons of energy e = hu, incident on the 
target per unit area, in a collision with impact parameter b, in analogy with the results obtained 
in section 2.1.2. 

Again we stress the usefulness of the concept of virtual photon numbers, especially in rela- 
tivistic collisions. In these collisions the momentum and the energy transfer due to the Coulomb 
interaction are related by Ap = AE/v ~ AE/c. This means that the virtual photons are almost 
real. One usually explores this fact to extract information about real photon processes from the 
reactions induced by relativistic charges, and vice-versa. This is the basis of the Weizsacker- 
Williams method, commonly used to calculate cross sections for Coulomb excitation, particle 
production, Bremsstrahlung, etc (see, e.g., Ref. 0). In the case of Coulomb excitation, even 
at low energies, although the equivalent photon numbers should not be interpreted as (almost) 
real ones, the cross sections can still be written as a product of them and the cross sections 
induced by real photons, as we have shown above. 



2.3 Comparison of Coulomb excitation of GR's at low energies and 
at relativistic energies 



Inserting the non- relativistic orbital integrals into Eq. (P^, we get the following relation for 
the non-relativistic equivalent photon numbers (A^-R) 

where 6 = for electric, and 6 = —1 for magnetic multipolarities, and ( = uoao/v. The non- 
relativistic Coulomb excitation functions fnxi'&X) ^i-re very well known and, e.g., are tabulated 
in Ref. [|l^], or maybe calculated numerically. 

Using the Eqs. (|90D , (|9T|) and (0), we make an analysis of Coulomb excitation extending 
from low to high energy collisions. As an example, we study the excitations induced by ^'^'^Pb 
in ^^O + 208p]^ collisions. Since the expression (|89|) is quite general, valid for all energies, 



under the assumption that the nuclei do not overlap, the equivalent photon numbers contain 
all information about the differences among the low and the high energy scattering. In Figs. 
|, I and I we show dn^x,s, for the El (Fig. D, E2 (Fig. §, and Ml (Fig. |) multipolar- 
ities, and for the collision ^^O + ^osp]^ with an impact parameter 6 = 15 fm. They are the 
equivalent photon numbers with frequency uj = 10 MeV/h incident on ^°^Pb. The dotted lines 
are obtained by using the non-relativistic Eq. (P^), while the dashed lines correspond to the 
relativistic expressions (^ , ^^31) . One observes that the relativistic expressions overestimate 
the equivalent photon numbers at low energies, while the non-relativistic expressions underes- 
timate them at high energies. The most correct values are given by the solid lines, calculated 
according to Eqs. (|90D and (pT]). They reproduce the low and the high energy limits, giving an 
improved interpolation between these limits at intermediate energies. These discrepancies are 
more apparent for the El and the E2 multipolarities. In the energy interval around lOOA MeV 
neither the low energy theory nor the high energy one can reproduce well the correct values. 
This energy interval is indeed very sensitive to the effects of retardation and of Coulomb recoil. 

At these bombarding energies, the differences between the magnitude of the non-relativistic 
and the correct relativistic virtual photon numbers are kept at a constant value, of about 20%, 
for excitation energies e = hu < 10 MeV. However, they increase sharply when one reaches the 
excitation energy of e = huj > 10 MeV. The reason is that, for such excitation energies, the 
adiabaticity factor becomes greater than unity (^ > 1). This means that excitation energies 
of order of 10 MeV (like in the case of giant resonance excitation) are in the transition region 
from a constant behavior of the equivalent photon numbers to that of an exponential (~ e~'"^) 
decay. This is more transparent in Fig. ^ where we plot the equivalent photon numbers for 
Eiab =100A MeV, 6 = 15 fm, and as a function of hu. One also observes from this figure 
that the E2 multipolarity component of the electromagnetic field dominates at low frequencies. 
Nonetheless, over the range of hu up to some tens of MeV, the E2 matrix elements of excitation 
are much smaller than the El elements for most nuclei, and the E2 effects become unimportant. 
However, such effects are relevant for the excitation of the isoscalar E2 giant resonance (GQRjs) 
which have large matrix elements. 

As an application of the semiclassical approach to Coulomb excitation in intermediate energy 
collisions, we study the excitation of giant isovector dipole resonances (El) and of giant isoscalar 
quadrupole resonances {E2) in ^osp]^ ^gans of the Coulomb interaction with a ^^O projectile. 
At lOOA MeV the maximum scattering angle which still leads to a pure Coulomb scattering 
(assuming a sharp cut-off at an impact parameter b = Rp + Rt) is 3.9°. The cross sections are 
calculated by assuming a Lorentzian shape for the photonuclear cross sections: 

= (e2 _ Ely + e2r2 ^^^^ 
with Cm chosen to reproduce the Thomas-Reiche-Kuhn sum rule for El excitations, 

NZ 

^\e) ~ 60— - MeVmb (96) 
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and the energy-weighted sum rule for the quadrupole mode, 

de 



J (e) ^ ~ 0.22 /ib/MeV . (97) 



The resonance energies are approximately given by -Egdr 

~ 77 ■ A-^/^ MeV and E, 
63 ■ MeV. We use the widths Toun = 4 MeV and Fgqr = 2.2 MeV for ^ospb. 



GQR 



We will discuss the differential cross sections as a function of the scattering angle later, 
when we introduce the effects of strong absorption. To obtain the total cross sections, one has 
to integrate the equivalent photon numbers in ( pOD and (pi]) from 0° to a maximum scattering 
angle 6max, where the nuclear absorption sets in, or equivalently, one can integrate over the 
impact parameter, from a minimum value bmim up to infinity. In Fig. ^ is shown the total 
cross section for the excitation of giant dipole and of giant quadrupole resonances in ^°^Pb in 
a collision with ^^O as a function of the laboratory energy per nucleon. The same average 
behavior of the photonuclear cross sections, as assumed in Eqs. (^) and (|96D , is used. 

Only for the El multipolarity the angular integration can be performed analytically. One 
obtains^ 

TT L ^ 

X 




{aeKl + Kl-Kl--\K,,[^] -K[,[^-^] \\ , (98) 



where 



I ^1 foi' 2ct > bmim i fr\r\\ 

' " ' - for 2a < bmim , 




and ^ = eoC = i^bmim/iv. 

It is easy to see that this equation reduces to Eq. (^51) in the relativistic limit, when C, — > 0, 
Co — > oo. 

The cross sections increase very rapidly to large values, which are already attained at inter- 
mediate energies. A salient feature is that the cross section for the excitation of giant quadrupole 
modes is very large at low and intermediate energies, decreasing in importance (about 10% of 
the El cross section) as the energy increases above lA GeV. This occurs because the equivalent 
photon number for the E2 multipolarity is much larger than that for the El multipolarity at low 
collision energies. That is, ^ ^ei? for v <^ c. This has a simple explanation. Pictorially, as 
seen from an observator at rest, when a charged particle moves at low energies the lines of force 
of its corresponding electric field are isotropic, diverging from its center in all directions. This 
means that the field carries a large amount of tidal {E2) components. On the other hand, when 
the particle moves very fast its lines of force appear contracted in the direction perpendicular 
to its motion due to Lorentz contraction. For the observator this field looks like a pulse of 
plane waves of light. But plane waves contain all multipolarities with the same weight, and 
the equivalent photon numbers become all approximately equal, i.e., uei — — nMi, and 
increase logarithmically with the energy for 7^1. The difference in the cross sections when 
7^1 are then approximately equal to the difference in the relative strength of the two giant 
resonances a^'^/a^^ < 0.1. The excitation of giant magnetic monopole resonances is of less 
importance, since for low energies umi <^ ^bi (?t-mi — {v/c)^nEi), whereas for high energies, 
where Umi — "^ei, it will be also much smaller than the excitation of electric dipole resonances 
since their relative strength a^'^^/a^^ is much smaller than unity. 

At very large energies the cross sections for the Coulomb excitation of giant resonances 
overcome the nuclear geometrical cross sections. Since these resonances decay mostly through 
particle emission or fission, this indicates that Coulomb excitation of giant resonances is a 
very important process to be considered in relativistic heavy ion collisions and fragmentation 



^We observe that the original formula for the dipole case appearing in Q has a misprinted sign in one of its 
terms. 



processes, especially in heavy ion colliders. At intermediate energies the cross sections are also 
large and this offers good possibilities to establish and study the properties of giant resonances. 



2.4 Quantum description of Coulomb excitation at high energies 

Inelastic scattering of heavy ions at intermediate energy collisions is an important tool to 
investigate the structure of stable and unstable nuclei. Laboratories like GANIL/France, 
GSl/Germany, RIKEN/ Japan, and MSU/USA, frequently use this technique. The angular 
distribution of the inelastically scattered fragments are particularly useful to identify unam- 
biguously the multipolarity of the interaction, and consequently the spin and parities of the 
excited states. In previous sections we have shown that recoil and retardation effects, are im- 



portant at this energy regime. However, as shown by Bertulani and Nathan ||2^, in order to 
describe correctly the angular distribution, absorption and diffraction effects have to be in- 
cluded properly. Next we show how quantum mechanical effects show up in the differential 
cross sections. 



2.4.1 Inelastic amplitudes and virtual photon numbers 

Defining r as the separation between the centers of mass of the two nuclei and r' to be the 
intrinsic coordinate of the target nucleus to first-order the inelastic scattering amplitude is given 
by 

f{9) = ^ Jd'r d'r' < ^^tV) ) I V^r, r') | <|.(+V) 0.(r') > , (100) 

where (r) and (r) are the incoming and outgoing distorted waves, respectively, for the 
scattering of the center of mass of the nuclei, and 0(r') is the intrinsic nuclear wavefunction of 
the target nucleus. 

At intermediate energies, AE/Eiab <^ 1, and forward angles, 6 <^ 1, we can use eikonal 
wavefunctions for the distorted waves; i.e., 

4'^*(r) 4''^(r) = exp|-zq.r + zx(&)| , (101) 

where 

X{b) = - / UT{z', b) dz' + ii,c{h) (102) 
nv J-oo 

is the eikonal-phase, q = k' — k, [/^* is the nuclear optical potential, and ipcib) is the Coulomb 
eikonal phase. We have defined the impact parameter b by b = |r x z|. 

For light nuclei, one can assume Gaussian nuclear densities, and the Coulomb phase is given 

by 

Mb) = 2 [Hkb) + Ie,{^)\ , (103) 

G 

with R^Q equal to the size parameter of each Gaussian matter density, Rq = [Rg^]"^ + [-Rg^]^) 
and 



oo 



Ei{x) = I ^ dt. (104) 



X 



The first term in Eq. (|103|) is the contribution to the Coulomb phase of a point-like charge 
distribution. It reproduces the elastic Coulomb amplitude when introduced into the eikonal 
expression for the elastic scattering amplitude. The second term in Eq. (|103|) is a correction 
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due to the extended Gaussian charge distribution. It eliminates the divergence of the Coulomb 
phase at 6 = 0, so that 
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(105) 



where C = 0.577 is the Euler constant. 

For heavy nuclei a "black-sphere" absorption model is more appropriate. Assuming an 
absorption radius Rq, the Coulomb phase is given by 
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(106) 



Again, the first term inside the parentheses is the Coulomb eikonal phase for pointlike charge 
distributions. The second term accounts for the finite extension of the charge distributions. 

For high energy collisions, the optical potential U{r) can be constructed by using the t-pp 
approximation p^. One gets 
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where aNN is the nucleon-nucleon cross section, and unn the real-to-imaginary ratio of the 
forward {6 = 0°) nucleon-nucleon scattering amplitude. A set of the experimental values of 
these quantities, useful for our purposes, is given in Table 0. 

In Eq. ( |100| ) the interaction potential, assumed to be purely Coulomb, is given by 



Vintir, r ) = — j^(r') 
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where = (c, v), with v equal to the projectile velocity, k = u/c, and j/i(r') is the charge 
four-current for the intrinsic excitation of nucleus 1 by an energy of hu. Inserting (|102|) and 



( p.08| ) in ( [Loop and following the same steps as in Ref. 0], one finds 

Zick 

'jhv , 

' TrAm 



X ^-Am(^) <If Mf\ MinX, -m) \h Mi > (109) 

where irXm denotes the multipolarity, G-j^xm are the Winther- Alder relativistic functions [Q, 
and < // Mf\Ai{7TX, — m)|Jj Mi > is the matrix element for the electromagnetic transition of 
multipolarity TrAm from |/j Mi > to |// Mf >, with Ef — Ei = huj. The function Qrn{(l) is given 

by 

= / rf6 6 Jmiqb) kJ — ) exp{zx(&)} , (110) 

Jo X'JV/ 

where q = 2k sin(^/2) is the momentum transfer, 6 and 4> are the polar and azimuthal scattering 
angles, respectively. 

For the El, E2 and Ml multipolarity, the functions GirXmic/v) are given by |]18 



GEn{x) = —Gei-i{x) = xVSn/S; Geio{x) = — i4^7r(x^ — l)/3; 
G'A/ii(a;) = Gmi-i{x) = -iy/Sn/S; GMioix) = -iA^Jn{x'^ - l)/3; 
G£;22(a;) = 6*^2-2(3;) = -2x\J-k{x'^ - l)/6; 



G'i52i(a;) = -G'i,2-i(a:) = -^2^^(2x2 - l)/5; Gi^2o(x) = 2^Ti{x^ - l)/5 . (Ill) 

Using the Wigner-Eckart theorem, one can calculate the inelastic differential cross section 
from ( |1U9D , using techniques similar to those discussed in previous sections. One obtains 

(^,)=i^i:^<MJ5,) (112) 



dVtdE^ ' ^7 

where a'!^^{E^) is the photonuclear cross section for the absorption of a real photon with energy 
E^ by nucleus 2, and driT^x/dVL is the virtual photon number, which is given by 

^^-A _ o (uk\^ x[i2X + m' 2 10 (,M2 nni 

where a = /Tic. 

The total cross section for Coulomb excitation can be obtained from Eqs. ( |112| ) and ( [113| ), 
using the approximation dVt ~ 2TTqdq/k'^, valid for small scattering angles and small energy 
losses. Using the closure relation for the Bessel functions, we obtain 

^ iE,)= ^ Y.n.x{E,) {E,) , (114) 

where 

, , ^2 A[(2A + 1)!!]' 

and 
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2-(:^) /^^^^^(^) exp{-2xK&)}, (116) 

where X-f(^) is the imaginary part of x{^)i which is obtained from Eq. (|102|) and the imaginary 
part of the optical potential. 



Before proceeding further, it is worthwhile to mention that the present calculations differ 
from those of previous sections by the proper inclusion of absorption. To reproduce the angular 
distributions of the cross sections, it is essential to include the nuclear transparency. In the 
limit of a black-disk approximation, the above formulas reproduce the results presented in Ref. 
0. One also observes that the Coulomb phase in the distorted waves, which is necessary 
for the quantitative reproduction of the experimental angular distributions, is not important 
for the total cross section in high energy collisions. This fact explains why semiclassical and 
quantum methods give the same result for the total cross section for Coulomb excitation at 
relativistic energies 0. At intermediate energies, however, it is just this important phase which 
reproduces the semiclassical limit for the scattering of large-Z ions, as we shall see next. Using 
the semiclassical terminology, for Eiat ~ lOOA MeV or less, the recoil in the Coulomb trajectory 
is relevant. At the distance of closest approach, when the Coulomb field is most effective at 
inducing the excitation, the ions are displaced farther from each other due to the Coulomb 
recoil. As we discussed before, this effect can be accounted for approximately by using the 
effective impact parameter fog// = b -\- 7iZiZ2e^ /4Eiab in the semiclassical calculations. This 
recoil approximation can also be used in Eq. replacing b by beff in the Bessel function 

and the nuclear phase, in order to obtain the total cross section. Since the modified Bessel 
function is a rapidly decreasing function of its argument, this modification leads to sizable 
modifications of the total cross section at intermediate energy collisions. 

Finally, we point out that for very light heavy ion partners, the distortion of the scattering 
wavefunctions caused by the nuclear field is not important. This distortion is manifested in 
the diffraction peaks of the angular distributions, characteristic of strong absorption processes. 
If ZiZio. ^ 1, one can neglect the diffraction peaks in the inelastic scattering cross sections 
and a purely Coulomb excitation process emerges. One can gain insight into the excitation 
mechanism by looking at how the semiclassical limit of the excitation amplitudes emerges from 
the general result ( |113| ). We do this next. 

2.4.2 Semiclassical limit of the excitation amplitudes 

If we assume that Coulomb scattering is dominant and neglect the nuclear phase in Eq. ( |102| ), 

we get 

n„,{q) ^ dbb JrrMb) kJ — ] exp{#c(&)} • (117) 
Jo K'jv/ 



This integral can be done analytically by rewriting it as 

'ujb^ 
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QUq) = / db 61+*^" Jmiqb) kJ — ), (118) 
Jo X'yv/ 



where we used the simple form ipc{b) = 21] ln{kb), with rj = ZiZ2e^ /hv. Using standard 
techniques found in Ref. pOf, we find 



ml 

X A" (^^^ F(^l + m + ir]; 1 + iri; l + m;-A'^^ , (119) 

where 

A = ^ , (120) 

to 



and F is the hypergeometric function |30 . 



The connection with the semiclassical results may be obtained by using the low momentum 
transfer limit 



Jmi^b) ^ ^— cos(g6-^-^ 



1 



^igb g-j7r{m+l/2)/2 _^ ^-igb ^i7r{m+l/2)/2 I ^^21) 



and using the stationary phase method, i.e., 

G{x) e''^(") dx ~ ( ^;7^ j <^(^o) e'-^^'o) , (122) 
where 

d(j) 11/ \ d'^4> , 

— (xo) = and (xq) = ^(^o) • (123) 

This result is valid for a slowly varying function G{x). 

Only the second term in brackets of Eq. ( |121| ) will have a positive {b = bo > 0) stationary 
point, and 

where 

0(6) = -g6 + 2ri ln{kb) . (125) 

The condition 0'(^o) = implies 

*. = ^ = ^. (126) 

g sm(fc'/2) 

where oq = ZiZ2e^ / ^v"^ is half the distance of closest approach in a classical head-on collision. 

We observe that the relation ( |126| ) is the same [with cot(6'/2) ~ sin~^(^/2)] as that between 
impact parameter and deflection angle of a particle following a classical Rutherford trajectory. 
Also, 

which implies that in the semiclassical limit 



■^''^-o) = -TJ = -ll . (127) 



P UnjRuth '"V7^sin(^/2)y • ^ ^ 

Using the above results, Eq. ( p.l3|) becomes 

- ''^) ^l- i-i'W^r^. E ^) . (129) 



da USljRutk ' Vt"/ (2ir)3(A + l) tr "■V7i'sin(e/2) 

If strong absorption is not relevant, the above formula can be used to calculate the equivalent 
photon numbers. The stationary value given by Eq. (|126|) means that the important values of 



b which contribute to flrn{(l) are those close to the classical impact parameter. Dropping the 
index from Eq. (|126|) , we can also rewrite ( p.29| ) as 



27rbdb ' (27r)3(A + l) ^ ' ""K-fv 



which is equal to the semi-classical expression given in Ref. pSj' Eq. (A. 2). 

For very forward scattering angles, such that A << 1, a further approximation can be made 
by setting the hypergeometric function in Eq. ( |119| ) equal to unity and we obtain 



fi^(g) = 22^'' — r(l + m + zr/)r(l + zr/)A'- H- . (131) 



m! ' V \ u / 

The main value of m in this case will be m = 0, for which one gets 



Qo{q) ^ 22^T(l+zr/) r(l + zr7) 

= -V 

and 



' 2"" r(i,,) r(,,) fray*'"' _ (132) 
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which, for ?7 ^ 1, results in 



|fio(g)r = [^-\ , . , (133) 



|Qo(g)r = 47rS' (^)%-^-^. (134) 



This result shows that in the absence of strong absorption and for t] ^ 1, Coulomb excitation 
is strongly suppressed at 6* = 0. This also follows from semiclassical arguments, since ^ — > 
means large impact parameters, 6^1, for which the action of the Coulomb field is weak. 

2.5 Singles spectra in Coulomb excitation of GDR 

In this section, we apply the formalism developed in previous sections in the analysis of the 
data of Ref. |]HT]], in which a projectile of ^^O with an energy of Eiab = 84A MeV excites the 
target nucleus ^osp]^ ^^le GDR. We first seek parameters of the optical potential which fits 
the elastic scattering data. We use the eikonal approximation for the elastic amplitude in the 
form given by 

Ue) =ikj Mqb) |l - exp[zx(&)]} bdb , (135) 

where Jq is the Bessel function of zeroth-order and the phase x{b) is given by Eq. ( |102| ). In Fig. 
^ we compare the calculated elastic scattering angular distribution to the data from Ref. [|T^]. 
The calculation utilized Eq. (|135|) , with x(6) obtained from an optical potential of a standard 
Woods-Saxon form with parameters 

Vo = 50 MeV , 1^0 = 58 MeV , Ry = Rw = 8.5 fm and 

ay = aw = 0.85 fm . (136) 

In order to calculate the inelastic cross section for the excitation of the GDR, we use a 
Lorentzian parameterization for the photoabsorption cross section of ^°^Pb , assumed to be 



all El, with i^^GDR = 13.5 MeV and F = 4.0 MeV. Inserting this form into Eq. (|1 14|) and doing 



the calculations implicit in Eq. (|113|) for duEi/dQ, we calculate the angular distribution and 
compare it with the data in Fig. p. The agreement with the data is excellent, provided we 
adjust the overall normalization to a value corresponding to 93 % of the energy weighted sum 
rule (EWSR) in the energy interval 7 — 18.9 MeV. Taking into account the ±10% uncertainty 
in the absolute cross sections quoted in Ref. |T2|, this is consistent with photoabsorption cross 
section in that energy range, for which approximately 110 % of the EWSR is exhausted. 

To unravel the effects of relativistic corrections, we repeat the previous calculations unplug- 
ging the factor 7 = (1 — f ^/c^)^^ which appears in the expressions (|115|) and ( |116|) and using 
the non- relativistic limit of the functions Geitu, as given in Eq. ( |111| ). These modifications 



eliminate the relativistic corrections on the interaction potential. The result of this calculation 
is shown in Fig. |10| (dotted curve). For comparison, we also show the result of a full calcula- 
tion, keeping the relativistic corrections (dashed curve). We observe that the two results have 
approximately the same pattern, except that the non-relativistic result is slightly smaller than 
the relativistic one. This fact may explain the discrepancy between the fit of Ref. [|I^] and ours 
as due to relativistic corrections not properly accounted for in the ECIS code [Q. In fact, if 
we repeat the calculation for the excitation of GDRjt, using the non-relativistic limit of Eqs. 
( [L15| ) and (|116|) , we find that the best fit to the data is obtained by exhausting 113% of the 
EWSR. This value is very close to the 110 % obtained by Barrette et al ||T2|. 



In Fig. |T0| we also show the result of a semiclassical calculation (solid curve) for the GDRj^, 
excitation in lead, using Eq. ( 129 ) for the virtual photon numbers. One observes that the 
semiclassical curve is not able to fit the experimental data. This is mainly because diffraction 
effects and strong absorption are not included. But the semiclassical calculation displays the 
region of relevance for Coulomb excitation. At small angles the scattering is dominated by large 
impact parameters, for which the Coulomb field is weak. Therefore the Coulomb excitation is 
small and the semiclassical approximation fails. It also fails in describing the large angle data 
(dark-side of the rainbow angle), since absorption is not treated properly. One sees that there 
is a "window" in the inelastic scattering data near = 2 — 3° in which the semiclassical and 
full calculations give approximately the same cross section. 

In Fig. we perform the same calculation, but for the excitation of the GDR, the isoscalar 
giant quadrupole resonance (CQRi^), and the isovector quadrupole resonance (GQRjt;), in Pb 
for the colhsion ^ospb + ^osp^ ^t 640A MeV. The sohd (dotted) [dashed-dotted] line is the 
differential cross section for the excitation of the GDR (GQRjs) [CQR^]. The dashed line is 
the result of a semiclassical calculation. Here we see that a purely semiclassical calculation, 
using Eq. (|9^) is able to reproduce the quantum results up to a maximum scattering angle 
9mi at which strong absorption sets in. This justifies the use of semiclassical calculations for 
heavy systems, even to calculate angular distributions. The cross sections increase rapidly with 
increasing scattering angle, up to an approximately constant value as the maximum Coulomb 
scattering angle is neared. This is explained as follows. Very forward angles correspond to large 
impact parameter collisions in which case lohj^v > 1 and the excitation of giant resonances in 
the nuclei is not achieved. As the impact parameter decreases, increasing the scattering angle, 
this adiabaticity condition is fulfilled and excitation occurs. 

As discussed above, the semiclassical result works for large Z nuclei and for relativistic 
energies where the approximation of Eq. ( |117] ) is justified. However, angular distributions 
are not useful at relativistic energies since the scattering is concentrated at extremely forward 
angles. The quantity of interest in this case is the total inelastic cross section. If we use a 
sharp-cutoff model for the strong absorption, so that Xi{b) = 00 for b < bmin and otherwise. 



then Eqs. ( |115| ) and ( p,16| ) yield the same result as an integration of the semiclassical expression, 
Eq. (|130|) , from bmin to oo. In fact, this result has been obtained earlier in Ref. 0. 



2.6 Excitation and photon decay of the GDR 

We now consider the excitation of the target nucleus to the giant dipole resonance and the 
subsequent photon decay of that excited nucleus, leaving the target in the ground state. Ex- 
perimentally, one detects the inelastically scattered projectile in coincidence with the decay 
photon and demands that the energy lost by the projectile is equal to the energy of the detected 
photon. To the extent that the excitation mechanism is dominated by Coulomb excitation, with 
the exchange of a single virtual photon, this reaction is very similar to the photon scattering 
reaction, except that in the present case the incident photon is virtual rather than real. In this 
section, we investigate whether the connection between these two reactions can be formalized. 

We first review the excitation mechanism. The physical situation is that of a heavy ion of 
energy E incident on a target. The projectile loses an energy AE while scattering through an 
angle 6. We have shown that, under the conditions AE/E <^ 1, the cross section for excitation 
of the target nucleus partitions into the following expression (we assume that the contribution 
of the E'l-multipolarity is dominant): 



dndK, ^^^^ k; In 



{E,)= — — ^ {E,) a, (E,), (137) 



where a^{E^) is the photonuclear cross section for the absorption of a real photon with energy 
E^ = AE by the target nucleus, and the remaining terms on the right-hand-side are collectively 
the number of virtual photons per unit energy with energy E^. This latter quantity depends 
on the kinematics of the scattered heavy ion and on the optical potential but is otherwise inde- 
pendent of the target degrees of freedom. This partitioning allows one to relate the excitation 
cross section to the photoabsorption cross section. 

Now, the usual way to write the cross section d'^acj/dfldEj for the excitation of the target 
followed by photon decay to the ground state is simply to multiply the above expression by a 
branching ratio R^, which represents the probability that the nucleus excited to an energy E^ 
will emit a photon leaving it in the ground state ||13[: 



d ac^ 1 dn. 



{E,)= — — ^ (E,) a, (E,) R, (E,). (138) 



dndE^ ' E^ dn 



Instead, we propose the following expression, in complete analogy with Eq. (|137|) : 



where a^^{E^) is the cross section for the elastic scattering of photons with energy E^. Formally, 
these expressions are equivalent in that they simply define the quantity R^. However, if one 
treats R^ literally as a branching ratio, then these expressions are equivalent only if it were 
true that the photon scattering cross section is just product of the the photoabsorption cross 
section and the branching ratio. In fact, it is well-known from the theory of photon scattering 
that the relationship between the photoabsorption cross section and the photon scattering cross 
section is more complicated In particular, it is not correct to think of photon scattering 
as a two-step process consisting of absorption, in which the target nucleus is excited to an 



intermediate state of energy E^, followed by emission, in which the emitted photon has the 
same energy E^. Since the intermediate state is not observable, one must sum over all possible 
intermediate states and not just those allowed by conservation of energy. Now, if the energy 
E^ happens to coincide with a narrow level, then that level will completely dominate in the 
sum over intermediate states. In that case, it is proper to regard the scattering as a two-step 
process in the manner described above, and the two expressions for the cross section will be 
equal. However, for E^ in the nuclear continuum region (e.g., in the region of the GDR), this 
will not be the demonstrated in the following discussion. 

We again consider the inelastic scattering of ^^O projectiles of energy Eiab = 84 MeV/nu- 
cleon from a ^^'^Pb nucleus at an angle of 2.5°. We use Eq. ( |113| ) to calculate the El virtual 
photon number and we use a Lorentzian parameterization of the GDR of ^°^Pb. We calculate 

and cr^-y according to the prescriptions of Ref . [|1^ and Ref . ||3^ , respectively; in both cases 
we neglect the statistical contribution to the photon decay. The results are compared in Fig. 
|T2|, where is is very evident that they make very different predictions for the cross section, 
especially in the wings of the GDR. 

We next use our expression to compare directly with the data of Ref. |TB|. For this pur- 



pose, we again calculate a^^ using the formalism of Ref. [34|, which relates cr^^ to the total 



photoabsorption. For the latter, we use the numerically-defined data set of Ref. |^ rather than 
a Lorentzian parameterization. The effect of the underlying compound nuclear levels (i.e., the 
statistical contribution to the photon scattering) is also included. The calculation is compared 
to the data in Figs. ^ and [1^. Fig. ^ shows the cross section for the excitation of the GDR 
without the detection of the decay photon. The agreement with the data is excellent, giving 
us confidence that our calculation of the virtual photon number as a function of E^ is correct. 
Fig. shows the cross section for the excitation-decay process as a function of E^. Although 
the qualitative trend of the data are well described, the calculation systematically overpredicts 
the cross section on the high-energy side of the GDR (solid curve). If the Thompson amplitude 
is not included in a^^, the calculation is in significantly better agreement with the data (dashed 
curve) . 



2.7 Nuclear excitation and strong absorption 

Up to this point we have only considered the Coulomb excitation of the nuclei, without account- 
ing for nuclear excitation. But, in peripheral collisions, the nuclear interaction between the ions 
can also induce excitations. This can be easily calculated in a vibrational model. The ampli- 
tude for the excitation of a vibrational mode by the nuclear interaction in relativistic heavy 
ion collisions can be obtained assuming that a residual interaction U between the projectile 
and the target exists, and that it is weak. According to the Bohr-Mottelson particle-vibrator 
coupling model, the matrix element for the transition i — > f is given by 

Vf,^'''\r) IfMj\U\m, >= ^/=^ < IfMf\Y,,\IM > Y,,ir) f/,(r) (140) 

where 6x = f3xR is the vibrational amplitude, or deformation length, R is the nuclear radius, 
and U\{r) is the transition potential. 

The deformation length S\ can be directly related to the reduced matrix elements for electro- 
magnetic transitions. Using well-known sum-rules for these matrix elements one finds a relation 
between the deformation length and the nuclear masses and sizes. For isoscalar excitations one 



gets |35 



^0^ 



2vr ^ — , 5l>o = A (2A + 1) — - (141) 

where A is the atomic number, < > is the r.m.s. radius of the nucleus, and E^ is the 
excitation energy. 

The transition potentials for nuclear excitations can be related to the optical potential in 
the elastic channel. The basic idea is that the interaction between the projectile and the target 
induces surface vibrations in the target. Only the contact region between the nuclei in grazing 
collisions is of relevance. One thus expects that the interaction potential is proportional to the 
derivatives of the optical potential in the elastic channel, which peak at the surface. This is 
discussed in details in Ref. ||3^. The transition potentials for isoscalar excitations are 
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for monopole, and 
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(142) 
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for quadrupole modes. 

For dipole isovector excitations 
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where Z (N) the charge (neutron) number. The transition potential in this case is p5 
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(145) 



where the factor x depends on the difference between the proton and the neutron matter radii 
as 



X- 



2{N -Z) Rr,- Rp 



np 



(146) 



3A 2 (Rn + Rp) Ro 

Thus, the strength of isovector excitations increases with the difference between the neutron 
and the proton matter radii. This difference is accentuated for neutron-rich nuclei and should 
be a good test for the quantity ARnp which enters the above equations. 

The time dependence of the matrix elements above can be obtained by making a Lorentz 
boost. Since the potentials U\ r{t) peak strongly at t = 0, we can safely approximate 6{t) ~ 
9{t = 0) = 7r/2 in the spherical harmonic of Eq. ( |140|) . One gets 



< IfMf\U\IiMi > 

7 ^7=^ < IfMf\Y,,\m, > Y,,{e = ^) Ux[r{t)] 



(147) 



where r(t) = y/b'^ + ^'^v'H'^. 

Using the Wigner-Eckart theorem, the matrix element of the spherical harmonics becomes 
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For high energy coUisions, the optical potential U{r) can be constructed by using the t-pp 
approximation [|28|, as given by Eq. ( |107| ). 

We are not interested here in diffraction and refraction effects in the scattering, but on 
the excitation probabilities for a given impact parameter. The strong absorption occurring in 
collisions with small impact parameters can be included. This can be done by using the eikonal 
approximation and the optical potential, given by Eq. ( |1U7D . The practical result is that the 
excitation probabilities for a given impact parameter b, including the sum of the nuclear and 
the Coulomb contributions to the excitation, are given by 



(b) = \a%{b) + af,{b)\' exp {-a^N j dz j Srp^ij') p^ij - v')} 



(149) 



where r = y/W + z^. The corresponding excitation cross sections are obtained by an integration 
of the above equation over impact parameters. 



2.8 Nucleon removal in peripheral relativistic heavy ion collisions 



In Fig. |I5| we plot the GDR excitation probability in Pb as a function of the impact parameter, 
for the system ^ospb + ^ospb at 640A MeV. We use 100% of the sum rule to calculate the B{E1)- 
value for the electromagnetic excitation of an isolated GDR state at 13.5 MeV. In the solid 
line, we consider absorption according to Eq. (|149|) . In the construction of the optical potential 
we used the g.s. densities calculated from the droplet model of Myers and Swiatecki ||36| in 
accordance with Shen et al. [0. We will call it by soft-spheres model. 

As shown in Ref. this parametrization yields the best agreement between experiment 
and theory. The dashed line does not include absorption. To simulate strong absorption at low 
impact parameters, we use bmm = 15.1 fm as a lower limit in the impact parameter integration 
of Eq. (j^). This value was chosen such as to lead to the same cross section as that obtained 
from the solid line. However, a more detailed comparison of the soft-sphere model for strong 
absorption and a simple semiclassical calculation, based on a single parameter bmm, is described 
next . 

In Fig. 0, we plot the nuclear contributions to the excitation probability as a function 
of the impact parameter. We study the excitation of the isoscalar giant monopole resonance 
(GMRi,), the GDRi^, and the GQRi,, in lead for the collision ^ospb + sosp^ 543^ MeV. 
The GMRjs in ^^'^Pb is located at 13.8 MeV. As discussed previously, isovector excitations are 
suppressed in nuclear excitation processes, due to the approximate charge independence of the 
nuclear interaction. We use the formalism of this section, with the deformation parameters 
such that 100% of the sum rule is exhausted. This corresponds to the monopole amplitude 

= 0.054. The GDRj^, and GQRj^ deformation parameters are 61 = 0.31 fm and 62 = 0.625 
fm, respectively. The GQRjt, excitation probability is much smaller than the other excitation 
probabilities and is, therefore, not shown. The nuclear excitation is peaked at the grazing 
impact parameter and is only relevant within an impact parameter range of ~ 2 fm. Comparing 



to Fig. we see that these excitation probabilities are orders of magnitude smaller than those 
for Coulomb excitation. Consequently, the corresponding cross sections are much smaller. We 
get 14.8 mb for the isoscalar GDR, 2.3 mb for the GQRjs, and 2.3 mb for the GDR,^. The 
interference between the nuclear and the Coulomb excitation is also small and can be neglected. 

Since they are high lying states above the continuum, giant resonances mostly decay by 
particle emission (mainly neutron emission in heavy nuclei). Therefore data on neutron removal 
in relativistic heavy ion collisions is an appropriate comparison between theory and experiment. 



As we have seen, above, nuclear excitation of GR's contribute very little to the cross section, 
as compared to Coulomb excitation. However, strong interactions at peripheral collisions also 
contribute to "direct" knockout (or stripping) of neutrons, and also should be considered. It has 
been observed 0, however, that neutron removal cross sections induced by strong interactions 

1/3 1/3 

scale with + A2 , while the Coulomb excitation cross sections scale with the projectile's 
charge as Z|, approximately. One can thus separate the nuclear contribution for the nucleon 
removal of a target (or projectile) by measuring the cross sections for different projectiles (or 
targets). 

In the semiclassical approach the total cross section for relativistic Coulomb excitation 
is obtained by integrating the excitation probabilities over impact parameter, starting from a 
minimum value 6mm- It is assumed that below this minimum value the interaction is exclusively 
due to the strong interaction ("sharp-cutofF' approximation). It has been found that with this 
approximation the Coulomb cross sections are very sensitive to the parameterization of the 
minimum impact parameter [R^, HO, ETI, 



One commonly used parameterization at relativistic energies is that of Benesh et al.[^], 
fitted to Glauber-type calculations and reading 

6^f„^ = 1.35 ■ (A^^ + A'J' - 0.75 ■ «3 + A'/') fm (150) 



which we refer to hereafter as " BCV" . In Ref . a detailed study has been made concerning 
the parameterization procedure of the minimum impact parameter. It was also found that the 
nuclear contribution to the neutron removal channels in peripheral collisions has a negligible 
interference with the Coulomb excitation mechanism. This is a very useful result since the 
Coulomb and nuclear part of the cross sections may be treated separately. 



The other parametrization is that of Kox et al. |Q which reproduced well measured total 



reaction cross sections of light and medium-mass systems: 

b^^: = 1.1 (Af + A]/' + 1.85^^^^-^ - 1.9 ) fm . (151) 



This parametrization has been used previously and a reasonable agreement with the 
measured data for In cross sections was found. It should be noted, however, that the Kox 
parametrization of total interaction cross sections has been derived mainly from experiments 
with light projectiles and that its application to heavy systems involves an extrapolation into 
a region where no data points are available. 

To achieve a good a comparison with experimental data on neutron-removal cross sections 
we will use the experimental photo- neutron emission cross sections from Refs. ||32|, A 



Lorentzian fit to the (7, ?T,)-data is used to parameterize the GDR in gold. The parameters are 
13.72 MeV excitation energy, a width of 4.61 MeV, and a strength of 128% of the TRK sum rule. 
The Lorentz parameters for the isoscalar (isovector) GQR are taken as 10.8 (23.0) MeV for the 
excitation energy, 2.9 (7.0) MeV for the width, and we assume 95% exhaustion of the respective 
sum rules [|]. With these parameters we calculate the excitation cross sections da{E)/dE for 
dipole- and quadrupole-excitations. The respective neutron emission cross sections are given 
by 

<yn = j ^^fn{E)dE, (152) 

where fn{E) is the probability to evaporate one neutron at excitation energy E. fn{E) is 
taken from the experimental (7, n)-data at low E and from a statistical decay calculation with 



the code HIVAP for excitation energies above 20 MeV. Since the three-neutron emission 
threshold in gold is above the energy of the GDR state, this channel is fed mainly by the two- 
phonon excitation mechanism, while the In cross section is dominated by the excitation of the 
GDR. 

We expect that the BCV parametrization of bmin should yield similar results as the soft- 
spheres calculation since it was derived in fitting the complementary process, the nuclear in- 
teraction, calculated also with Glauber theory. Fig. |l^ shows that this expectation could be 
verified: the soft-spheres calculation for In-removal from ^^''Au by ED processes (upper full 
curve) is almost indistinguishable from a sharp-cutoff calculation using (upper dotted 

curve). This remarkable agreement tells us that for practical purposes we can avoid the extra 
numerical complication connected with the use of a soft spheres model and corroborates the use 
of in sharp-cutoff calculations in earlier works |^ ^ . We also think that the soft-spheres 
calculation (and the sharp-cutoff calculation using b^^) is physically better justified than the 
Kox parametrization since the former is derived from realistic nuclear density distributions, 
whereas the latter is an extrapolation of measured total reaction cross sections into a region 
where no data points are available. 



We will return later to discuss the other data points of Fig. when we treat the problem 
of the excitation of multiphonon states. However, it is worthwhile noticing that a point in the 
above curve, for uranium targets, is not well reproduced by the theory. In fact, this has been 



observed in other experiments |J9|, and deserves a special treatment. 



2.9 Excitation by a deformed nucleus 

Either by using the soft-sphere model, or by means of a semiclassical calculation. Coulomb 
excitation by a relativistic projectile, or target, is well described theoretically if the charge 
distribution of the projectile is spherically symmetric |]SD[. 

However, there was found a discrepancy between theory and experiment with data with 
deformed projectiles, as measured by Justice et al. |^ for uranium projectiles. This problem 



was studied theoretically by Bertulani in [^, and we will briefly discuss it here. To obtain a 
qualitative insight of the effects we shall consider a prolate deformed projectile with a variable 
deformation. 

In the frame of reference of the projectile the Coulomb field at a position r with respect to 
the center-of-charge of the distribution is given by 

'^(r) = 4vr E ^,{0, 0) M{E\^^) , (153) 

where 

M{E\^^) = J p(r') r" Yx^{r') d\' , (154) 

with p(r) equal to the ground state charge distribution of the projectile. For simplicity we will 
consider a uniform spheroidal charge distribution with the z-axis along the symmetry axis. The 
charge distribution drops to zero for distances to the center greater than the angle- dependent 
radius 

R{e) = Ro{l + f3Y2o{e)) . (155) 
In lowest order in the multipole expansion, Eq. ( |153| ) becomes 



+ i y^oie) gf) (156) 



where Q'q^ is the quadrupole moment of the charge distribution, 



Q^^ = Z^e Rlp{l + 0-16/?) + 0{p^) . (157) 
V 57r 



To obtain the (time-dependent) field in the frame of reference of the target we perform a 
Lorentz transformation of Eq. (|156|) . For a straight-hne trajectory one finds 

0(r',t) = ^ + 7/f ^r,o(^')*^ (158) 
where r' = y/b'^ + 'y'^v'H''^, with b equal to the impact parameter, v the projectile velocity, and 

7 = (1 - t;2/c)-^/2. 

The first term in the above equation is the well-known Lienard-Wiechert potential of a 
relativistic charge. It gives rise to monopole-multipole excitations of the target, which we have 
discussed so far. The second term accounts for quadrupole-multipole excitations of the target 
and is a correction due to the deformation of the projectile. This field will depend on the 
orientation of the projectile with respect to its trajectory (see Fig. 0). We can separate the 
orientation angles from the angular position of the projectile (along its trajectory) with respect 
to the target by using the identity 

Y2oie') = J—Y^Y^UOA) Y^mix'.O) , (159) 

where [6, (p) denotes the orientation of the projectile symmetry axis with respect to the bom- 
barding axis and x' = ivt/r'{t). 

The dipole excitation of the target is the most relevant and we shall restrict ourselves to this 
case only 0. At a point r = (x, z) from the center of mass of the target the field is obtained 
by replacing r' = {b, 0,'-yvt) by [b — x,y,'-y{vt — z)] in Eq. ( |158|) . The excitation amplitude to 
first order is given by Eq. (|62|). Using the continuity equation and expanding (|62D to lowest 
order in r we find 

af, = af^+af^ (160) 

where 

a 



and 



- %\ ^ {^2(0^}? + ^ ^ E^2^(^'</') >^2™(|,0) (162) 



where ^ = ujbj^v and Ki is the modified Bessel function of order i. To simplify the notation, 
we have used the Cartesian definition of the matrix elements. The dipole matrix elements for 
the nuclear excitation are given by 

Dft^ =<j\x{z)\i> . (163) 

In terms of spherical coordinates, Df^^ = Mf^{El, 1) + Mfi{El, -1) and Df^ = Mfi{El,0). 
Thus, Eq. ([iBlD is equal to Eq. (|3l]). 

In the expression (|162D we have used the approximation l2m(X) 0) — ^2m(f , 0) which is valid 
for high energy collisions since the quadrupole field is strongly peaked at t = 0, corresponding 
to the distance of closet approach. Eqs. (|161|) and (|162|) allow us to calculate the dipole 



excitation cross section by integrating their absolute squares over impact parameter, starting 
from a minimum impact parameter for which the strong interaction sets in. Neglecting the 
diffuseness of the matter distribution of the nuclei we can write (see Fig. |1^) 

h^,n{e) ~ i?i + i?2 [l + /? l20(^ + e)] (164) 

1 /3 

with the nuclear radii given by Ri = 1.2 ■ A/. The total cross section is 

2n f dbb < |a/i(6,^])p > (165) 



a 



where the < ■ ■ ■ > sign means that an average over all the possible orientations of the projectile, 
i.e., over all angles Q = {9, 0), is done. 

We will apply the above formalism to the Coulomb excitation of ^osp]^ 238jj projectiles. 
We will give the ^ss-q artificial deformation in the range /5 = — 1 to check the dependence 
of the cross sections with this parameter. The cross section given above contains three terms: 
a = 0"! + (72 + (Tu- 0"! is due to the monopole-dipole excitation amplitude, a2 is due to the 
quadrupole-dipole excitation amplitude, and ai2 is the interference between them. 

In Fig. we present the results for the numerical calculation of the quantity 

(3=0 

A = 100 X (166) 

which is the percent correction of dipole excitations in ^osp]^ ^ uranium projectile due the 
average over the orientation of the projectile. (Ti~° is the cross section for P = 0. We present 
results for three bombarding energies, lOA GeV, lA GeV and lOOA MeV, and as a function of 

The quantity defined by Eq. ( |166| ) is independent of the nature of the state excited, since 
the dipole matrix elements cancel out. They depend on the energy of the state. In order to see 
how the effect depends qualitatively on the energy of the state we used three different excitation 
energies Eji = 1, 10 and 25 MeV, respectively. These correspond to the dotted, dashed and 
solid lines in Fig. |19|, respectively. 

One observes from Fig. ^ that the deformation effect accounted for by an average of the 
minimum impact parameter which enters Eq. ( |165[ ) increases the magnitude of the cross section. 
Thus the average is equivalent to a smaller "effective" impact parameter, since the cross sections 
increase with decreasing values of &mm- The effect is larger the greater the excitation energy is. 
This effect also decreases with the bombarding energy. For very high bombarding energies it 
is very small even for the largest deformation. These results can be explained as follows. The 
Coulomb excitation cross section at very high bombarding energies, or very small excitation 



energies, is proportional to In 



an average of \n{bmin) due to the additivity 



. Averaging over orientation of the projectile means 
aw of the logarithm. One can easily do this average 
and the net result is a rescaling of bmin as fbmin, with f smaller, but very close to one. 

For high excitation energies, or small bombarding energies, the cross section is proportional 
to exp|— 2co'6mm(^)/7'i^} due to the adiabaticity condition [0. Thus, in these situations, the 
cross section is strongly dependent on the average over orientation due to the strong variation 
of the exponential function with the argument. 

Now we consider the effect of the second term of Eq. ( [158| ), namely of the quadrupole- 
dipole excitations. In Fig. |2^ we show the excitation of a giant resonance dipole state in lead 
{Efi = 13.5 MeV) due to the second term Eq. (|158|) , as a function the deformation parameter 
f3 and for a bombarding energy of lOOA MeV. We assume that the giant dipole state exhausts 



Table 2: Cross sections (in mb) for Coulomb excitation of the giant dipole resonance in ^''^Pb 
by 238-[j projectiles at lOOA MeV. In the second (third) column the cross sections are due to 
the monopole (quadrupole)- dipole interaction. The last column is the total cross section. An 
average over the orientation of the projectile was done. A realistic value of the deformation of 
^^^U corresponds to /3 ~ 0.3. But, a variation of (3 is used to obtain an insight of the magnitude 
of the effect. 
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fully the TRK sum-rule, Eq. (|96|) , in lead. Now the average over orientation also includes the 
dependence of the quadrupole- dipole interaction on i7 = (6*, 0). As expected the cross section 
increases with /3. But it is small as compared to the monopole- dipole excitations even for a 
large deformation. At this beam energy the monopole-dipole excitation is of order of 1 barn. 

The total cross section contains an interference between the amplitudes a~^} and a~Jl . This is 
shown in Table ^ for lOOA MeV for which the effect is larger. The second column gives the cross 
sections for monopole-dipole excitations of a giant resonance dipole state in lead. The effect 
of the orientation average can be seen as an increase of the cross section as compared to the 
value in the first row (zero deformation). For /3 = 0.3 which is approximately the deformation 
parameter for ^ss-q correction to the cross section is negligible. In the third column the 
cross section for quadrupole- dipole excitation are given. They are also much smaller than those 
for the monopole-dipole excitations. The total cross sections, given in the last column, are also 
little dependent on the effect of the deformation. For /? = 0.3 it corresponds to an increase of 
3% of the value of the original cross section (first row ). This effect also decreases with the 
bombarding energy. For lA GeV, cr^^° = 5922 mb, while a = 5932 mb for (3 = 0.3, with all 
effects included. 

In conclusion, the effect of excitation by a deformed projectile, which can be studied by 
averaging over the projectile orientation, is to increase slightly the cross sections. The inclusion 
of the quadrupole-dipole interaction increases the cross section, too. However, these corrections 
are small for realistic deformations. They cannot be responsible for the large deviations of the 
experimental values of the Coulomb fragmentation cross sections from the standard theory 
[^, |18[, as has been observed [^, for deformed projectiles. 



3 Heavy ion excitation of multiphonon resonances 



3.1 Introduction 

Much of the interest on multiphonon resonances rehes on the possibihty of looking at exotic 
particle decay of these states. For example, in Ref. a hydrodynamical model was used 
to predict the proton and neutron dynamical densities in a multiphonon state of a nucleus. 
Large proton and neutron excesses at the surface are developed in a multiphonon state. Thus, 
the emission of exotic clusters from the decay of these states are a natural possibility. A more 
classical point of view is that the Lorentz contracted Coulomb field in a peripheral relativistic 
heavy ion collision acts as a hammer on the protons of the nuclei [^]. This (collective) motion 
of the protons seem only to be probed in relativistic Coulomb excitation. It is not well known 
how this classical view can be related to microscopic properties of the nuclei in a multiphonon 
state. 

Since there is more energy deposit in the nuclei, other decay channels are open for the 
multiphonon states. Generally, the GR's in heavy nuclei decay by neutron emission. One 
expects that the double, or triple, GDR decays mainly in the 2n and 3n decay channel. In fact, 
such a picture has been adopted by |3^, ^ with success to explain the total cross sections for 
the neutron removal in peripheral collisions. The method is the same that we used to explain 
the one-neutron removal cross sections, i.e., by replacing /„ by f2n, and /3„, in Eq. ( |152|) . 

Although the perspectives for an experimental evidence of the DGDR via relativistic Cou- 
lomb excitation were good, on the basis of the large theoretical cross sections, it was first found 
in pion scattering at the Los Alamos Pion Facility [Q. In pion scattering off nuclei the DGDR 
can be described as a two-step mechanism induced by the pion-nucleus interaction. Using the 
Axel-Brink hypotheses, the cross sections for the excitation of the DGDR with pions were shown 
to be well within the experimental possibilities |^3[. Only about 5 years later, the first Coulomb 
excitation experiments for the excitation of the DGDR were performed at the GSI facility in 



Darmstadt /Germany [^, In Fig. ^ we show the result of one of these experiments, which 
looked for the neutron decay channels of giant resonances excited in relativistic projectiles. 
The excitation spectrum of relativistic ^^^Xe projectiles incident on Pb are compared with the 
spectrum obtained in C targets. A comparison of the two spectra immediately proofs that 
nuclear contribution to the excitation is very small. Another experiment ||3^ dealt with the 
photon decay of the double giant resonance. A clear bump in the spectra of coincident photon 
pairs was observed around the energy of two times the GDR centroid energy in ^°'^Pb targets 
excited with relativistic ^°^Bi projectiles. 

The advantages of relativistic Coulomb excitation of heavy ions over other probes (pions. 



nuclear excitation, etc.) was clearly demonstrated in several GSI experiments ^ 



A collection of the experimental data on the energy and width of the DGDR is shown in 



Fig. ^ The data points are from a compilation from pion (open symbols), and Coulomb 
excitation and nuclear excitation (full symbols) experiments 0. 

The dashed lines are guide to the eyes. We see from Fig. |2^(a) that the energy of the DGDR 
agrees reasonably with the expected harmonic prediction that the energy should be about twice 
the energy of the GDR, although small departures from this prediction are seen, especially in 
pion and nuclear excitation experiments. The width of the DGDR seems to agree with an 
average value of -\/2 times that of the GDR, although a factor 2 seems also to be possible, as 
we see from Fig. ^(b). Fig. ^(c) shows the ratio between the experimentally determined 
cross sections and the calculated ones. Here is where the data appear to be more dispersed. 



The largest values of Ccxp/o"^;, come from pion experiments, yielding up to a value of 5 for this 
quantity. 

We now discuss many features of the double GDR excitation theoretically and some attempts 
to solve the discrepancies between theory and experiment observed in Fig. 



3.2 Perturbation theory and harmonic models 
3.2.1 Sum rules for single and double resonances 

The simplest way to determine the matrix elements of excitation of giant resonances is by 
means of sum rules under the assumption that those sum rules are exhausted by collective 
states. We have done this when we used the sum rules ( p6| , p7D . Let us look at these with more 
details, since they will be useful for the determination of the matrix elements for multiphonon 
excitations. The conventional sum rules for the dipole and quadrupole transitions, derived 
without exchange and velocity-dependent corrections, are {h = 1) 



/ 



3 1 iVZ 2 



isoscalar excitations, 



Tujf.Qr ='^^-^(?^{ ',\ ^^^'^^ , (168) 

^ ■' 2mN 47r [ NZ/A, isovector excitations. 

where D^""^ = M{Elm) and Q(™) = M{E2m). 

We explain our procedure on the example of the dipole sum rule ( [L67| ). The right hand 
side Sd of (|167|) being calculated for the fixed initial state |i > in fact does not depend on 



the choice of \i >. (This dependence is rather weak even if the exchange terms are taken into 
account). Since Sd does not depend on the projection m of the dipole operator D^^^ as well, 
it is convenient to introduce in usual way the reduced matrix elements of multipole operators, 

< f-IfMf\0'r\t-.hM, >=< IfMflIMm > {f;If\\Oi\\t;I,), (169) 

where / stands now for all quantum numbers except angular momentum ones, / and A/, and 
to perform the additional summation of Eq. ( 167 ) over m. In such a way one obtains 



^ ujj,{2If + 1)(/; If\\D\\z; k? = 3(2/, + 1)^^. (170) 

Now let us take the ground state |0 > of an even-even nucleus with angular momentum 
Jo = as an initial one \i\ It >. If we assume that the single GDR |1 >= |1; 1 > is an isolated 
state saturating the corresponding sum rule, we just divide the right hand side of ( |17(]| ) by the 
excitation energy ujiq to obtain the reduced matrix element 

(111^110)2 = ^. (171) 

In order to be able to calculate the cross section of excitation of the double GDR, we have 



to take the single GDR state |1 > as an initial one. The corresponding sum in Eq. (|17CI|) 



according to our assumption, is saturated by (i) "down" transition to the ground state |0 >, 
which has negative transition energy — tuio and, due to the symmetry properties of the Clebsch- 



Gordan coefficients, the strength which is 3 times larger than that of Eq. (|171|), and (ii) "up 



transitions to the double GDR states |2;/2 = L > where L can be equal to and 2. The 
resulting sum rule for the up transitions is 

J2 (2L + (2;L||D||1)2 = 12^fl. (172) 

L=0,2 

where u!2i^ = E2-l — Ei is the energy of the second excitation. Actually, considering, instead 
of the sum over m, the original dipole sum rule ( [I67| ) for fixed m, one can separate the two 
contributions to the sum (|172|) and find 

(2;L||Z}||1)2 = 2%. (173) 

^21 



Obviously, it is consistent with the sum rule ( p.72| ). 

Eqs. (|171|) and ( p.73| ) imply the relation between the strengths of sequential excitation 



processes, 

{2;L\\D\\l)' = 2^{l\\D\\0f. (174) 

For the equidistant vibrational spectrum this result is nothing but the standard Bose factor 
of stimulated radiation; our result is valid under more broad assumptions. The resulting en- 
hancement factor includes, in addition, the ratio of transition frequencies which, according to 
the data, is slightly larger than 1. The generalization for the third and higher order excitation 
processes is straightforward. 



3.2.2 Spreading widths of single and double resonances 

The above assumption of saturation certainly does not account for the fact that the resonances 
are wide. In fact, this might be also relevant for the calculation of total cross sections since 
the Coulomb excitation amplitudes given by may vary strongly with the excitation energy. 
Therefore they might be sensitive to the shape of the strength function. The widths of the 
resonances can be taken into account in a simplified approach, as we describe next. 

In a microscopic approach, the GDR is described by a coherent superposition of one-particle 
one-hole states. One of the many such states is pushed up by the residual interaction to 
the experimentally observed position of the GDR. This state carries practically all the El 
strength. This situation is simply realized in a model with a separable residual interaction. We 
write the GDR state as (one phonon with angular momentum IM) |1, IM) = A\j^,j |0) where 
A\^f is a proper superposition of particle-hole creation operators. Applying the quasi-boson 
approximation we can use the boson commutation relations and construct the multiphonon 
states (N-phonon states). A N-phonon state will be a coherent superposition of N-particle 
N-hole states. The width of the GDR in heavy nuclei is essentially due to the spreading width, 
i.e., to the coupling to more complex quasibound configurations. The escape width plays only 
a minor role. We are not interested in a detailed microscopic description of these states here. 
We use a simple model for the strength function We couple a state \a) (i.e. a GDR state) 
by some mechanism to more complex states \a), for simplicity we assume a constant coupling 
matrix element Vaa = {(^iVla) = {a\V\a) = v. With an equal spacing of D of the levels \a) 
one obtains a width ^ 

r = 2vr^ , (175) 

for the state |a) . We assume the same mechanism to be responsible for the width of the N- 
phonon state: one of the N- independent phonons decays into the more complex states la) 



while the other (N-l)-phonons remain spectators. We write the couphng interaction in terms 
of creation (destruction) operators cj, (cq,) of the complex states |a) as 



V 



V ( A\j., Ca + Aim 4 



(176) 



For the coupling matrix elements v^v, which connects an N-phonon state |A^) to the state 
— (N-1 spectator phonons) one obtains 

yn = {N -l,a\V\N) {N -I \Aim\ N) = y . \^ , (177) 

i.e., one obtains for the width F^r of the N-phonon state 



2'kN- 
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NT 



(178) 



where F is given by Eq. ( |175| ). 

Thus, the factor N in (|178|) arises naturally from the bosonic character of the collective 
states. For the DGDR this would mean F2 = 2Fi. The data points shown in Fig. ^(b) seem 
to favor a lower multiplicative factor. 

We can also give a qualitative explanation for a smaller F2/F1 value. We again assume 
that the damping of the collective modes is mostly due to the coupling to the background of 
complex configurations in the vicinity of the resonance energy. Then the resonance state |A > 
gets fragmented acquiring the spreading width F^. The stationary final states |/ > in the 
region of the GR are superpositions (with the same exact quantum numbers as the collective 
mode) of the form 

\f>=d/^ |A>+EcyV> , (179) 



where |A > is a pure GR state and {u > are complex many particle-many hole states. If the 
resonance component dominates in the excitation process as it should be for the one-body 
multipole operator, we find the first order amplitude a^^^ of the excitation of the individual 
state 1/ > in the fragmentation region 



C 



(180) 



Here a\''^ stands for the original first order excitation amplitude. As a function of the transition 



energy, the probability for the one-phonon excitation is 



p. 
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where we introduced the strength function J-'\{u). 

The traditional derivation of the strength function (see Ref. ||5^) is based on the rough 
assumptions concerning mixing matrix elements and the equidistant spectrum of complex states. 
The matrix elements Vxi, which couple the collective mode to the background states are assumed 
to be of the same average magnitude for all remote states |z/ > from both sides of the resonance. 
Under those conditions the resulting strength function has the Breit-Wigner (BW) shape 

J_ Ta 

27r 



^a(^) 



+ TI/4 



(182) 



where F^ is the spreading width of the collective resonance, 

'Xiy >f 
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(183) 



d is the mean level spacing of complex states, coupling matrix elements are averaged over the 
states \i> > and ux is the energy centroid. We will use in our numerical calculations the BW 
strength function (|182|) with the empirical parameters Ux and Tx- However, the same procedure 
can be applied to any specific form of J^x{^)- Later we come back to the question of justification 
of the model leading to Eqs. ( |I8^ ) and ( [1831) . 

The multiphonon states could also be reached by a direct excitation. Quite similarly, we 
can repeat the above arguments to calculate the probability for the direct excitation of a 
multiphonon state, with the appropriate spreading width and energy centroid of that state. 
The direct (or first-order) probabilities are then given by 

P^^\u) = \al^\u;)\\ (184) 

Let us now treat the case of the two-step excitation of GR (double-phonon). For simplicity, 
we denote the single-phonon state by |1 > and the double-phonon state by |2 >, the correspond- 
ing centroids being at ui and UJ2 respectively. The total probability to excite the double-phonon 
state is obtained by 

f 

= P^'\u;) + P2«rf(^) + p*"*(cj) , (185) 

where P^'** is the direct (or first-order) excitation of the double-phonon state, P^nd ^j^g 
two-step (or second-order) excitation term, and the last term in Eq. ( |185|) is the interference 
between the two. 



3.2.3 Second-order perturbation theory 

To second-order, the amplitude for a two-step excitation to a state |2 > via intermediate states 
|1 > is given by 



00 



^2o-E71^/ dt e^--^' V2,{t) dt'e^-»*yio(t'), (186) 
^ [tn)'^ J-00 J-00 

where V2i(t) is a short notation for the interaction potential inside brackets of the integral of 

Eq. ( |186D for the transition |1 >^ |2 >. 

Using the integral representation of the step function 

1 roD p-iq(t-t') c 1 if f > f' 

e{t-t') = -\im — -dq= ' / , ' (187) 

^ ' s^o+ 27Ti J-00 q + i5 1 0, if t < t\ ^ ' 



one finds 



^ 1 



dq 

-00 q 

where V stands for the principal value of the integral. For numerical evaluation it is more 

appropriate to rewrite the principal value integral in Eq. ( p.88| ) as 

dq 

-00 q 
dq 

1 

(189) 



Zn ^ J-00 q 

3rincipal value of the integ 
3 principal value integral in 

V r ^4f(u;2i-g)aJ^*(u;io + g) = 
J —00 

a\Wu2i - q) a}o*(^io + q) - ^21* (^21 + q) aio*(^io - q) 



To calculate a^'^^{uj) for negative values of u, we note that the interaction potential can be 

r n * 

written as a sum of an even and an odd part. This implies that a^'^*(— cu) = — a^^*{uj) . 

For three-phonon excitation we use the third term of the time-dependent perturbation 
expansion, and the same procedure as above (Eqs. ( |187| - 1189| )). 



3.2.4 Harmonic vibrator model 

A simplified model, often used in connection with multiphonon excitations, is the harmonic 
vibrator model. In this model, the resonance widths are neglected and the Coupled-Channel 
equations can be solved exactly, in terms of the first-order excitation amplitudes [0. The 
excitation amplitude of the n-th harmonic oscillator state, for any time t, is given by 

aS.W = i^^SL exp{ - |ai.,(t)|V2} , (190) 



where aisj(t) is the excitation amplitude for the (g.s.) — > 1 {one phonon) calculated with 
the first-order perturbation theory. 

For the excitation of giant resonances, n can be identified with the state corresponding to 
a multiple n of the single giant resonance state. This procedure has been often used in order 
to calculate the cross sections for the excitation of multiphonon giant resonances. Since this 
result is exact in the harmonic vibrator model, it accounts for all coupling between the states. 
However, this result can be applied to studies of giant resonance excitation only if the same 
class of multipole states is involved. I.e., if one considers only electric dipole excitations, and 
use the harmonic oscillator model, one can calculate the excitation probabilities, and cross 



sections, of the GDR, double-GDR, triple-GDR, etc. Eq. ( p.90| ) is not valid if the excitation 
of other multipolarities are involved, e.g., if the excitation of dipole states and quadrupole 
states are treated simultaneously. In Ref. JSD| a hybrid harmonic oscillator model has been 
used. In this work, it is assumed that the difference between the amplitudes obtained with the 
harmonic oscillator model and with n-th order perturbation theory is due to the appearance 
of the exponential term on the r.h.s. of Eq. (|190[) . This exponential takes care of the decrease 
in the occupation amplitude of the ground state as a function of time. As argued in Ref. 0], 
the presence of other multipole states, e.g., of quadrupole states, together with dipole states, 
may be accounted for by adding the first order excitation amplitudes for the quadrupole states 
to the exponent in Eq. (|190|) . This would correct for the flux from the ground state to the 
quadrupole states. In other words, Eq. ( |190| ) should be corrected to read 



■in 

a^i(7rA, t) = L ' J- exp{ - ^ \autin'\', t)\ /2} . (191) 

The harmonic oscillator model is not in complete agreement with the experimental findings. 
The double-GDR and double-GQR states do not have exactly twice the energy of the respective 
GDR and GQR states [0, p, 0. Apparently, the matrix elements for the transition from the 
GDR (GQR) to the double-GDR (double-GQR) state does not follow the boson-rule [0. This 
is borne out by the discrepancy between the experimental cross sections for the excitation of 
the double-GDR and the double-GQR with the perturbation theory, and with the harmonic 
oscillator model 0, H, P]. Thus, a Coupled-Channels calculation is useful to determine which 
matrix elements for the transitions among the giant resonance states reproduce the experimental 
data. 



Assuming that one has a'!^^{E) somehow (either from experiments, or from theory), a simple 
harmonic model following the discussion above can be formulated to obtain the include the 
widths of the states. As we have mentioned, in the harmonic oscillator model the inclusion 
of the coupling between all multiphonon states can be performed analytically [|1^]. One of 
the basic changes is that the excitation probabilities calculated to first-order, P^^{E,b), are 
modified to include the flux of probability to the other states. That is, 

P^xiE, b) = Pl^E, h) exp {-Pltih)] , (192) 

where P^l^{b) is the integral of over the excitation energy E. In general, the probability to 
reach a multiphonon state with the energy i?'-"^ from the ground state, with energy IS 
obtained by an integral over all intermediate energies 

Pi?l.(E("\ b) = ^ exp {-P^i\b)} J dE^""-^^ dE^""-^^ ... dE^^^ (193) 



3.2.5 Comparison with experiments 

The reactions ^^^Xe + ^"^Pb at 0.69A GeV and ^^^Bi + ^ospb at lA GeV have been measured 
at GSI ||3^, 1^]. We apply the formalism developed in the preceding sections to calculate the 
excitation probabilities and cross sections for these systems. 

Cross sections (in mb) for the Coulomb excitation of the GDRj^, GQRjs and GQRj„ in ^^^Xe 
incident on Pb at 0.69A GeV are given in Table ||. We have assumed that the GDRj^, GQRjs 
and the GQRjt, are located at 15.3, 12.3 and 24 MeV, and that they exhaust 100%, 70% and 
80% of the corresponding sum rules, respectively [Q. We used hmin = 1.2 ■ (^i''^ + A^2^) fm 
= 13.3 fm as a lower hmit guess and bmin = 15.6 fm suggested by the parameterization as 
an upper limit (number inside parentheses). The parameterization yields an intermediate 
value for this quantity. The contributions to various angular momentum projections of each 
state are shown separately. In the last column the total cross sections are calculated with 
the widths of the states taken into account. We use for the GDRj^,, GQRjs and GQRj^ the 
BW strength functions (|182|) with the resonance widths F = 4.8, 4 and 7 MeV, respectively 



56| . We see that states with higher angular momentum projections are more populated. The 
inclusion of the widths of the resonances in the calculation increases the cross sections by about 
10-20%. The experimental value 1110 ± 80 mb for the GDR is much smaller which made 
the authors of [^] to claim that the GDR absorbs only 65% of the sum rule (this number 
apparently contradicts to the systematics of data for real monochromatic photons |Q). Using 
this value, our result reduces to 1613 (1183) mb which seems to prefer the upper value of 6mm- 
The numbers in parentheses are also in rough agreement with the data for the GQRjs and 
GQRto. 

Using the formalism developed in Sect. 2.8 we have also calculated the cross sections for 
the nuclear excitation of the GQRjs in the same reaction. The cross sections for the excitation 
of isovector modes are reduced by a factor (A^ — Z)/A since the isovector mode is excited 
due to the difference in strength of the nuclear interaction between the target and the protons 
and neutrons of the projectile [^. This implies that the isovector excitations are strongly 
suppressed in nuclear excitations. Therefore, we do not consider them here. For the excitation 
of the GQRjs we find = 5.3 mb, if we use the deformation parameter PR = 0.7 fm for ^^^Xe. 
In the calculation of the nuclear potential we used Fermi density distributions with parameters 



Table 3: Cross sections (in mb) for the Coulomb excitation of the GDRjt,, GQRjs and GQRi„ 
in ^^^Xe incident on ^osp]^ 0.69A GeV. The cross sections in the last column are calculated 
with the widths of the states taken into account. The values outside (inside) parentheses use 
6mm=13.3 (15.6) fm. 





m = ±2 


m = ±1 


771 = 


0" total 


^width 


GDRj^ 




949 (712) 


264 (201) 


2162 (1630) 


2482 (1820) 


GQRjs 


90 (64) 


8.4 (6.09) 


14.3 (10.6) 


211 (150) 


241 (169) 


GQRto 


29.7 (25.6) 


6.1 (5.46) 


14 (12.4) 


84.1 (74.5) 


102 (93) 



Table 4: Excitation cross sections (in mb) of the GDR^, and of the [GDR]"- states in the 
reaction ^osp]^ _|_ 208p|^ 640A MeV. A comparison with first order perturbation theory and 
the harmonic oscillator is made. 



State 


1st pert. th. 


harm. osc. 


c.c. 


GDRjt, 


3891 


3235 


3210 


[GDR„]2 


388 


281 


280 


[GDR,,]3 


39.2 


27.3 


32.7 


[GDRj^,]"^ 


4.2 


2.4 


3.2 



Po = 0.17 fm~^ and R =5.6 (6.5) fm, a =0.65 (0.65) fm for Xe (Pb). The nucleon-nucleon cross 
section used was 40 mb. Again we see that the nuclear contribution to the total cross section 
is very small. 

The double dipole phonon state can couple to total angular momentum or 2. As we 
mentioned in Sect. 2, for the state with L = 2 there is the possibility of a direct quadrupole 
Coulomb excitation (L = states cannot be Coulomb excited [§). For simplicity, we do not 
consider here the physics of the isospin coupling of the two GDR. 

We calculated the direct and the two-step probabilities for the excitation of the double- 
phonon state according to the approach discussed in the previous sections. The total cross 
sections obtained are shown in Table ^ We found that the principal value term in Eq. (|n 



contributes very little (less than 1%) to the GDR x GDR cross section via a two-step process. 

From Table ^ we see that the inclusion of the widths of the final (GDR x GDR) and 
the intermediate (GDR) state increase the cross sections by 10-20%. For the position and 



width of the GDR x GDR state we took E = 28.3 MeV and F = 7 MeV, respectively [40 



which corresponds to ujiq = 15.3 MeV and u)2i = 13 MeV, both for L = 2 and L = 0. 
For the calculation of the direct excitation we assumed that the resonance would exhaust 
20% of the GQRjs sum rule. It is based on the hypotheses that the missing strength of the 
low-lying GQRjs could be located at the double dipole phonon state as a consequence of the 
anharmonic phonon coupling of the (Q-D-D) -type. Obviously, it should be considered as highly 
overestimated upper boundary of the direct excitation. In Ref. the reduced transition 
probability for the excitation of double-phonon states within the Quasiparticle-Phonon Model 
have been calculated. The value B{2~^, E2) = 4.2 fm'* has been obtained. Using this value 
we get that the cross section for the direct excitation of the L = 2 state is 12 /ib, much smaller 
than what we quote above. We conclude that even in the more optimistic cases the contribution 



Table 5: Cross sections (in mb) for the Coulomb excitation of the double GDR in ^^^Xe incident 
on Pb at 0.69A GeV. The cross sections in the last column are calculated with the widths of 
the states taken into account. The values outside (inside) parentheses use 6min=13.3 (15.6) fm. 



DGDR state 


m = ±2 


m = ±1 


171 = 


C total 


^width 


L=0 (two-step) 






22.8 (10.7) 


22.8 (10.7) 


28.4 (13.3) 


L=2 (two-step) 


23.3 (11.2) 


13.4 (6.6) 


51.4 (26.8) 


124.8 (62.4) 


154 (77) 


L=2 (direct 












- 20% of SR) 


3.27 (2.85) 


0.86 (0.77) 


2.12 (1.88) 


10.3 (9.12) 


11.8 (10.8) 



of the direct mechanism to the total cross section for Coulomb excitation of the double-phonon 
state is much less than that of the two-step process. 

Another conclusion drawn from the numbers of Table ^ is that the excitation of the L = 2 
double-phonon state is much stronger than for the L = state. Adding the two contributions 
we find that the total cross section for the excitation of the double-phonon state (excluding 
the direct mechanism) in the reaction above is equal to 182 (101) mb. The experimental value 
of Ref. 1^ is about 215 ± 50 mb. As stated above, the nuclear contribution to the (direct) 
excitation of the double-phonon state is not relevant. If we assume again that about 20% of 
the sum rule strength is exhausted by this state (using e.g. /?R=0.1 fm), we get 1.1 mb for the 
nuclear excitation of the L = 2 double-phonon state. Contrary to the single phonon case, the 
appropriate value of bmin for the double GDR experiment is 6mm=13.3 fm. 

We also compare our results with the experiment of Ritman et al. [^. They measured the 
excitation of a ^'^'^Pb target by means of 209Bi projectiles at lA GeV and obtained 770 ± 220 
mb for the excitation cross section of the double resonance. We calculate the cross sections for 
the same system, using Ei = 13.5 MeV, Fi = 4 MeV, E2 = 27 MeV and = 6 MeV for the 
energy position and widths of the GDR and the GDR x GDR in ^°^Pb, respectively. Using 
the formalism developed in Sects. 3.2.2 and 3.2.3 and including the effects of the widths of the 
states, we find ai = 5234 b for the excitation of the GDR and a2 = 692 mb for the excitation 
of the GDR x GDR, using = 1. 2{A]^^ + A]^^) fm = 14.2 fm. Thus, while the cross section 
for the excitation of single phonons is a factor 2.8 larger than that of the experiment of Ref. 
5D[, the cross sections for the excitation of double phonons is larger by a factor 3.8. This is 



due to the larger value for the excitation probabilities caused by a larger B{E1) value for this 
reaction. The parameterization |^4| with bmin = &mm=16.97 fm would lead to smaller cross 
sections cxi = 4130 mb and (T2 = 319 mb. 

We found the ratio of (Pm=+i + Pm=-i) / Pm=o = 9.4 for the excitation of the GDR in the 
experiment of Ref. They quote the value 28 in their calculations and fit the gamma-ray 

angular distribution according to this value. We think that this result could somewhat change 



the extracted value of the GDR x GDR cross section which is quoted in Ref. . 

Using the formalism shown of section 3.2.4 we find that the cross sections for the excitation 
of three-phonon states in the experiment of Schmidt et al. |^D[ is equal to 19.2 mb (with 
bmin = 13.3 fm) while it is equal to 117 mb (with bmin = 14.2 fm) for the experiment of Ritman 
et al. |39|. The identification of these resonances is therefore more difficult, but possible with 
the present experimental techniques. Using the same arguments leading to Eq. (|174|) we find 
for the reduced matrix elements, in obvious notations, = 3(a;io/a;32)|£'ioP, which we 



used in our calculation. We assumed that ujiq/uj2,2 — ojiq/uj2i. These enhancement factors 
for the excitation of higher phonon states are very important to explain the magnitude of the 
cross sections. The anharmonic effects, suggested in [|0| to explain the large excitation of 
double GDR, are expected to be small since the mixing of single- and double- phonon states 
is forbidden by the angular momentum and parity. The main anharmonic effect, apart from 
the weak coupling of the double GDR with L = 2 to GQR, is the IBM-like scattering of 
dipole phonons which splits L = and L = 2 states but hardly changes excitation and decay 
properties. 

Another important question is related to the expected width of the multiphonon states. 
Early estimates |]15[ presented in Sect. 3.2.2 indicated that these widths should scale as r„ = 
nTi. The experiments show however that a scaling as r„ = ^/n^l is more appropriate, at least 
for the double GDR. We next address in detail different aspects of physics responsible for the 
width of the double phonon state. 



3.3 General arguments on the width of the double phonon state 

Here we discuss in qualitative terms the problem of the width of a collective state which can be 
thought of as being created by the excitation of two quanta in a complex many-body system. 
We assume that the genuine decay to continuum is of minor importance at the given excitation 
energy. Therefore we focus on the damping width which comes from the fact that the collective 
mode is a specific coherent superposition of simple configurations (for instance, of a particle-hole 
character) rather than a pure stationary state. 

In the actual excitation process the predominant mechanism is that of the sequential one- 
phonon excitation. Under our assumption that the sum rule is saturated by the GR the 
intermediate states contribute to this process as far as they contain a significant collective 
component. Therefore the interference of many incoherent paths can be neglected so that we 
are interested in the shape P{E) of the excitation function at a given energy E — Ei -\- E2 
which can be obtained as a convolution of the single phonon excitation functions, 

P{E) = J dE^dE2Pi{E^)P2{E2)5{E - E^ - E2). (194) 

The same shape should be revealed in the deexcitation process. 

In this formulation the problem is different from what is usually looked at when one is 
interested, for example, in sound attenuation. In such classical problems the conventional 
exponential decrease of the wave intensity does not correspond to the decay of the state with 
a certain initial number of quanta. Contrary to that, here we have to compare the damping 
rates of individual quantum states with the fixed number of quanta, single- and double-phonon 
states in particular. 

We have to mention also that in the nuclear GR case quantum effects are more pronounced 
since the temperature corresponding to the relevant excitation energy is less than huo whereas 
in the measurements of the attenuation of the zero and first sound in the macroscopic Fermi- 
liquid |^8[ the situation is always inverse and the quantum limit is hardly attainable. (In nuclear 
physics the classical case can be studied with low-lying quadrupole vibrations). 

Independently of specific features of nuclear structure (level density, A-dependence, shell 
effects, finiteness of the system leading to the linear momentum nonconservation and, therefore, 
to the estimate of the available phase space which could be different from that for infinite matter, 
and so on) we can try to make several comments of general nature. 



If the anharmonic effects could be considered to be small we can assume that the phonons 
decay independently by what can be described, using the language of stationary quantum 
mechanics, as mixing to complex background states. The decay rate ri(e) of an individual 
quasiparticle (elementary excitation) with energy e depends on the background level density 
and, whence, on the excitation energy. The decay of a state with n quasiparticles occurs 
as far as one of the constituents decays. It implies the simple estimate of the width r„ of 
the n-quantum state, ~ nTi{E/n). For the decay of typical many particle- many hole 
configurations ||59|, ISO, ^ one usually takes the Fermi-liquid estimate Fi(e) oc which leads to 



F„ oc oc E"^/"^ since the average number of quasiparticles in a typical thermal configuration 
at temperature T is n oc T. This estimate agrees with data. In the case of the pure n-phonon 
state E /n = hu which results in the ratio r„ = F„/Fi ~ n. 

Thus, the simplest line of reasoning favors the width of the double GR to be twice as big 
as the width of the single GR. At the first glance, this estimate is especially reasonable for 
the giant dipole since here the anharmonic effects, determining the whole pattern of low-lying 
vibrations, are expected to be very weak. Angular momentum and parity conservation forbids 
cubic anharmonicity which would mix single- and double-quantum states and influence both 
excitation cross sections and spreading widths. The main anharmonic term, apart from the 
mentioned in Sect. 3.2.2 weak mixing of the giant quadrupole to the double dipole state with 
L = 2, probably corresponds to the phonon scattering similar to that in the IBM. It results in 
the shift of the double-phonon state from 2huj and splitting of L = and L = 2 states hardly 
changing the decay properties. Experimentally, the energy shift seems to be rather small. 

There are also other arguments for the width ratio r2 = 2. In our calculation of cross sections 
we assumed the BW shape (|182| ) of strength functions ( p,81|) . If the sequential excitation is 
described by the BW functions Pi{Ei) with the centroid at e and the width F, and P2{E2j 
with corresponding parameters e' and F', the convolution ( p.94| ) restores the BW shape with 
the centroid at e + e' and the total width F + F'. For identical phonons it means that the width 
ratio r2 = 2. 

As we mentioned in Sect. 3.2.2, the BW shape of the strength function is derived analytically 



within the simple model ||5^ of coupling between a phonon and complex background states. 
One diagonalizes first the Hamiltonian in the subspace of those complex states and get their 
energies ev If the underlying dynamics is nearly chaotic, the resulting spectrum will show 
up level repulsion and rigid structure similar to that of the Gaussian Orthogonal Ensemble 
(GOE), with the mean level spacing d. Roughly speaking, one can assume the equidistant 
energy spectrum. The collective phonon |1 > at energy Ei is coupled to those states and 
corresponding matrix elements Viy are assumed to be of the same order of magnitude (much 
larger than the level spacing d) for all states |z/ > in the large energy interval around the 
collective resonance. Then the energies of the stationary states (final states |/ > in the notations 
of previous sections) are the roots E = Ej of the secular equation 

FiE)^E-E^-J2^^ = 0, (195) 
and the distribution of the collective strength, Eq. ( |179| ), 

ICF^P = [dF/dE]-^l^^. = [1 + E jEj^^~' ^^^^^ 



reveals the BW shape (|182|) and the "golden rule" expression ( 



I83|) for the width Fi. 



We can repeat the procedure for the double phonon state. Phonons of different kind would 
couple to different background states with different level spacing and coupling matrix elements. 
It corresponds to independent decay leading as we discussed above to r2 = F + T'. For the 
identical phonons, we should take into account that the double phonon state |2 > is coupled to 
the states "single phonon + background" and the background states here are the same as those 
determining the width of the single phonon state |1 >. This picture is in accordance with the 
famous Axel-Brink hypotheses. Therefore the expression for the width, Eq. (|183|) , contains the 
same level density whereas all coupling matrix elements for the transition to a complex state 
\u > (plus a remaining phonon) have to be multiplied by the Bose factor, = v^Vi^,. Thus, 
we come again to r2 = 2. 

The approach of the proceeding paragraph can be slightly modified by introducing explicitly 
coupling via a doorway state ||6^ or GOE internal dynamics [Q. In both cases the Bose factor 
^/2 leads to the same result r2 = 2. 

In addition, the collective resonance might be further broadened by the coupling to low- 
lying collective vibrational or rotational modes. For example, in the simplest model where the 
dipole phonon radiates and absorbs low energy scalar quanta, it is easy to show that, in the 
stationary cloud of scalar quanta, their average number, which determines the fragmentation 
region of the dipole mode, is proportional to the squared number of dipole phonons. Hence it 
gives a large width ratio r2 = 4. For the nuclei where actual data exist, this is not important 
since they are rather rigid spherical nuclei with no adiabatic collective modes. 

On the other hand, one can present some arguments in favor of the width ratio r2 = -\/2 
which apparently is preferred by the existing data. 

First of all, this value follows from the convolution ( 194 ) of Gaussian distribution functions 
(instead of BW ones). Of course, this is the inconsistent approach since the experimentalists use 
BW or Lorentzian fit. But one can easily understand that the result r2 = is not restricted to 
Gaussian fit. For an arbitrary sequence of two excitation processes we have < E >=< E1+E2 > 
and < E^ > = < {Ei + E2Y >] for uncorrelated steps it results in the addition of fluctuations 
in quadrature, {AEY = {AEiY + {AE2Y. Identifying these fluctuations with the widths up 
to a common factor, we get for the identical phonons F^ = 2F^, or r2 = a/2. 

The same conclusion will be valid for any distribution function which, as the Gaussian one, 
has a finite second moment, contrary to the BW or Lorentzian ones with the second moment 
diverging. In some sense we may conclude that, in physical terms, the difference between r2 = 2 
and r2 = is due to the different treatment of the wings of the distribution functions which 
reflect small admixtures of far remote states. 

In the standard model of the strength function all remote states are coupled to the 
collective mode equally strong. This is obviously an unrealistic assumption. The shell model 
(more generally, mean field) basis is the "natural" one [^] for estimating a degree of complexity 
of various states in a Fermi system at not very high excitation energy. In this representation 
matrix elements of residual interaction couple the collective state (coherent superposition of 
particle-hole excitations found for example in the framework of the RPA) only to the states of 
the next level of complexity (exciton class). Those states, in turn, become mixed with more 
complex configurations. This process proliferates and each simple state acquires its spreading, 
or fragmentation, width 2a = Nd where stands for a typical number of stationary states 
carrying the noticeable weight of the ancestor state and the level spacing d is basically the same 
as in the mean field approximation. Inversely, A^ can be viewed as the localization length of a 
stationary complex state in the mean field basis. 



In the stochastic hmit the local background dynamical properties can be modeled by those 
of the GOE with the semicircle radius a. This intrinsic spreading width a, which is expected 
to be of the order of magnitude of typical matrix elements of the original residual interaction 
between simple configurations, is the dynamical scale missed in the standard model which 
corresponds to the limit a —>■ oo. The existence of this intrinsic scale can be associated with 
the saturation of the width of a single GR at high temperature. 

The standard model supposedly is valid for the spreading width F small in comparison with 
a. Because of the relatively weak interaction leading to the isospin impurity, this is the case for 



the isobaric analog states (IAS) |^ ^ where typical spreading widths are less than 100 keV. 
This approach allows one to explain, at least qualitatively, small variations of the spreading 
widths of the IAS. The tunneling mixing of superdeformed states with the normal deformed 
background presents an extreme example of the small spreading width. However, in the case 
of GR the situation might be different. 

To illustrate the new behavior in the opposite case of F > a, we can imagine the limit of the 
almost degenerate intrinsic states with very strong coupling to a collective mode. (The actual 
situation presumably is intermediate). Assuming that the unperturbed phonon state has an 
energy in the same region, one can easily see from Eqs. ( |195| ) and (|196|) that the coupling results 
in the appearance of the two collective states sharing evenly the collective strength and shifted 



symmetrically from the unperturbed region by AE = ±yJ2u K?- The physical reason is evident: 
the interaction of the background states through the collective mode creates a specific coherent 
superposition which is hybridized with and repelled from the original collective state. The 
similar effect was discussed in different context in |Q and observed in numerical simulations 



69 . The well known doubling of the resonance peak at the passage of a laser beam through a 



cavity containing a two-level atom is the simplest prototype of such a phenomenon. 



In this limit one gets the effective width of collective response 2AE = 2^ N <V^ > = 

2^aVs/T: where F^ is the standard spreading width ( p.83| ). This effective width is linearly 
proportional to the average coupling matrix element. Therefore it should increase by factor 
\/n when applied to a ra— phonon collective state. Thus, we anticipate in this limit r2 = \/2. 
One may say that the phonons do not decay independently being correlated via common decay 
channels. In the literature the similar result, due to apparently the same physical reasons, was 
mentioned in ||70| referring to the unpublished calculations in the framework of the second RPA. 



3.4 Coupled-channels calculations with inclusion of the GR widths 

We have seen that the excitation probabilities of excitation of single and double giant resonances 
are quite large. It is worthwhile to study the excitation process with a coupled channels 
calculation and compare to the other approximations. We will now study this effect by using the 
Coupled-Channels Born approximation. This approximation was used in Ref. ||71| to describe 
the excitation of the double giant resonance in relativistic heavy ion collisions. It is based on 
the idea that in such cases only the coupling between the ground state and the dominant giant 
dipole state has to be treated exactly. The reason is that the transitions to giant quadrupole and 
to the double-phonon states have low probability amplitudes, even for small impact parameters. 
However, an exact treatment of the back-and-forth transitions between the ground state and 
the giant dipole state is necessary. This leads to modifications of the transitions amplitudes 
to the remaining resonances, which are populated by the ground state and the GDR. In Ref. 



71 the application of the method was limited to the use of an schematic interaction, and the 



magnetic substates were neglected. These deficiencies are corrected here. The method allows 
the inclusion of the width of the giant resonances in a very simple and straightforward way. It 
will be useful for us to compare with the Coupled-Channels calculations with isolated states, as 
we described in the previous sections. Fig. ^ represents the procedure. The GDR is coupled 
to the ground state while the remaining resonances are fed by these two states according to 
first order perturbation theory. The coupling matrix elements involves the ground state and 
a set of doorway states {V^xil where n specifies the kind of resonance and A/i are angular 
momentum quantum numbers. The amplitudes of these resonances in real continuum states 
are 



(197) 



where 0(e) denotes the wavefunction of one of the numerous states which are responsible for 
the broad structure of the resonance. In this equation e = — En, where E^ is the excitation 
energy and En is the centroid of the resonance considered. 

As we have stated above, in this approach we use the Coupled-Channels equations for the 
coupling between the ground state and the GDR. This results in the following Coupled- Channels 
equations: 

th do(t) = j de< 0(e)|P2 > < I?i;]|Vi,i,^(t)|0 > exp { - '-{E^ + e)t] a%{t) 

= E Jde a^'\e) V^^'^t) exp { - ^{E, + e)t} a^^lit) , (198) 



and 



aS.W = [("^'He) ^rn^)]* exp{z(i?i + e)tA} ao(t) . 



(199) 



Above, (n = 1) stands for the GDR, denotes the occupation amplitude of the ground state 
and a[^i^ the occupation amplitude of a state located at an energy e away from the GDR 
centroid, and with magnetic quantum number /i (/i = —1,0,1). We used the short hand 
notation V^^^\t) =< V^^hVEi,MO >■ 



Integrating Eq. ( |199|) and inserting the result in Eq. (|198|) , we get the integro-differential 
equation for the ground state occupation amplitude 



do(t) 



He)\' 



X / dt' [v^(°'Ht')]* exp { - t{Ei + e){t - i!)iri] ao{t') 



(200) 



where we used that a^^l^{t = — oo) = 0. To carry out the integration over e, we should use 
an appropriate parametrization for the doorway amplitude a^^\e). A convenient choice is the 
Breit-Wigner (BW) form 



|a«(e)P 



1 

2^ 



+ r?/4 



(201) 



where Fi is chosen to fit the experimental width. In this case, this integral will be the simple 
exponential 



/ de |a(^)(e)|' exp { 



— I 



{El + e)t 



} = exp { 



.{Ei-iTi/2)t 



}■ 



(202) 



A better agreement with the experimental hne shapes of the giant resonances is obtained 
by using a Lorentzian (L) parametrization for |a*^^''(e)p, i.e., 



2 

TT 



{El - Elf + TlEl 



(203) 



where Ex = Ei + e. The energy integral can still be performed exactly but now it leads to 
the more complicated result 



j de \a^^\t)\^ exp{- 



.{Ei + e)t_ 

h 



2Ei' 



.{E^-iTi/2)t 

h 



} 



+ AC{t) 



(204) 



where AC (t) is a non-exponential correction to the decay. For the energies and widths involved 
in the excitation of giant resonances, this correction can be shown numerically to be negligible. 
It will therefore be ignored in our subsequent calculations. After integration over e, Eq. ( pOOj ) 
reduces to 



ao(t) 



J—oo '- ■' ^ 



.{Er-iT,/2){t-t') 



} ao(t') 



(205) 



where the factor Si is 5i = 1 for BW-shape and Si = 1 — iTi/2Ei for L-shape. 

We can take advantage of the exponential time-dependence in the integral of the above 
equation, to reduce it to a set of second order differential equations. Introducing the auxiliary 
amplitudes Af^{t), given by the relation 



ao{t) = 1 + Y.Mt) 



(206) 



with initial conditions A^{t = — oo) = 0, and taking the derivative of Eq. ( P05| ), we get 



Ti 

El — I — 
2 



h 



Mi) + Si 



. 



(207) 



Solving the above equation, we get aoit). Using this amplitude and integrating Eq. ( [I99| ), 
one can evaluate ae^i^(t). The probability density for the population of a GDR continuum 
state with energy E^ in a collision with impact parameter b, Pi{b,Ex), is obtained trough the 
summation over the asymptotic {t oo) contribution from each magnetic substate. We get 



dt' exp{tEj} [Vi''\t')] ao(t') 



(208) 



where \a^^\Ex — Ei)\'^ is given by Eq. ( p01| ) or by Eq. ( |203| ), depending on the choice of the 
resonance shape. 

To first order, DGDR continuum states can be populated through i^^2-coupling from the 
ground state or through £'l-coupling from GDR states. The probability density arising from 
the former is given by Eq. ( p08| ), with the replacement of the line shape ja'-^-'p by its DGDR 
counterpart |a*^^^P (defined in terms of parameters E2 and r2) and the use of the appropriate 



coupling- matrix elements V^^'^^t) with the E2 time dependence given by (pQ]). On the other 
hand, the contribution from the latter process is 



P2{h,E,) 




We should point out that Eq. ( ^09| ) is not equivalent to second-order perturbation theory. This 
would be true only in the limit ao(t) — > 1. In this approach, a^it) ^ 1, since it is modified 
by the time-dependent coupling to the GDR state. This coupling is treated exactly by means 
of the Coupled-Channels equations. We consider that this is the main effect on the calculation 
of the DGDR excitation probability. This approach is justified due to the small excitation 
amplitude for the transition 1 — > 2, since ai(t) <^ ao(^)- 

Equations similar to ( |208| ) can also be used to calculate the GQRj^ and GQRj^ excitation 
probabilities, with the proper choice of energies, widths, and transition potentials (e.g., VB2(t), 
or VAr2(t), or both). 

In the next section we will apply the results of this section to analyze the effect of the widths 
of the GR's in a Coupled-Channels approach to relativistic Coulomb excitation. 



3.4.1 Zero-width calculations 

We consider the excitation of giant resonances in ^osp]^ projectiles, incident on ^osp]^ targets 
at 640A MeV, which has been studied at the GSI/SIS, Darmstadt 0, ^. For this system 
the excitation probabilities of the isovector giant dipole (GDRj^,) at 13.5 MeV are large and, 
consequently, high order effects of channel coupling should be relevant. To assess the importance 
of these effects, we assume that the GDR state depletes 100% of the energy-weighted sum-rule 
and neglect the resonance width. 

As a first step, we study the time evolution of the excitation process, solving the Coupled- 
Channels equations for a reduced set of states. We consider only the ground state (g.s.) and 
the GDR. The excitation probability is then compared with that obtained with first order 



perturbation theory. This is done in Fig. where we plot the occupation probabilities of the 
g.s., |ao(t)P, and of the GDR, |ai(t)p, as functions of time, for a collision with impact parameter 
6 = 15 fm. As discussed earlier, the Coulomb interaction is strongly peaked around t = 0, with 
a width of the order At ~ 6/7^. Accordingly, the amplitudes are rapidly varying in this 
time range. A comparison between the CC-calculation (solid line) and first order perturbation 
theory (dashed line) shows that the the high order processes contained in the former lead 
to an appreciable reduction of the GDR excitation probability. From this figure we can also 
conclude that our numerical calculations can be restricted to the interval —10 < r < 10, where 
r = {'yv/b) t is the time variable measured in natural units. Outside this range, the amplitudes 
reach asymptotic values. 

It is worthwhile to compare the predictions of first order perturbation theory with those of 
the harmonic oscillator model and the CC calculations. In addition to the GDR, we include 
the following multiphonon states: a double giant dipole state ([GDR^]^) at 27 MeV, a triple 
giant dipole state ([GDRj^,]'^) at 40.5 MeV, and a quadruple giant dipole state ([GDRji,]"^) at 54 
MeV. The coupling between the multiphonon states are determined by boson factors, i.e., for 
0-^1 and n — 1 — > n 



<n-l\\VE/NA\n>Y = n\<Q\\VE/NA\\^>\ ■ (210) 



Table 6: Transition probabilities at 6 = 14.3 fm, for the reaction ^^^Pb + ^°^Pb at 640A MeV. 
A comparison with first order perturbation theory is made. 



Trans. 




1st pert. th. 


c.c. 


g.s. — > 


g.s. 




0.515 


g.s. — > 




0.506 


0.279 


g.s. — > 


GQRjs 


0.080 


0.064 


g.s. — > 




0.064 


0.049 


g.s. — > 


[GDR,J2 


0.128 


0.092 



Direct excitations of the multiphonon states from the g.s. are not considered. The angular 
momentum addition rules for bosons yields the following angular momentum states: L = 
and 2, for the [GDR]^ state; L = 1, 2, and 3, for the [GDR]^ state; and L = 0, 1, 2, 3, 
and 4, for the [GDR]^ state. We assume that states with the same number of phonons are 
degenerate. In Table ^ we show the resulting cross sections. The excitation probabilities and 
the cross section were calculated with the formalism of section 3.4. The integration over impact 
parameter was carried out in the interval bmin < b < oo. As we discuss below, the low-6 cut-off 



value [0 bmin = 14.3 fm mocks up absorption effects. We have checked that the CC results are 
not significantly affected by the unknown phases of the transition matrix elements. Since the 
multiphonon spectrum is equally spaced, and the coupling matrix-elements are related through 
boson factors (as in Eq. ( |210| )), the harmonic oscillator and the CC cross sections should be 
equal. In fact the numerical results of these calculations given in the table are very close. We 
also see that the excitation cross sections of triple- and quart ic-phonon states are much smaller 
than that for the [GDR]^. Therefore, we shall concentrate our studies on the [GDR]^, neglecting 
other multiphonon states. 

Next, we include the remaining important giant resonances in ^°^Pb. Namely, the isoscalar 
giant quadrupole (GQRjs) at 10.9 MeV and the isovector giant quadrupole (GDRj„) at 22 MeV. 
Also in this case, we use 100 % of the energy-weighted sum rules to deduce the strength 
matrix elements. In Table |^, we show the excitation probabilities in a grazing collision, with 
b = 14.3 fm. We see that first order perturbation theory yields a very large excitation probability 
for the GDRj^ state. This is strongly reduced in a c.c. calculation, as we have already discussed 
in connection with Fig. The excitations of the remaining states are also influenced. They 
are reduced due to the lowering of the occupation probabilities of the g.s. and of the GDRj^ 
state in the c.c. calculation. As expected, perturbation theory and c.c. calculations agree 
at large impact parameters, when the transition probabilities are small. For the excitation of 
the [GDRj^,]^ state we used second-order perturbation theory to obtain the value in the second 
column. The presence of the GQRj^ and the GQRjt, influence the c.c. probabilities for the 
excitation of the GDR and the [GDRj^]^, respectively. 



3.4.2 Effect of resonance widths 

We now turn to the influence of the giant resonance widths on the excitation dynamics. We 
had considered this in Sect. 3.2. But, now we show that the coupled-channels effects lead to 
important quantitative modifications of the results. We use the CCBA formalism developed in 
Sect. 3.4. Schematically, the CC problem is that represented Fig. E3[ As we have seen above, 



Table 7: Centroid energies and widths of the main Giant Resonances in ^°^Pb. 





GDR 


DGDR 


GQRjs 


GQRto 


E [MeV] 


13.5 


27.0 


10.9 


20.2 


r [MeV] 


4.0 


5.7 


4.8 


5.5 



the strongest couphng occurs between the g.s. and the GDR. 

In Fig. we show the excitation energy spectrum for the GDR, the DGDR (a notation for 
the [GDRj^]^), GQRjs and GQRj^. The centroid energies and the widths of these resonances 
are hsted in Table |^. The figure show excitation spectra obtained with both Breit-Wigner 
(BW) and Lorentzian (L) line shapes. One observes that the BW and L spectra have similar 
strengths at the resonance maxima. However, the low energy parts (one or two widths below 
the centroid) of the spectra are more than one order of magnitude higher in the BW calculation. 
The reason for this behavior is that Coulomb excitation favors low energy transitions and the 
BW has a larger low energy tail as compared with the Lorentzian line shape. The contribution 
from the DGDR leads to a pronounced bump in the total energy spectrum. This bump depends 
on the relative strength of the DGDR with respect to the GDR. In Fig. |2^, we show the ratio 
o"dgdr/o"gdr as a function of the bombarding energy. We observe that this ratio is roughly 
constant in the energy range Eiab/A = 200 — 1000 MeV and it falls beyond these limits. This 
range corresponds to the SIS-energies at the GSI-Darmstadt facility. 

We now study the influence of the resonance widths and shapes on the GDR and DGDR 
cross sections. This study is similar to that presented in Ref. |]7T[, except that we now have 
a realistic three dimensional treatment of the states and consider different line shapes. In the 
upper part of Fig. denoted by (a), we show o"gdr as a function of Fgdr, treated as a free 
parameter. We note that the BW and L parameterizations lead to different trends. In the BW 
case the cross section grows with Fqdr while in the L case it decreases. The growing trend 
is also found in Ref. which uses the BW line shape. The reason for this trend in the 
BW case is that an increase in the GDR width enhances the low energy tail of the line shape, 
picking up more contributions from the low energy transitions, favored in Coulomb excitation. 
On the other hand, an increase of the GDR width enhances the doorway amplitude to higher 
energies where Coulomb excitation is weaker. In Fig. ^ (b) and (c), we study the dependence 
of o"gdr on Fgdr- In (b), the DGDR width is kept fixed at the value 5.7 MeV while in (c) it is 
kept proportional to ctgdr, fixing the ratio FnGDR/rGDR = V^- The first point to be noticed is 
that the BW results are systematically higher than the L ones. This is a consequence of the 
different low energy tails of these functions, as discussed above. One notices also that Cdgdr 
decreases with Fgdr both in the BW and L cases. This trend can be understood in terms of the 
uncertainty principle. If the GDR width is increased, its life-time is reduced. Since the DGDR 
is dominantly populated from the GDR, its short life-time leads to decay before the transition 
to the DGDR. 

To assess the sensitivity of the DGDR cross section on the strength of the matrix elements 
and on the energy position of the resonance, we present in Table |^ the cross sections for the 
excitation of the GDR, DGDR, GQRj^ and GQRj„, obtained with the CCBA approximation 
and 100% of the sum-rules for the respective modes. In this calculation we have included the 
strong absorption, as explained in Sect. 2.8. For comparison, the values inside parenthesis 
(and brackets) of the DGDR excitation cross section include a direct excitation of the L=2 



Table 8: Cross sections (in mb) for the excitation of giant resonances in lead, for the reaction 
208pb + 208pb at 640A MeV. See text for details. 



GDR 


DGDR 






2704 


184 (199) [198] 


347 


186 



DGDR state. We assumed that 20% of the E2 sum rule could be allocated for this excitation 
mode of the DGDR. The cross sections increase by less than 10% in this case. The value inside 
parenthesis (brackets) assume a positive (negative) sign of the matrix element for the direct 
excitation. 

Since the excitation of the DGDR is weak, it is very well described by Eq. ( P09| ) and the 
DGDR population is approximately proportional to the squared strength of V^^'^\ Therefore, 
to increase the DGDR cross section by a factor of 2, it is necessary to violate the relation 
EuGDK = '^Eqy^t^ by the same factor. This would require a strongly anharmonic Hamiltonian for 
the nuclear collective modes, which would not be supported by traditional nuclear models [|^. 
Arguments supporting such anharmonicities have recently been presented in Ref. |j73|, 0, ^ . 
Another effect arising from anharmonicity would be the spin or isospin splitting of the DGDR. 
Since the Coulomb interaction favors lower energy excitations, it is clear that a decrease of 
the DGDR centroid would increase its cross section. A similar effect would occur if a strongly 
populated substate is splitted to lower energies. To study this point, we have varied the energy 
of the DGDR centroid in the range 20 MeV < E^gdr < 27 MeV. The obtained DGDR cross 
sections (including direct excitations) are equal to 620 mb, 299 mb and 199 mb, for the centroid 
energies of 20 MeV, 24 MeV, and 27 MeV, respectively. Although the experimental data on 
the DGDR excitation [0, |, Q seem to indicate that i?DGDH ~ 2_E'gdr, a small deviation (in the 
range of 10%-15%) of the centroid energy from this value might be possible. However, the data 
are not conclusive, and more experiments are clearly necessary. We conclude, that from the 
arguments analyzed here, the magnitude of the DGDR cross section is more sensitive to the 
energy position of this state. The magnitude of the DGDR cross section would increase by a 
factor 2 if the energy position of the DGDR decreases by 20%, as found in Ref. |]73|, |7^, 
due to anharmonic effects. In Ref. [|T^ one obtained (Tdgdr = 620 mb, 299 mb, and 199 mb 
for the centroid energies of -Edgdr = 20 MeV, 24 MeV and 27 MeV, respectively. This shows 
that anharmonic effects can play a big role in the Coulomb excitation cross sections of the 



DGDR, depending on the size of the shift of Eugou- However, in Ref. the source for 
anharmonic effects were discussed and it was suggested that it should be very small, i.e., A^'^^E 
— it; —OK ~ n 

The anharmonic behavior of the giant resonances as a possibility to explain the increase 
of the Coulomb excitation cross sections has been studied by several authors (see also Ref. 
76 1 , and references therein). It was found that the effect is indeed negligible and it could be 
estimated [0 as A^^)^ < Egdr/(50.A) ~ A-^l^ MeV. 

One attempt to explain the larger experimental cross section is to include contributions 
from the excitation of a single coherent phonon on "hot" fine structure states (Brink- Axel 
mechanism). Recently [^, |7^ this has been done through two different approaches. In the 
first one |7^, the nucleus is described as a collective harmonic oscillator interacting with a 



set of oscillators representing statistical degrees of freedom. In the second |]78[, a statistical 



approach along the lines proposed by Ko ||79| (see also [ITH]) is used. These works indicate that 



the Brink- Axel mechanism should play an important role, being able to explain, in part, some 
discrepancies between theory and experimental cross sections. Further work along these lines 



were published by, Hussein and collaborators [BC, 31, B2|. In a recent publication [R3i, the 



influence of the isospin structure of the double giant resonance was studied in details. It was 
shown that this structure also leads to an enhancement of the calculated cross sections. 

The calculations discussed so far are based on macroscopic properties of the nuclei, sum 
rules, etc. Now we show that, in order to obtain a better quantitative description of double giant 
resonances it is necessary to include the internal degrees of freedom of the nuclei appropriately. 
We will discuss this next. But, we flrst describe the formalism that we will use for this purpose. 



4 Description of one- and multi-phonon excited states 
within the Quasiparticle-Phonon Model 

4.1 The model Hamiltonian and phonons 

The Hamiltonian, H, of the Quasiparticle-Phonon Model (QPM) (see Refs. p^ , p5| , p6| for 
more details) is introduced on basis of physical ideas of nucleons moving in an average fleld 
and interacted among each other by means of a residual interaction. Schematically it can be 
written in the form: 

H = Hg.p. + Hpair + Hr.i. ■ (211) 

We limit ourselves here only by the formalism for even-even spherical nuclei. The flrst term 
of Eq. ( 211 ), Hs.p.i corresponds to the average fleld for neutrons (n) and protons (p). In the 



second-quantized representation it can be written in terms of creation (annihilation) {o-jm) 
operators of particles on the level of the average fleld with quantum numbers j = [n,l,j] and 
m as following 

n,p 

Hs.p.=T.T.E,al^a^m, (212) 

T j,m 

where Ej is the energy of the single particle level degenerated in spherical nuclei by magnetic 
quantum number m. The second term of Eq. (|211|) , Hpair, corresponds to residual interac- 



tion responsible for pairing in non-magic nuclei. In the QPM this interaction is described by 
monopole pairing with a constant matrix element G^^^ 



n,p 



Hparr = E ^^r^ E \/(2j + 1) " (2j' + 1) J 00 K'-™'a,'^']oO , (213) 



[a'j'ay]xp — E ^jmj'm'^tm'^t'm' (214) 
m,m' 

where Cjl^j,^, is the Clebsch-Gordan coefficient. Since the QPM is usually applied for a descrip- 
tion of properties of medium and heavy nuclei with a filling of different subshells for neutrons 
and protons, the neutron-proton monopole pairing is neglected. The residual interaction, Hj.,i., 
is taken in the QPM in a separable form as a multipole decomposition. Its part in the particle- 
hole channel can be written as 

^^ir'^ = E E(4'^ + p4"^)MUr)MxM) , (215) 



where 

^^0(1) ^-^^ 

the model parameters which determine the strength of isoscalar (isovector) 
residual interaction. The multipole operator M;^(r) has the form 

MU^)= E OmhV;(r)n^(fi)|jW)a+„a,w (216) 

j,m,j',m' 

for the natural parity states and the form 

M^r) = ^ (jm I t^flir)[a ■ Y^rr^mx, I f^'Hm ^fm' (217) 

j,m,j' ,m' ,lmi 

for the unnatural parity states. The function fl{r) is a radial formfactor which in actual 
calculations is taken either as or as a derivative of the central part of the average field: 
/^(r) = dU'^{r)ldr. The value r = — 1(+1) corresponds to neutrons (protons). We will 
not consider here the residual interaction in the particle-particle channel which is the most 
important for the description of two-nucleon transfer reactions. 

The basic QPM equations are obtained by means of step-by-step diagonalization of the 
model Hamiltonian ( 211 ). On the first step its first two terms ( pi2| ) and ( |213| ) are diagonal- 
ized. For that the Bogoliubov's canonical transformation from particle creation (annihilation) 
operators to quasiparticle creation (annihilation) operators (a^m) is applied: 

4m = «i<m + ■ (218) 

The ground state of even-even nucleus, |)q, is assumed as a quasiparticle vacuum: ajm\)q = 0. 
Then the energy of the ground state is minimized: 

6{{\ Hs.p. + Hpair I), + E/^i(«? + ^1 - 1)} = , (219) 

j 

where Hj are Lagrange coefficients. The result of this minimization are the well-known BCS 
equations solving which one obtains correlation functions Cj- = G^^ UjVj and chemical po- 
tentials Xr for neutron and proton systems. The coefficients of the Bogoliubov transformation 
Uj and Vj can be calculated from these values as following: 



where ej is the quasiparticle energy: 



= + [E, - A.]2 . (221) 

In magic nuclei the BCS equations yield a zero value for the correlation function and the 
position of the chemical potential in the gap between particles and hole is uncertain. This 
results in vanishing of monopole pairing correlations and the Bogoliubov's coefficients Uj{vj) 
equal to 0(1) for holes and to 1(0) for particles, respectively. 

After diagonalization of the first two terms of the model Hamiltonian ( |211| ) they can be 
written as following: 

n,p 

Hg.p. ~\~ Hpaif ^ ^ ^ ^ ^j^jm^jiTT- (222) 
T j,m 

and the multiple operator (|216|) in terms of quasiparticle operators has the form 



fW { A+) 

MUr) = E -f- {[c^M'W + (-1)^-'^K'«,]a-,) - v^)BAjf- A/i)} ; (223) 



.;v2A 



where f^^) =< j\\i^fx{r)Yx{^)\\j' > is the reduced matrix element of the multipole operator. 



We also introduced the following combinations of the Bogoliubov's coefficients: u^^) = UjVji ± 
UjiVj and f j^'' = UjUji =F VjVf to be used below. 

We have determined the ground state of even-even nuclei as the quasiparticle vacuum. In 
this case, the simplest excited states of nucleus are two-quasiparticle states, ttjm'^/m'l)'?; which 
correspond to particle-hole transitions if monopole pairing vanishes. Two fermion quasiparticle 
operators couple to the total integer angular momentum corresponding to the Bose statistics. 
Thus, it is convenient to project the bi-fermion terms [oi'ja^i]\^ and [«j'aj]A-/i in Eq- (^2^) 
into the space of quasi-boson operators. Following this boson mapping procedure we introduce 
the phonon operators of the multipolarity A and projection yU as following: 

= ^ E E - (-i)'^'^v'';'K-'«.]a-4 • (225) 

^ ^ n' 

The total number of different phonons for the given multipolarity A should be equal to the sum 
of neutron and proton two-quasiparticle states coupled to the same angular momentum. The 
index i is used to number these different phonons. 

The coefficients ip'^^t and yjj'j, of the linear transformation ( p25| ) one obtains by diagonaliza- 
tion of the model Hamiltonian in the space of one-phonon states, Q~x^i\)ph- This can be done 
for example by applying again the variation procedure: 

'^{(1 Qx,^HQl^, \U - ^E{(^^;')' - is^'f] - 2]} = . (226) 



2 



3f 



It yields the well-known equations of the random-phase approximation (RPA) which for the 
case of the separable form of the residual interaction in p/i-channel may be written as 



(227) 



where the following notation have been used 



2A + _ 



The determinant equation (|227|) is a function of the nucleus excitation energy uj. Solving this 
equation for each multipolarity A'^ one obtains the spectrum of nuclei one-phonon excitation lOxi. 
The index i in the definition of the phonon operator ( P25| ) gets the meaning of the order number 



of the solution of Eq. ( |227| ). The fermion structure of phonon excitation, i.e. the amplitudes 
tl) and corresponding to the contribution of different two-quasiparticle components to the 
phonon operator, are obtained from the following equation 

^'^^^^^^^ (229) 




where the value 3^^' is determined from normalization condition for phonon operators: 



n,p 



{\Q^,^QU)pk = EE {(^i]')' - (V''}) } = 2 (230) 



and one obtains 



2A + 1 ^ [fe + ^.'P-^^ 



|2 



Equations (|227|) , ( P29|) and (|231|) correspond to natural parity phonons. Similar equations 



are valid for unnatural parity phonons by substituting the reduced spin-multipole matrix ele- 
ment /jy'^ and combination of coefficients of Bogoliubov transformation u^j^) for /^y and u^jjK 
respectively. Also, amplitude ip^j' changes the sign in Eq. (E^UF) for unnatural parity phonons. 



The RPA equations have been obtained under the assumption that the nucleus ground state 
is the phonon vacuum, Qx^i\)ph = 0. This means that the ground state correlations due to the 
last term of the model Hamiltonian, Hr.i, are taken into account. If they are not accounted for 
and the ground state is still considered as a quasiparticle vacuum \)q, one obtains the so-called 
Tamm-Dankov approximation (TDA). The TDA equations can be easily obtained from the 
RPA ones by neglecting backward going amplitudes in the definition of the phonon operator 
(iH), i.e. applying cp^', = 0. 



The relation between the wave functions of the phonon and quasiparticle vacuums is the 
following [Bl]: 



1).^ = l^nexp -^E E(^j:.J-VJ;,.(-l)^-^[«>,t]A.[a;3a+]A-. \ I). (232) 

A in 3132 

\ 3334 

where A/" is a normalization factor. 

For actual numerical calculations one needs to determine the model parameters. The aver- 
age field for neutrons and protons is described in the QPM by phenomenologic Woods-Saxon 
potential: 

The parameters of this potential for different A-mass regions are listed in Table |^ (see, also 
Ref. [^). We usually use RJ^ = R^, aj^ = aj, and Rc = Ro- All single particle levels from the 



bottom are included in calculation. The single particle continuum is approximated by narrow 
quasibound states. This approximation gives a good description of the exhaust of the energy 
weighted sum rules (EWSR) for low values of A in medium and heavy nuclei. For the lead 
region we use the single particle spectrum near the Fermi surface from Ref. which was 



adjusted to achieve a correct description of low-lying states in neighboring odd nuclei. The 
parameters of the monopole G^^ have been fitted to reproduce the pairing energies. 

The parameters of the residual interaction are obtained the following way. The strength of 
the residual interaction for A'^ = 2+ and 3~ is adjusted to reproduce the properties (excitation 
energy and B{EX) value, known from experiment) of the 2f and 3j~ states. Usually it is not 
possible within one-phonon approximation, discussed in this subsection, if sufficiently large 
single particle spectrum is used. When the energy of the lowest excitation is adjusted to the 
experimental value, the RPA equation yield an overestimated collectivity, B{EX) value, for 
this state. And vice versa, if the collectivity of this state is reproduced, the excitation energy 



Table 9: Parameters of Woods-Saxon potential Eq. (|233|) for different A-mass regions. 



A neutrons protons 











^Is 






p 


^Is 




[MeV] 


[fml 

L J 


[fml 

L J 


[MeV] 


[MeV] 


[fm] 

L J 


[fm] 

L J 


[MeV] 


49 


-41.35 


4.852 


0.6200 


-9.655 


-58.68 


4.538 


0.6301 


-9.506 


59 


-46.20 


5.100 


0.6200 


-9.540 


-53.70 


4.827 


0.6301 


-8.270 


91 


-44.70 


5.802 


0.6200 


-9.231 


-56.86 


5.577 


0.6301 


-9.609 


121 


-43.20 


6.331 


0.6200 


-8.921 


-59.90 


6.133 


0.6301 


-10.363 


141 


-45.95 


6.610 


0.6200 


-9.489 


-57.70 


6.454 


0.6301 


-10.069 


209 


-44.83 


7.477 


0.6301 


-8.428 


-60.30 


7.359 


0.6301 


-11.186 



is too high as compared to the experimental value. The situation sufficiently improves when 
the coupling of one-phonon states to more complex configurations is taken into account as 
will be discussed in the next subsection. For the lowest excited state the coupling to complex 
configurations results in the energy shift downwards. Thus, for nuclei not very far from a closed 
shell it becomes possible to achieve a good description of both, the excitation energy and the 
B{EX) value. The ratio between isoscalar and isovector strength of the residual interaction 
is usually fixed as aei^V^o^^ — —1-2 in calculation with the radial formfactor of the multipole 
operator as a derivative of the average field. With this ratio the best description of isovector 
multipole resonances with A > 1 is achieved although the experimental information on these 
resonances is still sparse. For the dipole-dipole residual interaction the strength parameter are 
adjusted to exclude the spurious center of mass motion and to obtain a correct position of the 
GDR centroid. For the phonons with the multipolarity A > 4 the same procedure of adjusting 
the strength parameters as for A'" = 2+ and 3~ cannot be applied. First, it is because the lowest 
states of high multipolarity are much less collective and their properties are more sensitive to 
description of single particle levels near the Fermi surface than to the strength of the residual 
interaction. Second, in many cases the lowest states with A > 4 are either two-phonon states 
or the states with a large admixture of two-phonon configurations, thus, their properties are 
determined by phonons of another multipolarity. For these reasons we use ssq^I = "* for 
even parity phonons and aeo^i = aeg^i for odd parity in calculation with fl{r) = dW^ (r) / dr . 
In fact, the difference between aeo^^ and aeQ^^ does not exceed a few percent with this radial 
formfactor of residual force. 



4.2 Mixing between simple and complex configurations in wave fun- 
ctions of excited states 

Diagonalization of the model Hamiltonian in the space of one-phonon states allows us to write 
it in the form 

H = T. ^xiQUQx^^^ + Hint. , (234) 

Xfii 

H^nt. = E ( [i-^)'-''QU + Qx-,^] E ^^BAjj'^ a - /i) + h.c] , (235) 

^ A^i [ Ji'r V23^r J 



where the origin of the second term in Eq. (|234| ) can be traced back to the last term of multipole 



operator ( P23|) which cannot be projected onto the space of the phonon operators. On the other 
hand, applying Marumori expansion technique , one may expand the operator A — 



/i) ~ a^a in an infinite sum of even-number phonon operators. Keeping only the first term of 
this expansion, the non-diagonal term of the model Hamiltonian, -ffi„i,, in the space of phonon 
operators may be re-written as 

iirni= E f^A't(^OQt^[QA.;..nQA.;...jA;. + /^.C. (236) 

Ai/iiJi 

A2M2*2 

where the matrix element of interaction between one- and two-phonon configurations, U^^l^i^''')^ 
can be calculated by making use of the internal fermion structure of phonons, i.e. ip and (p 
coefficients, and reduced matrix elements of the separable force formfactor, f^-^^j^- It has the 
form 



'(2Ai + l)(2A2 + l) 



n,p 



X 



E E 





\ Ai 


A2 


A 






Ji 


h 




\ Ai 


A2 


A 






32 


Ji 


3 3132^3132 ^ 


\ Ai 


A2 


A 






J3 


J2 
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The upper (lower) sign in each of three terms in Eq. (|237|) correspond to multipole (spin- 
multipole) matrix element //^j^, fjl)r, or /j^jj, respectively. 

Thus, we have completed a projection of the nuclear Hamiltonian into the space of phonon 
operators. Now we may assume that phonons obey boson statistics and work in the space of 
boson operators only. The presence of the term of interaction. Hint, in the model Hamiltonian 
means that the approximation, in which excited states of the nucleus are considered as pure 
one-, two-, multi-phonon states, is not sufficient. In fact, we have already mentioned above 
that it is not possible to describe the properties of the lowest collective vibrations in spherical 
nuclei in one-phonon approximation. It is also well-known that the coupling between one- and 
two-phonon configurations is the main mechanism for the damping of giant resonances. All 
this means that one needs to go beyond the approximation of independent phonons and take 
into account a coupling between them. To accomplish this task we write the wave function of 
excited states with angular momentum J and projection M in even-even nuclei in the most 
general form as a mixture of one-, two-, ■ ■ ■ phonon configurations: 

f D" ( J) 

^"(JM) = {j:s:,iJ)Qt, + E [Qt2QlhM (238) 
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^a3,/33,73 ^a3,33 ~^ ^03,73 ~^ '^/33,73 ~^ 2(5q,3 /33 5,^3^3 . (239) 



By greek characters we mean the phonon's identity, i.e. its multipolarity and order number, 
a = X'^i, the index z/ (= 1, 2, 3 . . .) labels whether a state J is the first, second, etc., state in 
the total energy spectrum of the system. It is assumed that any combination a, /3, 7 of phonons 
appears only once. The second and the third terms in Eq. ( |238| ) include phonons of different 
multipolarities and parities, they only must couple to the same total angular momentum J as 
the one-phonon term. 

Let us limit the wave function of excited states by three-phonon terms and diagonalize 
the model Hamiltonian of Eqs. ( p34|J236|) in the space of these states. We use for that a 
minimization procedure 

6 {< ^^(JM) I H I "^"{JM) > -E^ < "^^{JM) I "^"{JM) >} = , (240) 

which yields a set of linear equations over unknown wave function coefficients 5*^^(7), D'^^p^{J) 

{u^, - E;^)s:^{j) + y: k,pP)kip, = o , (241) 

02, /32 

ai 03/3373 

E D:,^M)uTJ:r. + (^"3 + ^/33 + ^73 - Ei)T::.^,^,^{j) = . 

02/32 

Applying boson commutation relations for phonons, the matrix element of interaction between 
two- and three-phonon configurations. 



^aSslS ~ \/l + '^02,/32\/l + '^03,/33,73 ([<5q2<5/32]jM I Hint. \ [Qa^Q^aQlsl J ^^"^ ' (^^2) 

can be expressed as a function of matrix elements of interaction between one- and two-phonon 
configurations. 



O2 = VI + ^02,/32(Qo. I H,nt. I [QtM^hM) = VI + 5o2,/32f/o1(«l) , (243) 

as following 
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and the value U^^{ai) is calculated according to Eq. ( p371) . Since we have used pure boson 
commutation relations for phonons the two-phonon configuration [02/32] j couples only to those 
three-phonon configurations [a3/5373]j where either as, or 73 are equal to a2 or ^2- This is 
governed by 5-functions in Eq. ( P44| ) . 

The number of linear equations ( 241 ) equals to the number of one-, two- and three-phonon 
configurations included in the wave function ( P38| ). Solving these equations we obtain the 
energy spectrum of excited states described by wave function ( |238| ) and the coefficients of 
wave function ( |241| ), S, D and T. 

It should be pointed out that within this approximation, in which phonons are considered 
as ideal bosons and nuclear Hamiltonian includes one-phonon exchange term, multi-phonon 
configurations of course posses no anharmonicity features. The strength of any one- or many- 
phonon configuration included in the wave function (|238|) fragments over some energy interval 



due to the interaction with other configurations. But the centroid of the strength distribution 
remains at the unperturbed energy. Thus, the energy centroid of two-phonon configuration 
[q;2/52] J equals exactly to the sum of energies of a2 and phonons for all values of J. To consider 
anharmonic properties of multi-phonon states one needs to go beyond pure boson features of 
excitations in even-even nuclei and take into account their internal fermion structure. 

Another reason to return back to the fermion origin of phonon excitations is two main 
problems in considering multi-phonon states associated with the boson mapping procedure. 
The first problem is an admixture of spurious npnh configurations which violate Pauli principle 
in the wave function of ra-phonon state. The second one is related to the fact that the set of 
pure n-phonon states is mathematically non-orthonormal if the internal fermion structure of 
phonons is taken into account (see Refs. 0, ^ for more details). To overcome these problems 



we will keep on using a phonon's imaging of nuclear excitation and use the same expression for 
the wave function of excited states ( |238| ) but in calculation of the norm of this wave function, 
< JM) I JM) >, and the energy of this state, < ^""{JM) \ H \ ^"{JM) >, we will use 
exact commutation relations between phonon operators: 

[Q\fii,Qy^'i']- = ^x^'^^l,^i'^i,i' - (^tm^^i'm' (245) 
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and exact commutation relations between phonon and quasiparticle operators: 

[oijmjQxni]- = Yi'jj'C^jmj'm''^t'm'y (246) 
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Also we will not expand the operator BT-{jj'; A/i) in Eq. ( p35| ) into a sum of phonon operators 



but use its exact fermion structure. The first term of Eq. (|245|) corresponds to the ideal boson 
approximation while the second one is a correction due to the fermion structure of phonon 
operators. The overlap matrix elements between different two-phonon configurations modifies 
as 

{[Qp'Qa']j\[QtQ^]j) = {[h'ba'UlKb^h) + K'{(3'a'\a(3) (247) 
where 6^ is the ideal boson operator and the quantity K, 

K^{/3'a'\a/3) = {X4iA3i3\XiiiX2i2) = ^{2Xi + 1)(2A2 + 1)(2A3 + 1)(2A4 + 1) 
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(248) 



is the Pauli principle correction coefficient. The experience of realistic calculations shows that 
usually \K-^ {[3a\a[3)\ ^ \K-^ {(3'a'\a[3)\ (where a ^ a' and/or [3 7^ (3') and that the so-called 
diagonal Pauli principle approximation, K-^ {l3'a'\a(3) = K-^ {a(3)da,a'^i3,i3', provides rather good 
accuracy and sufficiently simplifies the calculation. For these reasons we will use this diagonal 
Pauli principle approximation in what follows. 

The similar expression, as (|247|) , is valid for the overlap matrix elements between different 



three-phonon configurations. It can be used as a definition of the Pauli principle correction 



quantity Kj (7'/3'a'|a/?7) which we will also keep in diagonal approximation only. The relation 
between K-^^ajd) and Kj{aP'y) quantities is the following ||9^ : 



Kf{ap-f) =K\al3) ( 3 + ^ f/2(a/?J7; /, /')K^'(/37) 



(249) 



where f/^ stands for the Jahn coefficients . 

When internal fermion structure of phonons is taken into account and exact commutation 
relations ( |245| , |246D are applied the secular equation ( |241| ) transforms into 



(250) 
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+ Acj^g^g + Au;;^^^^ are anharmonicity shifts of 



two- and three-phonon configurations, respectively, due to the Pauli principle corrections. In 
diagonal approximation they can be calculated according to: 
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Another role of Pauli principle corrections is a somewhat renormalization of the interaction 
between n- and {n + l)-phonon configurations. We have used the same notations for these 

as in the case of the 'ideal boson approximation' (see, Eqs. (|242| , |243|) ). 



matrix elements U^l^^ 



But calculating the matrix elements {Qa^ \ Hint. \ [QajQ/jj]-^*'^) ^ske into account the fermion 
structure of phonons and nuclear Hamiltonian and obtain 
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(252) 



where the value U^'^{ai) is calculated again according to Eq. ( |237|) . A similar additional factor 



[1 + 1/2 X Kj (a^Ps'js)] receives the matrix element of interaction between two- and three-phonon 
configurations. 

The minimal value of the quantity K"^ {a(3) equals to —2. It corresponds to the case of the 
maximal Pauli principle violation, i.e. to a spurious multi-phonon configuration. It happens 
only when a2 and (^2 phonons are purely two-quasiparticle states. In such a case the matrix 
element of interaction U'^lp.^ = (see, Eq. ( |252| )) and the spurious state is completely separated 
from other states. While dealing with collective a2 and (32 phonons, when a possible admixture 
of the spurious four-quasiparticle configurations is small, the value of {a(3) is close to 0. 
Nevertheless the value of the anharmonicity shift Acu;^ ^ is not vanishingly small for the later 



because of the relatively large value of the ratio /y^ in Eq. ( |251| ). This shift is the largest 
one for the collective low-lying multi-phonon configurations. For non-collective multi-phonon 
states the shift is small because of the small value of the above mentioned ratio. 

Equations ( p50| ) have been obtained under two main assumptions. The first one is the 
already discussed diagonal Pauli principle approximation. The second assumption is the ne- 
glecting of the higher order terms of the interaction part of the nuclear Hamiltonian as compared 



to the one in Eq. ( |236| ) which couples n- and [n ± l)-phonon configurations. For Eqs. (|250|) it 



means that a direct couphng between one- and three-phonon configurations of the wave func- 
tion ( P38| ) which is possible due to non-zero matrix element {Qai I c^jm'^jm \ [Qt^QXQ'ysljM) ^ 
is neglected. In realistic calculation we will also use a selection of three-phonon configurations 
provided by Eq. ( |244| ) although now the matrix element U^lplyj^ 7^ even if one of 0:3, or 73 
is not necessarily equal to 02 or (32- These omitted matrix elements are orders of magnitude 
smaller as compared to the accounted for ones. 

Solving the system of linear equations (|25CI|) we obtain the spectrum of excited states, E^, 



described by the wave function ( p38| ) and coefficients 5'^^(J), D^^f^^i^J) and T^^p^^^{J) refiecting 
the phonon structure of excited states. Usually, in calculation of the properties of single giant 
resonances the three-phonon terms of the wave function (|238|) are omitted. Then it is possible 



to solve the system of linear equations ( |250|) with the rank of the 10^ — 10^ order by a direct 
diagonalization. But while considering the damping properties of two-phonon resonances, three- 
phonon configurations cannot be omitted. For this case instead of the diagonalization of the 
linear matrices of very high orders, an alternative solution is possible. We may substitute the 
first and last equations of ( |25(J| ) into the second equation and obtain the system of non-linear 
equations 

; / , , A 7 7-i7\r "2/32 aLfiL 

det 1 1 [uja^ + up^ + Aw„^^^ - E^ ) 5^2/32, - 
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, (253) 



03/3373 + ^/33 + ^73 + ^^4/3373 - Ei 

the rank of which equals to the number of two-phonon configurations included in the wave 
function ( p38| ). The solution of the system ( |253| ) by some iterative method yield again the 
spectrum of excited states E-^ and coefficients D^^f^^{.J). Other coefficients of the wave function 
( P38| ) are related to these coefficients as following 

QV ( 1\ _ ^02/?2 -Po2/?2('^) ^02/32 

V ( 1\ r7"2/32 

rpu ( J\ _ _ ^02/32 -^02/^2 V*^ / ^03/3373 (O^A\ 

^03/3373Wj - , , I , I , , I A, ,J pJ • 

t^03 + t^/33 + ^^73 + ^'^03/3373 ~ 

It may be argued that the boson mapping with keeping the fermion information of the 
phonons' images at all stages of transformations gives no advantage as compared to npnh ap- 
proach since, mathematically, a direct correspondence between two methods can be established 
only if the full basis of n-phonon states is used. However many npnh configurations interact 
very weakly with other ones and as a result practically do not mix with them. It allows a 



sufficient truncation of multi-phonon configurations in the wave function (|238|) based on their 
physical properties with keeping a good accuracy for the components important for the subject 
of research. From the point of view of the Pauli principle violation the most dangerous multi- 
phonon configurations are the ones made of non- collective RPA states. On the other hand, 
these configurations interact with the other ones much weaker than the multi-phonon configu- 
rations including at least one collective phonon. For these reasons the first are not accounted 
for in the wave function (|238|) in realistic calculation. As the criteria "collective/non-collective" 



we take the contribution of the main two-quasiparticle component to the wave function of the 
phonon operator. If the contribution exceeds 50-60% we will call the phonon non-collective. 



Let us consider now the electromagnetic excitation of pure one- and multi-phonon states 
from the ground state. The one-body operator of electromagnetic transition has the form 



".P ^ i\\F\\\i' ^ 

M{E\y.) = E el'^ E ^-^y'"''^ % C}:,^'^U^r-m' (255) 
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where the single particle transition matrix element < j||EA||j' >=< j| |j' > and e^^^ are 

effective charges for neutrons and protons. In calculations we use the following values of effective 
charges: e^^^ = —Z/A and e^^ = N/A to separate the center of mass motion and e^^^^ = and 
e^^^^^ = 1. Performing the transformation from particle operators to quasiparticle and phonon 
ones in Eq. ( [^55| ), this equation transforms into 



E C^::.fm'{-y'-"''c^t'm'C^r-m' \ (256) 
mm' 



where the first term corresponds to one-phonon exchange between initial and final states and 
the second one is responsible for "boson forbidden" electromagnetic transitions (see, for details 
Ref. 0]). Then the reduced matrix element of the electromagnetic excitation of the one-phonon 
state Xi from the ground state 0^^ in even-even nuclei may be calculated according to 



n,p 1 

< QM\\MiEX)\\0^,,, >= E4'^ E 2 < J^\\EM\J2 > (^^,, + ^i:,^ . (257) 
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Due to the ground state correlations the direct excitation of pure two-phonon states [Qx^i-^ x 
Qx2i2\\ from the ground state is also possible when we are dealing with the RPA phonons. The 
physical reason for that becomes clear if we remember that the ground state wave function 
includes a small admixture of four-, eight-, ■ ■ ■ quasiparticle configurations (see, Eq. (|232|) ). 
The second term of Eq. (|256|) is responsible for these transitions and the reduced matrix element 
can be obtained by applying the commutation relations (|246|) . It has the form 

< [Qx,^, X gA,n]A II ^(i^A)||0+,. >= v/(2Ai + l)(2A2 + l)E4'^ E ^^"^ 



J1J2 



X <jiii^Aib,>| • (258) 

Another type of boson forbidden 7-transitions which take place due to the internal fermion 
structure of phonons are the ones between one-phonon initial, Qt^^i-^Dph, and final, (5^242 1)?^' 
states. The reduced matrix element of such transitions can be calculated according to 

n.p 

,(-) 



<g,,,l|A^(EA)||g+,^> = v^2A2 + iE4'^ E ^U<^ill^^lb2> 

jlj2j3 

The matrix element for transitions between the two-phonon states [Q^^iii ^ Qt2i2\^'\)ph 
[Qx'^i'i X QxdiA]>^"\)ph is very complex and not presented here. Its first order term is very similar 



to the one for transitions between the one-phonon states Q^inDp^ Qt^iiDph and may be 
obtained by assuming that the fermion structure of one phonon is "frozen", i.e., assuming that 
A2«2 = Asis. 

When the couphng between one- and multi-phonon configurations is accounted for in the 
wave function of excited states, the reduced matrix element of the electromagnetic excitation 
of the states of Eq. (|238D may be written as 



<vl/'^(J)||A<(EA)||0+,.> = \Y.S:,iJ)<QM\\MiEX)\\Ol,> (260) 
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where we have neglected the direct excitation of three-phonon configurations from the ground 
state. Since an admixture of multi-quasiparticle configurations in the ground state wave func- 



tion is very small, the reduced matrix element Eq. ( p58| ) is typically about two orders of mag- 
nitude smaller as compared to the reduced matrix element Eq. ( |257|) . For this reason, in most 
of the cases keeping only the first term in Eq. (|260|) and neglecting the second one together 



with interference effects provides very good accuracy in calculation. Nevertheless, there are a 
few exceptional cases. 

The first one is the excitation of the lowest 1~ state in spherical nuclei. It is well known 
that no collective one-phonon 1~ configurations appear in the low energy region and the wave 
function of the Ij" state has the dominant two-phonon component [2^ x Sfji-- There are three 
main mechanisms to explain the i?l-excitation of this state observed in the experiment [p5[] . 
The first is an influence of the GDR. In microscopic theories it appears in a natural way due to 
the coupling of one- and two-phonon configurations. Since the GDR is located about 10 MeV 
higher, this coupling yields only a very small portion of the observed strength. The second 
mechanism is the excitation of non- and weakly-collective one-phonon 1~ configurations which 
have relatively small B{E1) values but are located in low energy region. The last mechanism is 
the direct excitation of two-phonon configurations from the ground state. Although the direct 
excitation of two-phonon configurations from the ground state is a second order effect, excitation 
of collective two-phonon configurations [2^ x 3^]i- play an essential role since the other two 
mechanisms yield much weaker El strengths. In this case, interference effects between the first 
and the third mechanisms are also important |M . 



The second term of Eq. ( |260| ), although very weak as compared to the first one, may also play 
some role at the excitation energies above 20 MeV where the density two-phonon configurations 
is a few orders of magnitude higher as compared to the density of one-phonon configurations. 
It will be discussed below. 

Considering the two-step mechanism of the DGDR excitation in second order perturba- 
tion theory we also need the reduced matrix element of the electromagnetic excitation of the 
two-phonon DGDR state [1~ x from the one-phonon GDR state 1~. In ideal boson ap- 
proximation this matrix element 



< [1.^ X l-]j\\MiEl)\\lr >= ^(1 + 5,,)^^^ < l,:;||A^(i?l)||^... > . (261) 
4.3 Comparison with other approaches 



The properties of the double giant resonances have been also microscopically studied with the 
Skyrme forces 171, and within the second- RPA approach ETl, p8|. 



The most close to the QPM approach is the one of the first group of papers. The main 
difference between these two approaches is that in calculations with the Skyrme forces the 
properties of the ground and Iplh excited states are calculated self-consistently. As within 
the QPM, in calculations with the Skyrme forces the Iplh basis is mapped into the phonon 
space. Multi-phonon states are obtained by folding of one-phonon states. The phonon basis in 
Refs. JT^ 1^ is restricted by only a few, the most collective, phonons for each multipolarity. 
Calculations are performed with the wave function including one- and two-phonon terms. The 
main attention is paid to the effects of anharmonicity and non-linearity. The later is an influence 
of taking into account of boson forbidden transition matrix elements, Eqs. ( p58|J259|) , on the 
absolute value of the DGDR excitation in heavy ion collisions. 

In the second-RPA approach [|99| the wave function of excited states is written as a mixture 
of Ip— Ih and 2p— 2h configurations: 

Qt\)g-s. = y^("^pfe'^p Q-fe ~ ^pfeO-feO-p) + {^pp'hh'(^v(^p'(^h'0-h — Ypp'hh'^h^h'^v'^p)\) a s.- 

ph pp'hh' 

(262) 

The operators Q'^ are assumed as bosons and the energy spectrum and coefficients X and Y 
are obtained by diagonalization of the model Hamiltonian in the space of states described by 
the wave functions of Eq. (|262|). 



5 Physical properties of the double giant resonances 

In the present chapter we will consider the properties of the DGDR as predicted by the QPM 
mainly in ^^^Xe and ^osp]^ ^^j. ^j^ich experimental data in relativistic heavy ion collision (RHIC) 
are available. Before proceeding with that let us briefly check an accuracy of the description of 
the properties of low-lying states and single giant resonances within this approach. It provides 
an estimate how good the phonon basis, to be used in the forthcoming calculation of the DGDR 
properties, is since no extra free parameters are used after this basis is fixed. The results of our 
calculations of the position and exhaust of the energy weighted sum rule (EWSR) of low-lying 
states and giant resonances as well as the width of resonances in ^^^Xe and ^osp]^ g^j.^ presented 
in Table |TU| in comparison with the experimental findings. The comparison indicate a rather 
good correspondence between calculated characteristics and experimental data. The calculation 
somewhat underestimate the width of resonances, especially of the isovector GQR. The main 
reason is related to the necessity of truncating of complex configurations included in the wave 
function of excited states in actual calculation. The density of multi-phonon configurations is 
rapidly increasing with the excitation energy. That is why the effect of the basis truncation the 
most strongly influences on the width of the GQRj^ located at higher energies. 

5.1 One-step excitation of two-phonon states in the energy region 
of giant resonances 

Let us consider a direct photoexcitation of the two-phonon states in the energy region of giant 
resonances from the ground state of even-even nuclei (see, Refs. ||57| , |10CI|] for more details). 
Since in RHIC experiments the Coulomb mechanism of excitation plays the most essential 
role, the cross sections of photoexcitation can be easily recalculated into RHIC cross sections 
for different energies and Z- values of target and projectile nuclei. In calculation of the BiEX) 



Table 10: Integral characteristics (position, Ex, exhaust of the energy weighted sum rule 
(EWSR) and width of resonances, F) of low-lying excited states and one-phonon giant res- 
onances in ^^^Xe and ^''^Pb. 
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values we use only the terms proportional to (see, Eq. ( p58| )). The complete set of diagrams 
corresponding to a direct transition to two-phonon states from the ground state is presented in 
Ref. | |101| . As one can see from the analytical expressions the main part of the contributions 
from different terms disappears due to the cancellation between particles and holes. 

The cross sections of the direct photoexcitation of the groups of two-phonon states made of 
phonons of definite multipolarities in ^^^Xe and ^o^Pb are presented in Fig. £'2-excitation of 
[1~(8)1~]2+ states is plotted in the top part of the figure. i?l-excitation of the two-phonon states 
[1~ (g) 2+]i- and [2+ ® 3^]i- is shown in the middle and the bottom parts, respectively. The 
integral characteristics of two-phonon states which are a single giant resonances built on top of 
either a low-lying state or another single resonance in the same nuclei are given in Table [T^. 

The main feature of the top part of Fig. ^ is that just all two-phonon states which form 
this double-phonon resonance are constructed of the one-phonon 1~ states belonging to the 
GDR in the one-phonon approximation. The structure of the [l~(8>2+]i- and [2"'"®3~]i- states 
is more complex. For example, among [1^ ® states the substructure in the energy range 

from 15 to 20 MeV in ^osp]^ (right middle part of Fig. |28|) is formed mainly by 1^^ phonons from 
the GDR region coupled to the 2f state. The small substructure above 32 MeV is due to the 
GDR 1~ phonons coupled to the 2^ phonons of the isovector GQR. As for the broad structure 
between 20 and 30 MeV not only [GDR® GQRj^]i- states but many other two-phonon states 
built of less collective 1~ and 2^ phonons, the role of which is marginal for properties of single 
resonances, play an essential role. The same conclusions are valid for the direct photoexcitation 
of [1~ ® 2"'"]i- states in ^^^Xe. The cross section of the direct photoexcitation of the two-phonon 
1~ states built of phonons of the higher multipolarities yield non-resonance feature. It is already 
seen for the case of [2+ (g) 3~]i- states (bottom part of Fig. |28|), especially in ^°^Pb. 

While dealing with electromagnetic, or with Coulomb excitation from a O"*" ground state, the 



Table 11: Integral characteristics (energy centroid, width and cross section of direct photoexci- 
tation from the ground state) of some groups of two-phonon states which are a giant resonance 
built on top of either a low-lying state or another single giant resonance in ^^^Xe and ^''^Pb. 



Nucl. Configuration Centroid Width o"^ 

[MeV] [MeV] [mb] 





24.0 


2.9 


4.3 


'''Xe [l-,,^„§5 


1gdr,J2+ 


30.2 


4.0 


0.33 




2gqr,:J2+ 


21.3 


0.5 


0.1 






17.4 


2.2 


1.7 


[■^GDRj„ 


® 2th- 


17.2 


2.1 


8.7-10-^ 




®2+]3- 


17.7 


3.4 


4.9 ■ 10"^ 


[-'-GDRii, 




15.3 


3.9 


5.2 ■ 10-2 


[1gdr,„ '8 


^GQRjJl- 


25.1 


3.8 


9.6 


[1gdr,„ ® 


1gdr,J2+ 


25.5 


4.4 


0.22 



priority attention has to be paid to the final states with the total angular momentum and parity 
J'^ = 1-. For that we have calculated the cross section for the photoexcitation of two-phonon 
states [A^^ ® ^2^]i-, where A^^ and A2^ are both natural A'^" (vr" = (—1)'^) and unnatural A'^" 
(vr" = (— 1)'^+^) parity phonons with multipolarity A from to 9. 

The results of the calculation for ^^^Xe and ^osp^ integrated over the energy interval from 
20 to 35 MeV are presented in Table |T2l Each configuration [A^^ (S> A2^] in the table means a 
sum over a plenty of two-phonon states made of phonons with a given spin and parity A^ ^ , A2^ , 
but different RPA root numbers ii, 12 of its constituents 

a([Ar ® A^^]) = ® . (263) 

11,12 

The total number of two-phonon 1^ states included in this calculation for each nucleus is about 
10^ and they exhaust 25% and 15% of the EWSR in ^^^Xe and ^ospb, respectively. The absolute 
value of the photoexcitation of any two-phonon state under consideration is negligibly small 



but altogether they produce a sizable cross section. Table demonstrates that different two- 
phonon configurations give comparable contributions to the total cross section which decreases 
only for very high spins because of the lower densities of such states. As a rule, unnatural parity 
phonons play a less important role than natural parity ones. For these reasons we presented 
in the table only the sums for [natural®unnatural] and [unnatural(8)unnatural] two-phonon 
configurations. 

The cross section for the photoexcitation of all two-phonon 1" states in the energy region 
20-35 MeV from the ground state equals in our calculation to 511 mb and 423 mb for ^'^^Xe 
and ^ospb^ respectively. It is not surprising that we got a larger value for ^^^Xe than for 
208p|-| xhis is because the phonon states in Xe are composed of a larger number of two- 
quasiparticle configurations due to the pairing. The same values for two-phonon states with 
angular momentum and parity J'^ = 2"*" are an order of magnitude smaller. We point out that 
the direct excitation of [1~®1"]2+ or [GDR®GDR]2+ configurations is negligibly weak (compare 
results in Tables |lT] and 0). The calculated values should be compared to the cross section 
for the photoexcitation of the single-phonon GDR which in our calculation equals to 2006 mb 



Table 12: Cross sections for the direct photoexcitation of different two-phonon configurations 
from the ground state integrated over the energy interval from 20 to 35 MeV in ^^^Xe and 
^°^Pb. The GDR cross section integrated over the energy of its location is presented in the last 
line for a comparison. 







U y 


lliU 


Configuration 




208 


[0+ ® l~ 


1 


4.4 


3.9 


[1- ®2+ 


1- 


36.6 


44.8 


[2+ ®3- 


1- 


82.8 


33.1 


[3- ® 4+ 


1- 


101.0 


56.7 


[4+ ® 5- 


1- 


68.9 


37.3 


[5- ® 6+ 


1- 


49.2 


46.2 


[6+ ® 7- 


1- 


31.9 


49.8 


[7- ® 8+ 


1- 


13.6 


12.5 


[8+ ® 9- 


1- 


4.9 


9.0 


Ai,A2=l 
Ai,A2=l 


1- 


1 1.4 


oo.o 


1- 


46.7 


71.1 


9 n,u n,u 

E [Xi' ®Xl' 

Ai,A2=0 


1- 


511.4 


422.9 


[GDR ® GDR 

E [Ai^ ® A2^ 

Ai,A2=l 


2+ 


0.33 


0.22 


2+ 


38.1 


21.7 


GDR 




2006 


2790 



and 2790 mb, respectively. A contribution of two-phonon 1^ states to the total cross section 
at GDR energies is weaker than at higher energies because of the lower density of two-phonon 
states and the lower excitation energy and can be neglected considering the GDR itself. It is 
clearly demonstrated in Fig. pp|b. In this figure the cross sections of the photoexcitation of 1~ 
states in ^^^Xe and ^osp]^ presented. The top part of the figure corresponds to a calculation 
performed in one-phonon approximation. The results of calculations with the wave function 
which includes a coupling between one- and two-phonon 1~ configurations are plotted in the 
bottom part of the figure. For a visuality the last calculations are also presented as strength 
functions 

= (264) 

with a smearing parameter A = 1 MeV, where al is a partial cross section for the state with 
the excitation energy plotted also by a vertical line. The i?l-transitions to one-phonon 
components of the wave function of excited 1~ states are plotted by dashed curve. It should 
be compared with the solid curve which is the sum of transitions to one- and two-phonon 1~ 
configurations in the GDR energy region. 

For ^°'^Pb photoexcitation cross sections are known from experimental studies in (7, n) 
reactions up to the excitation energy about 25 MeV ||102| , |103|| . It was shown that QPM provides 
a very good description of the experimental data in the GDR region [p.02|| (see, Fig. ROl), while 



theoretical calculations at higher excitation energies which account for contributions from the 



single-phonon GDR and GQRj^ essentially underestimated the experimental cross section ||103|] . 
The experimental cross sections above 17 MeV are shown in Fig. ^ together with theoretical 
predictions. The results of the calculations are presented as strength functions obtained with 
averaging parameter equal to 1 MeV. The contribution to the total cross section of the GQRj^ 
(short-dashed curve), the high energy tail of GDR (long-dashed curve), and their sum (squared 
curve), are taken from Ref. [|103|1 . The curve with triangles represents the contribution of the 
direct excitation of the two-phonon states from our present studies. The two-phonon states 
form practically a fiat background in the whole energy region under consideration. Summing 
together the photoexcitation cross sections of all one- and two-phonon states we get a solid 
curve which is in a very good agreement with the experimental data. 

Thus, from our investigation of photoexcitation cross sections we conclude that in this 
reaction very many different two-phonon states above the GDR contribute on a comparable 
level, forming altogether a flat physical background which should be taken into account in the 
description of experimental data. On the other hand. Coulomb excitation in relativistic heavy 
ion collisions provides a unique opportunity to excite a very selected number of two-phonon 
states by the absorption of two virtual 7's in a single process of projectile-target interaction 
0. Theoretically this process is described using a the second order perturbation theory of the 
semi-classical approach of A. Winther and K. Alder ^ and discussed in Sect. 3.2.3. Since 



excitation cross sections to second order are much weaker than to first order of the theory, 
two-phonon states connected to the ground states by two i?l-transitions are predominantly 
excited. These two-phonon states have the structure [!"(«) ® and form the DGDR. 

5.2 1+ component of the DGDR 

According to the rules of angular momentum coupling two one-phonon states with the spin 
and parity equal to 1~ may couple to the total angular momentum J'" = 0"*", 1"*" and 2+. Thus, 
in principle, three components of the DGDR with these quantum numbers should exist. In 
phenomenological approaches describing the single GDR as one collective state, the [1^ (S> 
component the DGDR is forbidden by symmetry properties. Taking into account the Landau 
damping this collective state splits into a set of different Ij" states distributed over an energy 
interval. In microscopic studies the Landau damping is taken into account by solving the RPA 
equations. Again, the diagonal components [1^ ® ^ire forbidden by the same symmetry 

properties but nondiagonal ones [1~ exist and should be taken into consideration. 

Consequently, the role of these nondiagonal components depends on how strong is the Landau 
damping. 

We produce here two-phonon DGDR states with quantum numbers = 0"^, 1^ and 2+ by 
coupling one-phonon RPA states with the wave function |1~ >„, to each other. The index m 
stands for different magnetic substates. The wave function of the two-phonon states has the 
form: ^ 

|[1- ® lr]j-=o+,2+ >A/= ^ E {lrnlm'\JM)\l- >m iK >m' , (265) 

V m,m' 

for two-phonon states made of two identical phonons while for other DGDR states it is: 

® l,v]j.=o+,i+,2+ >M= E {lrnlm'\JM)\l- >^ . (266) 

m,m' 

In the present calculation we do not include the interaction between DGDR states, of Eqs. 
(P65|J266|), and we do not couple them to states with different than two number of phonons (it 



will be considered below). Thus, our two-phonon states |[1~ ® lj~]jvr >m have excitation energy 
equal to the sum of one-phonon energies uoi + uoii and are degenerated for different values of the 
total spin J'^ and its projection M. 

Since the main mechanism of excitation in projectile ions at relativistic energies is the 
Coulomb part of interaction with a target, the nuclear part of interaction has been neglected 
in the present analysis. In a semi-classical approach |^ the two-phonon DGDR states can be 
excited in second-order perturbation theory via the two-step process g.s. — > GDR — > DGDR. 
The second order amplitude can be written as 

where assuming the Coulomb mechanism of excitation the first order amplitude a[j^ is 
proportional to the reduced matrix element of < J/||£^l||Jj >. The reduced matrix element 
< [lj7 (S) Ij'jjTrlli^^llllj^ > of electromagnetic excitation of two-phonon states, Eqs. ( p65| , p66| ), 
from the one-phonon state |1~ >m is related, in the boson picture of nuclear excitation, to 
the excitation of |1~ from the ground state according to (261). It should be noted that 
although for the two-phonon states, Eq. (|265|) , we have an extra factor \/2, the states of Eq. 
( P66| ) play a more important role in two-step excitations since they can be reached by two 
different possibilities: g.s. 1~ [1~ ® and g.s. 1^ [1^ ® 1^]. 

First of all, we point out that in second-order perturbation theory the amplitude for this 
process is identically zero in a semi-classical approach. This can be understood by looking 
at Fig. The time-dependent field Vei carries angular momentum with projections m = 
0, ±1. Thus, to reach the 1"*" DGDR magnetic substates, many routes are possible. The lines 
represent transitions caused by the different projections of Vei'. (a) dashed lines are for m = 0, 
(b) dashed-dotted lines are for m = —1, and (c) solid-lines are for m = +1. The relation 
VEi^m=o 7^ VEi,m=±i holds, SO that not all routes yield the same excitation amplitude. Since the 
phases of the wave functions of each set of magnetic substates are equal, the difference between 
the transition amplitudes to a final M, can also arise from different values of the Clebsch-Gordan 
coefficients (lmlm'|lM). It is easy to see that, for any route to a final M, the second-order 
amplitude will be proportional to (001m|lm)(lmlm'|lM) Vei^th' VEi,m + {m ^ m'). The 
two amplitudes carry opposite signs from the value of the Clebsch-Gordan coefficients. Since 
(001m|lm) = 1, the identically zero result for the excitation amplitude of the 1+ DGDR state 
is therefore a consequence of 

^(Imlm'IlM) = . (268) 



We have also performed a coupled-channels calculation ||104|| following the theory described 
in Ref. As shown in Ref. [|^, the coupling of the electric quadrupole (isovector and 

isoscalar) and the electric dipole states is very weak and can be neglected. We therefore include 
in our space only one-phonon 1~ and two-phonon [1^^ ® '^i'].]^ {J^ = 0^, 1^ and 2+) states. 
In coupled-channels calculation we take into account interference effects in the excitation of 
different GDR and DGDR states and obtain the occupation amplitudes by solving the coupled- 
channels equations. By solving these equations we thus account for unitarity and for multi-step 
excitations, beyond the two-step processes of Eq. ( ^67| ). The time dependent electric dipole 
field is that of a straight-line moving particle with charge Ze, and impact parameter h (we use 
Eqs. (25-26) of Ref. pf). 

Due to the large number of degenerate magnetic substates, to make our coupled-channels 
calculation feasible, we have chosen a limited set of GDR and DGDR states. We have taken 



six 1~ states which have the largest value of the reduced matrix element < >. 
These six states exhaust 90.6% of the classical EWSR, while all 1~ states up to 25 MeV in 
our RPA calculation exhaust 94.3% of it. This value is somewhat smaller than the 122% 
reported in Ref. ||3^. It is because the continuum in our RPA calculation was approximated 
by narrow quasibound states. From these six one-phonon 1~ states we construct two-phonon 
[1~ lj~]jvr states, Eqs. ( |265| , |266D , which also have the largest matrix element of excitation 
< [IJ' ® jvr I |1^ > for excitations starting from one-phonon states 0. The number of 
two-phonon states equals to twenty one for J'^ = 0^ and 2+, and to fifteen for J'^ = 1+. 
The cross section for the DGDR excitation was obtained by summing over the final magnetic 
substates of the square of the occupation amplitudes and, finally, by an integration over impact 
parameter. We have chosen the minimum impact parameter, b = 15.54 fm, corresponding to 
the parameterization of Ref. appropriate for lead- lead collisions. 

The electromagnetic excitation cross sections for the reaction ^°^Pb (640A MeV) -l-^'^'^Pb 
with excitation of all our basic 63 states is shown in Fig. BBI. The total cross sections for 
each multipolarity are presented in Table [T^, together with the results of first-order (for one- 
phonon excitations) and second-order (for two-phonon excitations) perturbation theory. The 
coupled-channels calculation yields a non-zero cross section for the 1"*" DGDR state due to other 
possible routes (higher-order), not included in second-order perturbation theory. One observes 
a considerable reduction of the DGDR cross sections, as compared to the predictions of the 
second-order perturbation theory. The GDR cross sections are also reduced in magnitude. 
However, the population of the 1^ DGDR states are not appreciable and cannot be the source 
of the missing excitation cross section needed to explain the experiments. In general, the 
coupled-channels calculation practically does not change the relative contribution of different 
one-phonon 1^ and two-phonon states [1^ ® Ijv] to the total cross section with given J'^ = 1~ , 
0"*" and 2+. But since the l"*" component of the DGDR, with its zero value of excitation 
cross section in second-order perturbation theory, has a special place among the two other 
components, the main effect of coupled-channels is to redistribute the total cross section between 
the J'^ = O"*", and J'^ = components. 

The calculated cross section in coupled-channels for both GDR and DGDR are somewhat 
smaller than reported in experimental findings [^4|, |105|| . This is not surprising since as men- 



tioned above our chosen six 1~ states exhaust only 90.6% of EWSR while the photo-neutron 
data 1^ indicate that this value equals to 122%. Due to this underestimate of exhaust of 
the EWSR the cross section for the DGDR excitation reduces more strongly than the one for 
the single GDR. This is because the GDR cross section is roughly proportional to the total 
B[E1) value while for the DGDR it is proportional to the square of it. We will return back in 
more details to the problem of absolute cross sections of the DGDR excitation in RHIC in the 
forthcoming subsection. 



^As demonstrated in the previous subsection the direct excitation of two-phonon configurations from the 
ground state is very weak. It allows us to exclude in our calculation matrix elements of the form < [1~ ® 
l~]2+(i+) I |-E2(M1)| l^.s. > which correspond to direct transitions and produce higher order effects in comparison 
with accounted ones. These matrix elements give rise to DGDR excitation in first order perturbation theory. 
Thus, to prove our approximation we have calculated such cross sections and got total values equal to 0.11 mb 
and <0.01 mb for the twenty one 2+ and the fifteen 1+ basic two-phonon states, respectively. These values have 
to be compared to 244.9 mb for the total DGDR cross section in the second order perturbation theory. 



Table 13: Cross section (in mb) for the excitation of the GDR and the three components 
with r = 0+,2+,l+ of the DGDR in ^ospb (640A MeV) +208pb colhsions. Calculations 
are performed within coupled-channels (CC) and within first (PT-1) and second (PT-2) order 
perturbation theory, respectively. 





CC 


PT-1 


PT-2 


GDR 


2830. 


3275. 


0. 


DGDRo+ 


33.0 


0. 


43.1 


DGDR2+ 


163.0 


0.11 


201.8 


DGDR1+ 


6.3 


<0.01 


0. 


DGDR/GDR 


0.071 


0.075 





5.3 Position, width and cross section of excitation in RHIC of the 
DGDR in ^^^Xe and ^o^pb 

To describe the width of two-phonon resonances it is necessary to take into account a coupling 
of two-phonon configurations, which form these resonances, with more complex ones. For that 



two types of calculations have been performed. In the first of them | 106 | the fine structure of 
the GDR calculated with the wave function which includes one- and two-phonon configurations 
and presented in Fig. p9| b has been used. The DGDR states have been constructed as a product 
of the GDR to itself. In other words, following the Axel-Brink hypotheses on top of each 1~ 
state in Fig. |29[b we have built the full set of the same 1~ states. 



The calculation has been performed for the nucleus ^^^Xe. In the dipole case, A'^ = 1 , the 
one-phonon states exhaust 107% of the classical oscillator strength and are displayed in the 



left part of Fig. ^a. Of these, 20 states have an oscillator strength which is at least 1% of 
the strongest strength and together exhaust 104% of the classical EWSR. We have used these 
states in the coupling to two-phonon states. We have included all the natural parity phonons 
A'^ = 1^ — 8+ with energy lower or equal to 21 MeV, obtaining 2632 two-phonon configurations. 
One obtains 1614 states described by the wave function which includes one- and two-phonon 
configurations, in the energy interval from 7 MeV to 19.5 MeV. Their photoexcitation cross 
sections are shown in Fig. |29|b. The B{E1) value associated with each mixed state is calculated 
through its admixture with one-phonon states, as | < i/||A1(i?l)||0 > p = \J2iSi{l~) < 
0\\Qi-iM.{El)\\0 > p. Also shown by dashed curve in the left part of Fig. p9| b is the result 
obtained adding an averaging parameter of 1.0 MeV. This parameter represents in some average 
way the coupling to increasingly more complicated states and eventually to the compound 
nuclear states. From the resulting smooth response it is easy to directly extract the centroid 
and the full width at half maximum of the GDR. The corresponding values are -Egdr = 15.1 
MeV and Tqj^^ = 4 MeV. They can be compared with the values extracted from experiment, 
EaoR= 15.2 MeV and Fgdr = 4.8 MeV. 

The isoscalar and the isovector giant quadrupole resonances (GQR) have also been calcu- 
lated. The centroid, width and percentage of the EWSR associated with the isoscalar mode 
are 12.5 MeV, 3.2 MeV and 75% respectively. The corresponding quantities associated with 
the isovector GQR are 23.1 MeV, 3.6 MeV and 80%. 

The differential Coulomb-excitation cross sections as a function of the energy associated 
with the one-phonon GDR and GQR resonances and the two-phonon DGDR in ^^^Xe (690A 



Table 14: Calculated (for two values of Tq) and experimental cross section (in mb) for the 
excitation of giant resonances in ^^^Xe in ^^^Xe (690A MeV) + ^ospb reaction. In the last row, 
the experimental cross sections for Coulomb excitation of one- and two-phonon states from 
Ref. are shown. The value of the integrated cross section reported in Ref. ||4^ is 1.85 ±0.1 
b. The nuclear contribution has been estimated in Ref. to be about 100 mb, while about 
3% (50 mb) of the cross section is found at higher energy. Subtracting these two contributions 
and the two-phonon cross section, leads to the value 1485 ±100 mb shown in the Table. 





GDR 


GQRjs 


GQRto 


GDR ± GQR 


DGDR 


To = 1 , 2 f m 


2180 


170 


120 


2470 
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To = 1, 5 fm 


1480 


110 


60 


1650 


50 


Experiment 


1024±100 






1485±100 


215 ± 50 



MeV) ± ^°^Pb reaction are displayed in Fig. 133. It is seen that the centroid of the two-phonon 
dipole excitation falls at 30.6 MeV, about twice that of the one-phonon states, while the width 
is r ^ 6 MeV, the ratio to that of the one-phonon excitation being 1.5. 

The associated integrated values are displayed in Table |l^, in comparison with the exper- 
imental findings. The cross sections depend strongly on the choice of the value of 6mm = 
ro{Ap^^ + A]^^). In keeping with the standard "safe distance", that is, the distance beyond 
which nuclear excitation can be safely neglected, we have used = 1.5 fm. Because their val- 
ues essentially do not depend on the width of the GDR, we view the calculated cross section of 
1650 mb as a rather accurate value and if anything an upper limit for the one-phonon Coulomb 
excitation cross section. It is satisfactory that the measured cross section is rather close to this 
value. Also shown in Table H are the predictions associated with the sequential excitation of 
the DGDR. This result is essentially not modified evaluating the direct Coulomb excitation of 
the double GDR. In fact, the cross section associated with this process is a factor 10~^ smaller 
than that associated with the two-step process. The calculated value of 50 mb is a factor of 
0.25 smaller than experimentally observed. 

Two other processes are possible within the sequential excitation of the giant modes which 
can lead to an excitation energy similar to that of the two-phonon GDR. They are the excitation 
of the isoscalar GQR mode followed by a GDR mode and vice versa. The resulting cross 
section is estimated to be an order of magnitude smaller, cfr. Table [T^, and does not change 
qualitatively this result. In order to make clearer the seriousness of this discrepancy, we have 
recalculated all the cross sections using = 1.2 fm, namely with a much smaller radius than 
that prescribed in order to respect the safe Coulomb excitation distance of closest approach. 
The calculated value of 130 mb is still a factor of 0.6 smaller than the reported experimental 
cross section. At the same time the cross section of the one-phonon states has become a 
factor 1.7 larger than the empirical value. This factor becomes 1.5 when the coupling to higher 
multiphonon states is included according to the standard Poisson distribution for the excitation 
probabilities p!7| . 

The main shortcoming of the above discussed theoretical scheme to treat the DGDR, when 
the DGDR states are obtained by folding of the fine structure of two GDR's, is the fact that 
the DGDR states obtained this way are not eigenstates of the used microscopic Hamiltonian. 
To overcome this shortcoming another calculations have been performed in which two-phonon 
[1^ (g) 1^] DGDR states are coupled directly to more complex ones ||107|, |108|, |109||. From rather 



general arguments [Q, the most important couplings leading to real transitions of the double 
giant resonances and thus to a damping width of these modes are to configurations built out 
by promoting three nucleon across the Fermi surface. That is, configurations containing three 
holes in the Fermi sea and three particles above the Fermi surface {3p3h configurations). We 
use the wave function ( |238| ) to describe the DGDR states and their coupling to Iplh and to 
3p3h doorway configurations. 

The spectrum of excited states which form the DGDR is obtained by solving the secular 
equation (|253|) and the wave function coefficient S, D and T are calculated from Eq. ( p54| ). 



Pauli principle corrections, the coefficients i^j(/32a2|a2/^2) anharmonicity shifts Ac<j 



were omitted in calculations presented in Ref. ||107|| and accounted for in Ref. 



While they 

are small, they produce shifts in the energy centroid of the double giant resonance. Similar co- 
efficients appear also in connection with the term arising from the "doorway states" containing 
three phonons in Eq. ( |238| ). We have neglected them because they again are small and further- 
more act only in higher order as compared to the previous term, in defining the properties of 
the double giant dipole resonance. Finally, the corresponding i^'-coefficient associated with the 
first term in Eq. ( ^38| ) is proportional to the number of quasiparticles present in the ground 
state of the system, a quantity which is assumed to be zero within linear response theory. 

In keeping with the fact that the Q- value dependence of the Coulomb excitation amplitude is 
rather weak at relativistic energies fl^ , the cross section associated with the two-step excitation 
of the double giant dipole resonance is proportional to 



[BiEl) X BiEl)] = \J2{%,,,,.\MiElW,l) . (VP- |A^(El)|v[/^,, 



(269) 



"2/32 



1 + ^a2,l32 a'/3' 



where Ma = (Qa||A^(-El)||0^^ ) is the reduced matrix element of the _El-operator which acting 
on the ground state \)ph excites the one-phonon state with quantum numbers a = i). 

Making use of the elements discussed above we calculated the distribution of the quantity 
Eq. ( |269|) over the states Eq. ( |238|) in ^^^Xe. We considered only J'" =0^^ and components 
of the two-phonon giant dipole resonance. As already discussed above its =1+ compo- 
nent cannot be excited in the second order perturbation theory and is sufficiently quenched in 
coupled-channels calculation. The fifteen configurations = {0(2,^2} displaying the 

largest [B{E1) x B{E1)] values were used in the calculation. They are built up out of the 
five most collective RPA roots associated with the one-phonon giant dipole resonance carrying 
the largest B[E1) values and exhausting 77% of energy weighted sum rule (EWSR). Two- 
phonon states of collective character and with quantum numbers different from 1~ lie, as a 
rule, at energies few MeV away from the double giant dipole states and were not included in 
the calculations. The three-phonon states {a3/?373} were built out of phonons with angular 
momentum and parity 1~, 2"'', 3~ and 4+. Only those configurations where either 03, j3z or 73 
were equal to 02 or /?2 were chosen. This is because other configurations lead to matrix elements 
Ualpl^^{J) of the interaction, which are orders of magnitude smaller than those associated with 
the above mentioned three-phonon configurations, and which contain in the present calculation 
5742 states up to an excitation energy 38 MeV. The single particle continuum has been approx- 
imated in the present calculation by quasibound states. This approximation provides rather 
good description of the single GDR properties in ^^^Xe. This means that our (2p2/i)fi-xi-i 



spectrum is also rather complete for the description of the DGDR properties although it is 
located at higher energies. 

If one assumes a pure boson picture to describe the phonons, without taking into account 
their fermion structure, the three-phonon configurations omitted in the present calculation do 
not couple to two-phonon states under consideration. Furthermore, although the density of 3p3h 
configurations is quite high in the energy region corresponding to the DGDR, a selection of the 
important doorway configurations in terms of the efficiency with which configurations couple 
to the DGDR, can be done rather easily. The above considerations testify to the advantage of 
employing a microscopic phonon picture in describing the nuclear excitation spectrum, instead 
of a particle-hole representation. One can more readily identify the regularities typical of the 
collective picture of the vibrational spectrum, and still deal with the fermion structure of these 
excitations. As far as the one-phonon term appearing in Eq. ( P38|) is concerned, essentially all 
phonons with angular momentum and parity O"*" and 2"*" were taken into account within the 
energy interval 20-40 MeV. 

A rather general feature displayed by the results of the present calculation is that all two- 
phonon configurations of the type building the DGDR in the "harmonic" picture 
are fragmented over a few MeV due to the coupling to 3p3h "doorway states" . Fragmentation 
of the most collective one is presented in the bottom part of Fig. For a comparison the 
fragmentation of the most collective one-phonon 1~ configuration due to the coupling to 2p2h 
"doorway states" is plotted in the top part of the same figure. The results have been averaged 
with the aid of a Breit-Wigner distribution of width 0.2 MeV. The maximum amplitude with 
which each two-phonon configuration enters in the wave function ( |238| ) does not exceed a 
few percent. Two-phonon configurations made out of two different 1" phonons are fragmented 
stronger then two-phonon configurations made out of two identical 1~ phonons. This in keeping 
with the fact that, as a rule, states of the type l^i'} with i ^ i' are less harmonic 
than states with i = i' and consequently are coupled to a larger number of three-phonon 
configurations. 

In Figs. [36|b-c, the [B{E1) x B{E1)] quantity of Eq. (|269|) associated with Coulomb ex- 
citation of the almost degenerate J'^ = and J'" = 2"*" components of double giant dipole 
resonance are shown. For comparison, the B{E1) quantity associated with the Coulomb ex- 
citation of the one-phonon giant dipole resonance is also shown in Fig. p6|a. The reason why 
the two angular momentum components of the DGDR are almost degenerate can be traced 
back to the fact that the density of one-phonon configurations to which the DGDR couple and 
which are different for = 0^ and J'^ = 2+ type states is much lower than the density of 
states associated with 3p3h "doorway states", density of states which is the same in the present 
calculation for the two different angular momentum and parity. Effects associated with the 
J-dependence of the Kj and Aj coefficients are not able to remove the mentioned degeneracy, 
because of the small size of these coefficients. These coefficients can also affect the excitation 
probability with which the J'^ = 0+ and J'" = 2+ states are excited (cf. Eq. ( p69| )). The 
effect however is rather small, leading to a decrease of the order of 2-3% in both cases. The 
J-degeneracy would be probably somehow broken if one goes beyond a one-boson exchange 
picture in the present approach of interaction between different nuclear modes. The next order 
term of interaction would couple the DGDR states to many other 3p3h configurations, not in- 
cluded in the present studies, some of these 3p3h configurations would be different for different 
J'^ values. Unfortunately, such calculation is not possible the present moment. 

The calculated excitation functions displayed in Figs. |36|b-c yield the following values for 



Table 15: Position, width and the ratio values R, Eq. ( p70|) , for of J = 0+ and 7 = 2+ 
components of the DGDR in respect to the ones of the single GDR in ^^^Xe. The third row 
corresponds to pure harmonic picture. 
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the centroid and width of the DGDR in ^^exe: < Eo+ > = 30.68 MeV and ro+ = 6.82 MeV for 
the 0+ component of the DGDR and < E2+ > = 30.71 MeV and T2+ = 6.84 MeV for the 2+ 
component. These values have to be compared to < Ei- > = 15.40 MeV and Fi- = 4.72 MeV 
for the single GDR in this nucleus from our calculation. The correspondence between these 
values is presented in Table O in comparison with the prediction of the harmonic model. Also 
shown is the ratio 

E.j:.A%n2^)\MiEi)\^';i) . (VP- |^(Ei)|vi/^ 



R 
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(270) 



between the two-step excitation probability of the DGDR normalized to the summed excitation 
probability of the one-phonon GDR. The numerical results lie quite close to the predictions of 
the harmonical model (see also a discussion of this problems in Ref ||110|| ). While the on- 
the-energy-shell transitions are easier to identify and calculate properly, off-the-energy shell 
corrections are considerably more elusive. In fact, it may be argued that the calculated shift 
of the energy centroid of the DGDR with respect to that expected in the harmonic picture is 
somewhat underestimated, because of the limitations used in selecting two-phonon basis states 
used in the calculation. Our calculated value AE = 2{Eqj^^) — (-Edgdr) shown in Table |15] can be 
compared to the ones in ^°Ca ||111|| and ^°^Pb in calculations with Skyrme forces. One of the 
purposes of the last calculations was to consider the anharmonic properties of the DGDR with 
the wave function which includes collective Iplh and 2p2h states. Thus, an interaction not only 
between the two-phonon DGDR states, [1^ x 1^] among themselves, but with other two-phonon 
states made up of collective 2+ and 3^ phonons was taken into account. The reported value of 
AE in these studies is of the order of —200 keV in consistency with our results. It should be 
pointed out that the calculation with Skyrme forces also yield somewhat larger anharmonicity 
shifts for low- lying two-phonon states as compared to the QPM calculations [|112|| . The most 
complete basis of the 2p2h configurations has been used in the second-RPA calculations of 
the DGDR properties in ^°Ca [97| and ^''^Pb |98] which includes not only "collective phonons" 
but non-collective as well. The authors of Refs. ||9^, ^ obtained the values of AE equal to 
—670 (—40) and —960 (—470) keV for 0^ and 2+ components of the DGDR, respectively, in 
^°Ca (^^'^Pb). Recently, the problem of anharmonicity for the DGDR has been also studied 
within macroscopic approaches in Refs. [^, |113|| . In Ref. [j76| it has been concluded that the 
^-dependence of it should be as A~^^^ while in Ref. [ p.l3| | it is A~'^^^ in consistency with Ref. [Q. 

The fragmentation of the DGDR due to the coupling to three-phonon configurations has 
been also calculated in ^°^Pb 111 0911 • The fine structure of the GDR as a result of interaction with 



two-phonon 1~ configurations in this nucleus is presented in Fig. [37|a. The [B{E1) x B{E1)] 
values for the DGDR states described by the wave function which includes two- and three- 
phonon configurations are plotted in Fig. |37|b. In this calculation we have used the same basis 



of six the most collective one-phonon 1~ states for the GDR and twenty-one the most collective 
two-phonon [1^ x 1~] states for 0"'' and 2"*" components of the DGDR as in the coupled- channels 
calculation in the previous subsection (see, Fig. For a description of the GDR width a 
coupling to 1161 two-phonon 1~ configurations was taken into account. In calculation of the 
DGDR strength distribution we have neglected the interaction with one-phonon configurations 
and Pauli principle corrections since these effects are weaker in double-magic nucleus ^osp]^ 
as compared to the ones in ^^^Xe. As a results the 0"*" and components of the DGDR are 
completely degenerated in this calculation. The DGDR width is determined by the coupling of 
the selected 21 two-phonon configurations with 6972 three-phonon ones and is very close to 
times the width of the GDR. This is a natural result for a folding of two independent phonons in 
microscopic treatment of the problem. As already discussed in Sect.. 3.3, when damping width 
of giant resonances is described phenomenologically by Breit-Wigner strength distribution one 
obtains the value 2 for the quantity r = rDQDj,/rGDR- On the other hand, when the Gaussian 
strength distribution is used, it yields the value r = ^/2. This is due to the different behavior 
of the wings of the abovementioned strength functions at the infinity. In a microscopic picture, 
collective resonance state(s) couples to some finite number of doorway configurations and the 
strength distribution, as a result of this coupling, is always concentrated in a definite energy 
region. It results in r = ^/2. 

The square of the amplitude cl'P^^j^j-j^^mXE) of the Coulomb excitation of one-phonon res- 
onances in RHIC in the first order perturbation theory has a smooth exponential energy de- 
pendence. This rather simplifies a calculation of the excitation cross section in RHIC of the 
states of Eq. ( p38| ) which form single giant resonances. Although a large number of the states 



of Eq. (238), the giant resonance excitation cross section in this reaction can be easily calcu- 



lated as a product of the B{E1) values of each state, presented in Figs. p6| a and |37|a, and an 
interpolated value of the tabulated function jvj.(-E)p aX E = Eu^ . The cross sections 

of the GDR and GQR excitation in ^^^Xe (see. Fig. |3^) have been calculated this way. A sim- 
ilar procedure may be applied for calculation of the cross section of the DGDR(z/o+_2+) states 
excitation via the GDR(z/i-) states. In the second order perturbation theory it equals to 



Y,A{E,^_,E,^^^^^) < l-(z/i-)||El||0+,. > 



X < [r® l-](z/o+,2+)||^l||r(^^i-) > I' (271) 

where A{Ei, E2) is the reaction amplitude which has a very smooth dependence on both argu- 
ments. This function was tabulated and used in the final calculation of the DGDR Coulomb 
excitation cross section in relativistic heavy ion collisions. 

Let us consider the excitation of the DGDR in the projectile for a ^osp]^ (640A MeV) 
_j_208p|^ collision, according to the experiment in Ref. 153, and use the minimum value of 



the impact parameter, b= 15.54 fm, corresponding to the parameterization of Ref. [Q. The 
cross section for Coulomb excitation of the DGDR is presented in Fig. ^ by the short-dashed 
curve as a strength function calculated with an averaging parameter equal to 1 MeV. The 
contribution of the background of the two-phonon 1~ states to the total cross section is shown 
by a long-dashed curve in the same figure. It was calculated in first order perturbation theory. 
The role of the background in this reaction is much less important than in photoexcitation 
studies. First, it is because in heavy ion collisions we have a special mechanism to excite 
selected two-phonon states in the two-step process. Second, the Coulomb excitation amplitude 
is exponentially decreasing with the excitation energy, while the ii^l-photoexcitation amplitude 



is linearly increasing. Nonetheless, Fig. ^ shows that the direct excitation of two-phonon 1~ 
states cannot be completely excluded from consideration of this reaction. Integrated over the 
energy interval from 20 to 35 MeV these states give a cross section of 50.3 mb which should 
be compared to the experimental cross section in the DGDR region for the ^°^Pb (640A MeV) 
_l_208p]-, reaction which is equal to 380 mb Q . 

The solid line in Fig. PSI is the sum of DGDR and two-phonon background excitations in 
relativistic heavy ion collisions. The first and second moments of excitation functions, displayed 
by the short-dashed and solid curves in Fig. indicate that the centroid of the total strength 
is 200 keV lower and the width is 16% larger than the same quantities for the pure DGDR. We 
point out that this 200 keV shift is even somewhat larger than the one due to the anharmonicities 
studied in ^ospb 171. 



Direct excitation of the two-phonon 1~ states in ^°^Pb (640A MeV) +^°'^Pb reaction was also 
investigated in Ref. |^ in calculation with Skyrme forces. The reported effect (a difference 
between 5.07 and 3.55 mb for 22 < Ej. < 28 MeV) is much weaker than in our calculation 
because of a rather limited two-phonon space. Another source of the DGDR enhancement in 
7^ is due to anharmonicity effects. We also checked the last by coupling one-phonon GDR 



states to (the most important) 1200 two-phonon 1~ states in the DGDR region. Due to the 
constructive interference between one- and two-phonon states at DGDR energies we got an 
additional enhancement of 24 mb, which is again larger than the difference between 6.42 and 
3.55 mb obtained in Ref. [73| for the same reason. 

The absolute value of the total cross section of the DGDR excitation in RHIC in ^osp]^ jg 
somewhat small in our calculation (cf. Table |T3|) as compared to experimental findings. For 
example, the experimental value of the total DGDR excitation in the reaction ^osp]-, (^g40A MeV) 
_j_208p|^^ for which the calculations have been performed, equals to 0,38(4) b. As mentioned 



above our chosen six 1 states exhaust only 90.6% of EWSR while the photo-neutron data 



indicate that this value equals to 122%. Due to this underestimate of exhaust of the EWSR 
the cross section for the DGDR excitation reduces more strongly than the one for the single 
GDR. This is because the GDR cross section is roughly proportional to the total B{E1) value 
while for the DGDR it is proportional to the square of it. 

If we apply a primitive scaling to obtain the experimental value 122% of EWSR the ratio 
R = cr(DGDR)/c"{GDR), the last line of Table |T3|, changes into 0.096 and 0.101 for the coupled- 
channels calculation and for the perturbation theory, respectively. The experimental findings 
5^ yield the value Rexp = 0.116±0.014. The reported disagreement Rexpl Rcaic =1.33±0.16 



is the result of a comparison with Rcaic obtained within a folding model, assuming 122% of the 
EWSR. We get a somewhat larger value of Rcaic (taking into account our scaling procedure) 
because the B{E1) strength distribution over our six 1~ states is not symmetrical with respect 
to the centroid energy, -Egdr^ the lower part is enhanced. A weak energy dependence in the 
excitation amplitude, which is also squared for the DGDR, enhances the DGDR cross section 
for a non-symmetrical distribution with respect to the symmetrical one, or when the GDR is 
treated as a single state. The effect of the energy dependence is demonstrated for a single GDR 
in the top part of Fig.^ where the excitation cross sections are compared to the B{E1) strength 
distribution. It produces a shift to lower energies of the centroid of the GDR and the DGDR 
cross sections with respect to the centroid of the B{E1) and the [B{E1) x B{E1)] strength 
distribution, respectively. In our calculation this shift equals to 0.26 MeV for the GDR and to 
0.33, 0.28 MeV for the DGDR within coupled-channels and perturbation theory, respectively. 
Of course, this scaling procedure has no deep physical meaning but we have included this 



discussion to indicate that the disagreement between experiment and theory for the DGDR 
excitation cross sections in ^°®Pb reached the stage when theoretical calculations have to provide 
a very precise description of both the GDR and the DGDR to draw up final conclusions. 

The situation with the absolute values of cross sections of the DGDR excitation in ^'^^Xe 
in RHIC is much less clear than in ^°®Pb. The experimental value for the reaction ^^^Xe 
(700A MeV) +^''^Pb is reported to be equal to 215±50 mb [^. This value is sufficiently larger 
as compared to any theoretical predictions available (cf. Table But it should be pointed 
out that a comparison between of experimental data for xenon |^0| and lead reveal some 
essential contradiction. While for 208p]^ ^^le discussed above quantity of the ratio between the 
total cross sections of the DGDR and GDR excitation Rg^pC^^P^) = 0.116±0.014, its value 
for ^^^Xe: i?exp(^^^Xe) = 0.21±0.05 Taking into account that experiments for both nuclei 
have been performed at close projectile energies (per nucleon) and the cross section of the 
GDR excitation in ^^^Xe is about three times less as compared to the one in ^°^Pb, the ratio 
-Rexpl^^^Xe) should be sufficiently smaller than Rexp(^^^Ph) and not vice versa. Probably, the 
problem with the absolute value of the DGDR excitation in ^^^Xe is related to uncertainties in 
separating of the contribution of single resonances, the characteristics of which are unknown 
experimentally for this nucleus and the results of interpolation have been used in evaluating 
of the experimental data. Recently, the experiment for ^^^Xe has been repeated by LAND 
collaboration ||114|| . The analysis of the new data in the nearest future should clear up the 
situation. 



5.4 The role of transitions between complex configurations of the 
GDR and the DGDR 

In the previous subsection considering the excitation properties of the DGDR, [B{E1) x B{E1)] 
values or excitation cross sections in RHIC in second order perturbation theory, we have taken 
into account only the transition matrix elements between simple one-phonon 1^ GDR and 
two-phonon 0+ or 2+ DGDR configurations for the second step of the excitation process 
g.s. GDR DGDR. In fact, as already discussed above these configurations couple to 
more complex ones to produce the widths of single and double resonances and in principle, ad- 
ditional transitions between complex configurations of the GDR and the DGDR, together with 
interference effects, may alter the predicted values of excitation probabilities. This problem will 
be considered in the present subsection (see, also Ref. [ |115|| ). It will be concluded that their role 
is marginal in the process under consideration although a huge amount of the E'l-strength is 
hidden in the GDR— >DGDR transition. This negative result ensures that calculations, in which 
only transitions between simple components of the GDR and DGDR are taken into account 
and which are much easier to carry out, require no further corrections. 

In microscopic approaches the strength of the GDR is split among several one-phonon 1~ 
states (due to the Landau damping). The wave function |1~ > couples to complex configura- 
tions |1^ > yielding the GDR width. We use the index a for simple configurations and the 
index (3 for complex ones, respectively. Thus, the wave function of the z*^ 1~ state in the GDR 
energy region can be schematically written as: 

|ir >= E^r"(«)li^ > +E^r"(/?)li^ > (272) 

where coefficients S^°^{a) and C^°^{P) may be obtained by diagonalizing the nuclear model 



Hamiltonian on the set of wave functions (P72|) . 



The total i?l-strength of the GDR excitation from the ground state, 



Baun{El)=J2\<^7\\El\Ks.> 



|2 



remains practically the same as in the one-phonon RPA calculation because the direct excitation 
of complex configurations from the ground state is a few order of magnitude weaker as compared 
to excitation of one-phonon states. However these complex configurations play a fundamental 
role for the width of the GDR. 

The wave function of the 2"*" component of the DGDR states can be written in the similar 
fashion: 

12; >= E 5f xi;>> 

a={ai xa2} 

+ E'^f («")|2+">+Ecr™(/^')|2j'> • (273) 

a" 13' 

In this equation we separated in the first term the [1^ x 1^] DGDR configurations from other 
two-phonon configurations (second term) and complex configurations (the last term) . The same 
equation as ( p73| ) is valid for the 0"*" DGDR states. 

The total ii^l-transition strength between the GDR and DGDR, 

EEI<2no+)/i|i5;i||ir>P, 

is much larger as compared to that for the GDR excitation, \ < lj~ll-E'l||0^5 > p, from 
the ground state. This is because the former includes transitions not only between simple 
GDR and DGDR states but also between complex configurations as well. The enhancement 
factor should be the ratio between the density of simple and complex configuration in the GDR 
energy region. But in the two-step excitation process the sum over intermediate GDR states in 
Eq. ( p7^ ) reduces the total transition strength for g.s. GDR DGDR to ~ 2 ■ |5gdr(^1)P 
(the factor 2 appears due to the bosonic character of the two phonons which also holds if Landau 
damping is taken into account). To prove this let us consider the excitation probability of the 
DGDR 



(274) 



where the index i labels intermediate states belonging to the GDR, and 0\jl{wh)-^j2{M2) 
first-order El excitation amplitude for the transition Ji(Mi) — > ^2(^2) in a collision with 
impact parameter h. For each state, J and M denote the total angular momentum and the 
magnetic projection, respectively. 

The amplitude aJi^(Mi)^j2(M2) given by 

af[^^^^j^^,,^^{E, h) = iJ,MMJ2M2)x < J2\\E1\\J, > fEi^E, b) . (275) 

It is a product of the reduced matrix element < J2||ii^l|| Ji > for the i?l-transition between the 
states Ji(Mi) and 72(^^2) which carries nuclear structure information and the reaction function 
fEi(ij.){E, b). The latter depends on the excitation energy, charge of the target, beam energy, 
and is calculated according to Ref. [Ol. Except for the dependence on the excitation energy, it 



does not carry any nuclear structure information. The cross section for the DGDR is obtained 
from Eq. ( |274|) by integration over impact parameters, starting from a minimal value bmin to 



infinity. This minimal value is chosen according to Ref . [ 1^ . 



Now we substitute the wave functions of the GDR and DGDR states given by Eqs. ( |272| , p73| ) 



in expression ( |274| ). We obtain two terms. The first one corresponds to transitions between 
simple GDR and DGDR states (after the GDR is excited from the ground state through its 
simple component): 

= E E Sr^a)fEiUE,, b) < l;||i?l||0+,. > (276) 

i aa'a 

X Sr^a')S^^°^a)fEU,'){Ef - E.,, b) < [l'^ x l;J/||i?l||l;, > 

and the second one accounts transitions between complex configurations in the wave functions 
of Eqs. ( p7^ , p73| ): 

B,,' = E E ^r»/i^i(.)(^., &)<l;ll^l||0+..> (277) 

i aa'PP' 

X Cr^{(3)C]^-^{(3')fEii,'){Ef - E,, b) < [i;, X > V,[a'x/.] . 

The second reduced matrix element in the above equations is proportional to the re- 
duced matrix element between the ground state and the simple one-phonon configuration (see, 
Eq. (g63)). 

For a given impact parameter 6, the function fEi{ii.) {E, b) can be approximated by a constant 
value fEi{p) !§] for the relevant values of the excitation energies. Then the energy dependence 
can be taken out of summations and orthogonality relations between different components of 
the GDR wave functions can be applied [|107|| . The orthogonality relations between the wave 



functions imply that 

E^r»c^r"(/5) = o. (278) 

i 

This means that the term 5^^/ vanishes. The term A^^/ summed over projections and all 
final states yields a transition probability to the DGDR, Pogdr(-E/, b), which is proportional to 
2 - |Sgdr(£'1) p in second order perturbation theory. This argument was the reason for neglecting 
the term B^^r in previous calculations of DGDR excitation where the coupling of simple GDR 
and DGDR states to complex configurations was taken into account. 



In Fig. ^ we plot the value of xei{E) = 2tt J db bY^fj, |/bi(/x)(-E, b)\ as a function of energy 
calculated for the ^'^^Pb (640A MeV) + ^osp]^ reaction. This value corresponds to ctgdr if 
Bqj^^[E1) = 1. The square in this figure indicates the location of the GDR in ^'^^Pb. This 
figure demonstrates that the function xei{E) changes by 60% in the GDR energy region. The 



role of this energy dependence for other effects has been considered in Refs. |T9| , [7^ . Taking 



into account that one-phonon l^, configurations are fragmented over a few MeV |108|, when a 



sufficiently large two-phonon basis is included in the wave function given by Eq. ( |272| ), the role 



of the B^^i term in the excitation of the DGDR should be studied in more detail. 

To accomplish this task we have performed firstly a simplified calculation in which we used 
the boson type Hamiltonian: 

H = Y^^aQiQa + Y.^pQ^pQp + E U1{QiQp + h.c.) (279) 

a P a,l3 

where is the phonon creation operator and Ua is the energy of this one-phonon configuration; 
is the operator for creation of a complex configuration with energy a;^ and is the matrix 



element for the interaction between these configurations. We have assumed that the energy 
difference between two neighboring one-phonon configurations is constant and equals to Aa;. 
An equidistant spacing with the energy Auo was assumed for the complex configurations. We 
also have used a constant value U for the matrix elements of the interaction. The -Bgor(-EI) 
value was distributed symmetrically over one-phonon configurations. Thus, the free parameters 
of this model are: Ac^;, Aa), f/, the number of one-phonon and complex configurations, and the 
distribution of the Bqj^^{E1) value over the simple configurations. The only condition we want 
to be satisfied is that the energy spectrum for the GDR photoexcitation is the same as the one 
known from experiment. 

After all parameters are fixed we diagonalize the model Hamiltonian of Eq. ( |279|) on the set 
of wave functions of Eq. (|7|) for the GDR and on the set of Eq. (|273|) for the DGDR. The 
diagonalization procedure yields the information on eigen energies of the 1~ GDR states and 
on the coefficients Sf^^{a) and Cf°^{(3), respectively. One also obtains information on eigen 
energies of the 2j or Oj DGDR states and the coefficients S^'^^^{a) and Cj'~'°" (/?'), respectively. 
With this information we are able to study the role of the 5^^' term in the excitation of the 
DGDR in RHIC. 

The big number of free parameters allows an infinite number of suitable choices. In fact, 
not all of the parameters are really independent. For example, the increase in the number 
of simple or complex configurations goes together with the decreasing of the value of U. This 
is necessary for a correct description of the GDR photoabsorption cross section. This makes 
it possible to investigate the role of the B^^i term in different conditions of weak and strong 
Landau damping and for different density of complex configurations. In our calculations we 
vary the number of collective simple states from one to seven and the number of complex 
configurations from 50 to 500. The value of U then changes from about 100 to 500 keV. The 
results of one of these calculations for the excitation of the 2"*" component of the DGDR in 
208pj^ (640A MeV) + ^osp]^ collisions are presented in Fig. For a better visual appearance 
the results are averaged with a smearing parameter equal to 1 MeV. The dashed curve shows 
the results of a calculation in which a^^^^^^E) = <7ogdk{E) J db b\A^^'\'^ and the results of 
another one in which a^^j^^{E) = crooDRiE) J db b |v4^^/ + B^^/f are represented by a solid 
curve. 

Our calculation within this simple model indicates that the role of the -B^^/ term in second 
order perturbation theory is negligibly small, although the total B{E1) strength for transitions 
between complex GDR and DGDR configurations, considered separately, is more than two 
orders of magnitude larger than the ones between simple GDR and DGDR configurations. 
The value Aa = {ct^gt^^ - o-dgdhJ/c^dgdr^ where ctdgdr'^^ = ^ (Tuitt^\E)dE, changes in these 
calculations from 1% to 2.5%. The results practically do not depend on the number of complex 
configurations accounted for. The maximum value of Ao" is achieved in a calculation with 
a single one-phonon GDR state (no Landau damping). This is because the value of U is 
the larger in this case and the fragmentation of the one-phonon state is stronger. Thus, in 
such a situation, the energy dependence of the reaction amplitude modifies appreciably the 
orthogonality relations. But in general the effect is marginal. 

We also performed a calculation with more realistic wave functions for the GDR and DGDR 
states taken from our studies presented in the previous subsection. These wave functions include 
6 and 21 simple states for the GDR and DGDR, respectively. The complex configurations are 
two-phonon states for the GDR and three-phonon states for the DGDR. The value Act equals 
in this realistic calculation to 0.5%. This result is not surprising because realistic calculations 



with only two-phonon complex configurations, and a limited number of them, somewhat un- 
derestimate the GDR width which is crucial for the modification of the orthogonality relations. 

We have proved that the transitions between complex GDR and DGDR configurations 
within second-order perturbation theory for the DGDR excitation in RHI collisions play a 
marginal role in the process under consideration and it is sufficient to take into account only 
transitions between the ground state and one-phonon GDR and two-phonon DGDR configura- 
tions. 



5.5 The DGDR in deformed nuclei 

The possibility to observe two-phonon giant resonances in deformed nuclei with the present 
state of art experimental techniques is still questionable. This is mainly due to the fact that 
one has to expect a larger width of these resonances as compared to spherical nuclei. Also, the 
situation with the low-lying two-phonon states in deformed nuclei is much less clear than in 
spherical ones. 

The first experiment with the aim to observe the double giant dipole resonance (DGDR) in 
238^ in relativistic heavy ion collisions (RHIC) was performed recently at the GSI/SIS facility 
by the LAND collaboration ||1 14|| . It will take some time to analyze the experimental data and 



to present the first experimental evidence of the DGDR in deformed nuclei, if any. Thus, we 
present here the first theoretical predictions of the properties of the DGDR in deformed nuclei 
based on microscopic study [ |116|| . The main attention will be paid to the width of the DGDR 
and its shape. 

In a phenomenological approach the GDR is considered as a collective vibration of protons 
against neutrons. In spherical nuclei this state is degenerate in energy for different values of 
the spin J = 1~ projection M = 0, ±1. The same is true for the 2"*" component of the DGDR 
with projection M = 0, ±1, ±2. In deformed nuclei with an axial symmetry like ^^^U, the GDR 
is spit into two components I'^{K) = 1~(0) and I'^{K) = 1~(±1) corresponding to vibrations 
against two different axes. In this approach one expects a three-bump structure for the DGDR 
with the value K = 0, K = ±1 and K = 0, ±2, respectively (see. Fig. ^). Actually, the GDR 
possesses a width and the main mechanism responsible for it in deformed nuclei is the Landau 
damping. Thus, the conclusion on how three bumps overlap and what is the real shape of the 
DGDR in these nuclei, i.e., either a three-bump or a fiat broad structure, can be drawn out 
only from some consistent microscopic studies. 

We use in our calculations for 238u the parameters of Woods-Saxon potential for the average 



field and monopole pairing from Ref. ||117|| . They were adjusted to reproduce the properties 
of the ground state and low-lying excited states. The average field has a static deformation 
with the deformation parameters /?2 = 0.22 and /?4 = 0.08. To construct the phonon basis for 
the K = and K = ±1 components of the GDR we use the dipole-dipole residual interaction 
(for more details of the QPM application to deformed nuclei, see e.g. Ref. '^^^ strength 



parameters of this interaction are taken from Ref. ||118|| where they have been fitted to obtain 
the centroid of the B{E1, 0^^ 1~{K = 0, ±1)) strength distribution at the value known from 



experiment ||119|| and to exclude the center of mass motion. In this approach, the information on 



the phonon basis (i.e. the excitation energies of phonons and their internal fermion structure) 
is obtained by solving the RPA equations. For electromagnetic ii^l-transitions we use the values 
of the effective charges, e^jf^ = eN{—Z)/A to separate the center of mass motion. 

The results of our calculation of B{E1) strength distribution over 1a'=o(0/ \'^k=±i{^'] 



GDR states are presented in Fig. together with experimental data. The index i in the 
wave function stands for the different RPA states. All one-phonon states with the energy lower 
than 20 MeV and with the B{E1) value larger than 10~^ fm^ are accounted for. Their total 
number equals to 447 and 835 for the K = and K = ±1 components, respectively. Only 
the strongest of them with B{E1) > 0.2 fm^ are shown in the figure by vertical lines. Our 
phonon basis exhausts 32.6% and 76.3% of the energy weighted sum rules, 14.8 ■ NZ/A fm^ 
MeV, by the K = {] and K = ±1 components, respectively. For a better visual appearance 
we also present in the same figure the strength functions averaged with a smearing parameter, 
which we take as 1 MeV. The short (long) dashed-curve represent the K = {K = ±1) 
components of the GDR. The solid curve is their sum. The calculation reproduces well the 
two-bump structure of the GDR and the larger width of its -RT = ±1 component. The last is 
consistent with the experiment [p.l9|| which is best fitted by two Lorentzians with widths equal 
to Fi = 2.99 MeV and F2 = 5.10 MeV, respectively. The amplitudes of both maxima in the 
calculation are somewhat overestimated as compared to the experimental data. This happens 
because the coupling of one-phonon states to complex configurations is not taken into account 
which can be more relevant for the K = ±1 peak at higher energies. But in general the coupling 
matrix elements are much weaker in deformed nuclei as compared to spherical ones and the 
Landau damping describes the GDR width on a reasonable level. 

The wave function of the O"'' and 2"^ states belonging to the DGDR are constructed by 
the folding of two 1~ phonons from the previous calculation. When a two-phonon state is 
constructed as the product of two identical phonons its wave function gets an additional factor 
l/-\/2. The l"*" component of the DGDR is not considered here for the same reasons as in 
spherical nuclei. The anharmonicity effects which arise from interactions between different 
two-phonon states are also not included in the present study. 

The folding procedure yields three groups of the DGDR states: 



«) I[1k=o(^i) ® lx=o(^2)]o+^^„2+^o)> 

b) |[l7<=oW®lx=±i(^')]2+^^,) and 

c) \[lK=^^{^[)®lK=^^moi^^,2i^^^ J■ (280) 

The total number of non-degenerate two-phonon states equals to about 1.5 ■ 10®. The energy 
centroid of the first group is the lowest and of the last group is the highest among them. So, 
we also obtain the three-bump structure of the DGDR. But the total strength of each bump is 
fragmented over a wide energy region and they strongly overlap. 

Making use of the nuclear structure elements discussed above, we have calculated the ex- 
citation of the DGDR in ^^^U projectiles (0.5 GeV-A) incident on ^^°Sn and ^osp]^ targets, 
following the conditions of the experiment in Ref. ||114|| . These calculations have been per- 
formed in second order perturbation theory P|, in which the DGDR states of Eq. (|280|) are 
excited within a two-step process: g.s.— >GDR— »-DGDR. As intermediate states, the full set of 
one-phonon l~^^Q{i)'j and states was used. We have also calculated the GDR exci- 

tation to first order for the same systems. The minimal value of the impact parameter, which 

is very essential for the absolute values of excitation cross section has been taken according to 

1/3 



6™„ = i.28-(Ar + Ay3). 



The results of our calculations are summarized in Fig. ^ and Table |T6[ In Fig. ^ we 
present the cross sections of the GDR (part a) and the DGDR (part b) excitation in the ^^^U 
(0.5 GeV-A) + ^°^Pb reaction. We plot only the smeared strength functions of the energy dis- 
tributions because the number of two-phonon states involved is numerous. The results for ^^^U 



Table 16: The properties of the different components of the GDR and the DGDR in ^^^U. The 
energy centroid E^., the second moment of the strength distribution m2 in RHIC, and the cross 
sections a for the excitation of the projectile are presented for: a) ^'^^U (0.5A GeV) + ^^°Sn, 
and b) 238U (0.5A GeV) + ^o^Pb. 
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(0.5 GeV-A) + Sn reaction look very similar and differ only by the absolute value of cross sec- 
tions. In Table |l^ the properties of the GDR and the DGDR, and their different K components 



are given. The energy centroid and the second moment, m2 = y Z)fc • (-^fc — -^c)^/ Sfc ^k-, 
of the distributions are averaged values for the two reactions under consideration. 

The two-bump structure can still be seen in the curve representing the cross section of the 
GDR excitation in ^^^U in RHIC as a function of the excitation energy. But its shape differs 
appreciably from the B{E1) strength distribution (see Fig. |45| a in comparison with Fig. ^). 
The reason for that is the role of the virtual photon spectra. First, for the given value of 
the excitation energy and impact parameter it is larger for the K = ±1 component than that 
for the K = one (see also the first two lines in Table 0). Second, for both components 
it has a decreasing tendency with an increase of the excitation energy 0. As a result, the 
energy centroid of the GDR excitation in RHIC shifts by the value 0.7 MeV to lower energies 
as compared to the same value for the B{E1) strength distribution. The second moment m2 
increases by 0.2 MeV. 

The curves representing the cross sections of the excitation of the K = ±1 and K = ±2 
components of the DGDR in ^ss^j {{^hIC have typically a one-bump structure (see the curves 
with squares and triangles in Fig. ^3|b, respectively). It is because they are made of two- 
phonon 2+ states of one type: the states of Eq. ( p80| b) and Eq. (|280| c), respectively. Their 
centroids should be separated by an energy approximately equal to the difference between the 
energy centroids of the K = and K = ±1 components of the GDR. They correspond to 
the second and the third bumps in a phenomenological treatment of the DGDR. The K = 
components of the DGDR include two group of states: the states represented by Eq. ( |280| a) and 
those of Eq. ( p80| c). Its strength distribution has two-bumps (see the curve with circles for the 
2~^{K = 0) and the dashed curve for the 0~^{K = 0) components of the DGDR, respectively). 
The excitation of the states given by Eq. ( |280| a) in RHIC is enhanced due to their lower 
energies, while the enhancement of the excitation of the states given by Eq. ( p80| c) is related 
to the strongest response of the K = ±1 components to the external El Coulomb field in both 
stages of the two-step process. 

Summing together all components of the DGDR yields a broad one-bump distribution for 



the cross section for the excitation of the DGDR in ^^^U, as a function of excitation energy. 
It is presented by the sohd curve in Fig. |43|b. Another interesting resuh of our calculations is 
related to the position of the DGDR energy centroid and to the second moment of the DGDR 
cross section. The centroid of the DGDR in RHIC is shifted to the higher energies by about 
0.8 MeV from the expected value of two times the energy of the GDR centroid. The origin 
for this shift is in the energy dependence of the virtual photon spectra and it has nothing to 
do with anharmonicities of the two-phonon DGDR states. In fact, the energy centroid of the 
B[E1, g.s. 1^) X B{E1, Ij" ^DGDRj) strength function appears exactly at twice the 
energy of the centroid of the B{E1, g.s. — >GDR) strength distribution because the coupling 
between different two-phonon DGDR states are not accounted for in the present calculation. 
The same shift of the DGDR from twice the energy position of the GDR in RHIC also takes 
place in spherical nuclei. But the value of the shift is smaller there because in spherical nuclei 
the GDR and the DGDR strength is less fragmented over their simple configurations due to 
the Landau damping. But the larger value of the shift under consideration in deformed nuclei 
should somehow simplify the separation of the DGDR from the total cross section in RHIC. 

Another effect which also works in favor of the extraction of the DGDR from RHIC excita- 
tion studies with deformed nuclei is its smaller width than ^/2 times the width of the GDR, as 
observed with spherical nuclei. Our calculation yields the value 1.33 for the ratio roGDH/rGDR 
in this reaction. The origin for this effect is in the different contributions of the GDR K = 
and K = ±1 components to the total cross section, due to the reaction mechanism. It should 
be remembered that only the Landau damping is accounted for the width of both the GDR 
and the DGDR. But since the effect of narrowing of the DGDR width is due to the selectivity 
of the reaction mechanism it will still hold if the coupling to complex configurations is included 
in the calculation. 

It may be argued that the procedure of independent excitations of two RPA phonons applied 
here is not sufficient for a consistent description of the properties of the two-phonon giant 
resonances. This is true for the case of spherical nuclei where only the coupling of two GDR 
phonons to more complex, 3p3h, configurations allows one to describe the DGDR width as 
discussed above. But the typical matrix element of this coupling in deformed nuclei does not 
exceed the value of 200 keV ||120|| while in spherical nuclei it is an order of magnitude larger. It 
means that due to the coupling, the strength of each GDR RPA-phonon will fragment within 
the energy interval of 100-200 keV in deformed nuclei. The last value should be compared 
to the second moment, m2, presented in Table ^ which is the result of the Landau damping 
accounted for in our calculation. Taking into account that the reaction amplitude has very weak 
energy dependence and that mixing of different RPA phonons in the GDR wave function does 
not change the total strength ||115|| , the total cross sections of the GDR and DGDR excitation 
in RHIC will be also conserved. 
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Figure captions 



Figure 1: Experimental cross sections in arbitrary units for the excitation of ^°^Pb targets by 
^^O (22A MeV and 84A MeV) and by ^"^Ar (95A MeV), as a function of the excitation energy. 

Figure 2: A nuclear target is Coulomb excited by a fast moving projectile. The coordinates 
used in text are shown. 

Figure 3: Electric dipole number of equivalent photons per unit area d^b = 2TTbdb, with energy 
of 10 MeV, incident on ^osp]^ ^ collision with ^^O at impact parameter 6 = 15 fm, and as 
a function of the bombarding energy in MeV per nucleon. The dotted line and the dashed 
line correspond to calculations performed with the non-relativistic and with the relativistic 
approaches, respectively. The solid line represents a more correct calculation, as described in 
the text. 

Figure 4: Same as Fig. H, but for the E2 multipolarity. 

Figure 5: Same as Fig. ^ but for the Ml multipolarity. 

Figure 6: Equivalent photon numbers per unit area incident on ^''^Pb, in a collision with ^^O 
at lOOA MeV and with impact parameter 6 = 15 fm, as a function of the photon energy huo. 
The curves for the El, E2 and Ml multipolarities are shown. 

Figure 7: Total cross sections for the excitation of giant electric dipole (El) and quadrupole 
{E2) resonances in ^osp]^ means of the Coulomb interaction with ^^O, as a function of the 
laboratory energy. 

Figure 8: Ratio to the Rutherford cross section of the elastic cross section for the ^^O + ^osp]^ 
reaction at 84A MeV, as a function of the center-of-mass scattering angle. Data are from 



Ref. PI 



Figure 9: Differential cross section for the excitation of the isovector giant dipole resonance in 
Pb by means of ^^O projectiles at 84A MeV, function of the center-of-mass scattering 
angle. Data are from Ref. ||12| . 

Figure 10: Virtual photon numbers for the electric dipole multipolarity generated by 84A MeV 
^''O projectiles incident on ^"^Pb, as a function of the center-of-mass scattering angle. The solid 
curve is a semiclassical calculation. The dashed and dotted curves are eikonal calculations with 
and without relativistic corrections, respectively. 

Figure 11: Differential cross sections for the excitation of the giant dipole resonance (GDR), 
the isoscalar giant quadrupole resonance (GQRjs), and the isovector giant quadrupole resonance 
(GQRi^), in Pb for the collision 208p|3^208pb q40A MeV. The solid (dotted) [dashed-dotted] 
line is the differential cross section for the excitation of the GDR (GQRjs) [GQRjj;]. The dashed 
line is the result of a semiclassical calculation. 



Figure 12: Calculated cross section for the excitation followed by 7-decay of ^osp]^ induced 
by ^^O projectiles at 84A MeV. The photoabsorption cross section was parameterized by a 
simple Lorentzian representing the GDR, and the statistical component of the photon decay 
was neglected. The solid curve uses the formalism described in the text (Eq. p.39|) while the 
dashed curve uses a constant branching ratio for photon decay (Eq. |138|) . 



Figure 13: Cross section for the excitation of the GDR without the detection of the decay 



photon. Data are from Ref. |13 



Figure 14: Cross section for excitation followed by 7-decay of ^°^Pb by ^''O projectiles at 84A 
MeV. The solid (dashed) line includes (excludes) the Thomsom scattering amplitude. Data are 
from Ref. |13|]. 



Figure 15: The GDR Coulomb excitation probabilities as functions of the impact parameter, 
for sharp and smooth absorptions. The system is ^osp^ (640A MeV) + ^''^Pb. 

Figure 16: Nuclear excitation probabilities, as functions of the impact parameter, of the 
isoscalar giant monopole resonance (GMRj^), the GDRj^, and the GQRjs, in 208p|^ f^^. ^^^q 
collision ^ospb+aosp^ at 640A MeV. 

Figure 17: Experimental In- and 3n- removal cross sections for ^^^Au bombarded with rela- 
tivistic projectiles (from Ref. |3^) in comparison with theoretical calculations from this work 
(solid curve: "soft- spheres" model; dotted curve: " sharp-cutofF' model with from Eq. 

( 11501) ). For completeness, we also show a sharp-cutoff calculation with 6^°;^ (Eq. (|151|) ) used 
in Ref. Q. 

Figure 18: A nuclear target is Coulomb excited by a fast moving deformed projectile. Besides 
the angle 6, the orientation of the projectile also includes an azimuthal angle which can rotate 
its symmetry axis out the scattering plane. For simplicity, this is not shown, x is the angular 
position of the cm. of the projectile with respect to the target. 

Figure 19: Percent increase of the Coulomb excitation cross section of dipole states in ^°^Pb 
due to the dependence of the minimum impact parameter on the deformation. The effect is 
shown for 238]j projectiles at lOOA MeV, lA GeV and lOA GeV, respectively, and as a function 
of the deformation parameter /3. The solid (dashed) [dotted] line corresponds to an excitation 
energy of 1 (10) [25] MeV. For the actual deformation of ^°^U, /? ~ 0.3, the effect is small. 

Figure 20: Coulomb excitation cross section of a giant dipole resonance in 208p|^ ^j^g 
quadrupole-dipole interaction with lOOA MeV uranium projectiles, as a function of the de- 
formation parameter f3. These cross section are averaged over all possible orientations of the 
projectile. 

Figure 21: Experimental results for ^^^Xe projectile excitation (at 690A MeV) on a Pb target 
(squares) and a C target (circles). The spectrum for the C target is multiplied by a factor 2 
for better presentation. The resonance energies for one- and two-phonon giant resonances are 
indicated. The dashed curve reflects the results of a first-order calculation for the Pb target. 



Fig. is taken from Ref. [^ . 



Figure 22: Compilation of experimental findings with heavy ion (full symbols) and pion induced 
(open symbols) reactions for the energy, width, and cross sections of the double giant resonance. 
The data are compared to the energies and widths of the giant dipole resonance, respectively, 
and to the theoretical values of excitation cross sections. 

Figure 23: Schematic representation of the excitation of Giant Resonances, populated in heavy 
ion collisions. 

Figure 24: Time-dependence of the occupation probabilities |aoP and |aip, in a collision with 
impact parameter 6 = 15 fm. The time is measured in terms of the dimensionless variable 
r = iv-i/h) t. The system is ^ospb (640A MeV) + ^o^Pb. 

Figure 25: Excitation energy spectra of the main Giant Resonances for both Breit-Wigner and 
Lorentzian line shapes. The system is ^°^Pb (640A MeV) + ^^'^Pb. 

Figure 26: Ratio between the DGDR and the GDR cross sections in ^osp]^ _|_ 208p|^ collisions, 
as a function of the bombarding energy. 

Figure 27: Dependence of ctgdr and ctdgdr on the GDR width, treated as a free parameter. For 
details see the text. The system is ^"^Pb (640A MeV) + ^"^Pb. 

Figure 28: Cross sections of the direct photoexcitation of two-phonon configurations [l^(g'l^]2+, 
[1" (g) 2+]i- and [2+ ® 3"]i- from the groimd state in ^^^Xe and ^ospb. 

Figure 29: Photoexcitation cross section of the GDR in ^^^Xe and ^°®Pb. Calculations are 
performed: a) within one-phonon approximation and b) with taking into account of the coupling 
between one- and two-phonon configurations. Continues curves in the bottom part are the 
strength functions calculated with a smearing parameter A = 1 MeV; dashed curve corresponds 
to electromagnetic transitions to one-phonon 1~ states, solid curve - to one- and two-phonon 
1~ states. 

Figure 30: Photoneutron cross sections in ^°^Pb. Solid curve is the result of calculation with the 
wave function including one- and two-phonon terms presented with a smearing parameter A = 
1 MeV; vertical lines (in arbitrary units) - within one-phonon approximation. Experimental 
data are plotted by experimental error bars. 

Figure 31: Photo-neutron cross section for ^"'^Pb. Experimental data (dots with experimental 
errors) are from Ref. ||103|| . The long-dashed curve is the high energy tail of the GDR, the 



short-dashed curve is the CQR^ and the curve with squares is their sum. The contribution of 
two-phonon states is plotted by a curve with triangles. The solid curve is the total calculated 
cross section. 

Figure 32: The possible paths to the excitation of a given magnetic substate of the l"*" com- 
ponent of the DGDR are displayed. The transitions caused by the different projections of the 
operator Vei are shown by: (a) dashed lines for m = 0, (b) dashed-dotted lines for m = —1, 
and (c) solid-lines are m = -|-1. 



Figure 33: The electromagnetic excitation cross sections for the reaction 208p]^ (640A MeV) 
_j_208p|j calculated in coupled- channels. It is shown the excitation of the GDR (top) and the 
three components J'^ = 0+, 2+ and 1+ of the DGDR. The B{E1) strength distribution (in 
arbitrary units) over 1^ states is shown by dashed lines. For a visuality it is shifted up by 
100 keV. 

Figure 34: The cross section for Coulomb excitation of the one-phonon GDR (continuous 
curve), of the isoscalar GQR (dash-dotted), of the isovector GQR (long dashed) as well as for the 
double- phonon GDR (short dashed) are shown. They have been calculated at Eiab = 681 A MeV, 
taking into account the energy reduction of the beam in the target |4^. The one-phonon GDR 
cross section has been reduced in the figure by a factor 10. 

Figure 35: Fragmentation of the most collective a) one-phonon 1~ and b) two-phonon [1^ 
configurations in ^^^Xe due to the coupling to more complex configurations. The result are 
presented with a smearing parameter A = 0.2 MeV. 

Figure 36: a) B{E1) values for the GDR and b,c) [B{E1) x B{E1)] values Eq. ( p69| ) for the 
DGDR associated with Coulomb excitation in ^^^Xe in relativistic heavy ion collision, b) and 
c) correspond to J = 0^ and J = 2+ components of the DGDR, respectively. A smooth curve 
is a result of averaging over all states with a smearing parameter A = 0.5 MeV. See text for 
details. 

Figure 37: a) B{E1) values for the GDR and b) [B{E1) x B{E1)] values Eq. (|6|) for the 
DGDR (J = 0^ + 2+) associated with Coulomb excitation in ^osp]^ relativistic heavy ion 
collision. A smooth curve is a result of averaging over all states with a smearing parameter 
A = 0.5 MeV. 

Figure 38: The contribution for the excitation of two-phonon 1~ states (long-dashed curve) in 
first order perturbation theory, and for two-phonon 0"*" and 2+ DGDR states in second order 
(short-dashed curve). The total cross section (for 208p]^ (640A MeV) -|-^°^Pb) is shown by the 
solid curve. 

Figure 39: The energy dependence of the 208p]^ (640A MeV) -|- ^osp^ reaction function calcu- 
lated within first order perturbation theory. The square indicates the location of the GDR in 

208 pj^ 

Figure 40: The cross section for the excitation of the 2"^ component of the DGDR in the 
reaction 208p]^ (640A MeV) -|- ^°'^Pb, calculated within second order perturbation theory. The 
dashed curve shows the contribution of the £'l-transition between simple GDR and DGDR 
configurations only. The solid curve is a sum of the above result and the contribution of the 
i?l-transitions between complex GDR and DGDR configurations. See text for details. 



Figure 41: The possible paths to the excitation of a given magnetic substate of the O"*" and 2"*" 
components of the DGDR in spherical and deformed nuclei. The notations are the same as in 
Fig. |2. 



Figure 42: The B{E1) strength distribution over K = (short-dashed curve) and K = ±1 
(long-dashed curve) 1^ states in ^^^U. The sohd curve is their sum. The strongest one-phonon 
1~ states are shown by vertical lines, the ones with K = are marked by a triangle on top. 



Experimental data are from Ref. ||119 |. 



Figure 43: The strength functions for the excitation: a) of the GDR, and b) of the DGDR 
in ^^^U in the ^ss-y (^q_5A GeV) -f ^ospi^ reaction. In a), the short-dashed curve corresponds 
to the GDR {K = 0) and the long-dashed curve to the GDR {K = ±1). In b) the dashed 
curve corresponds to the DGDRo+ {K = 0), the curve with circles to the DGDR2+(-ft' = 0), the 
curve with squares to the DGDR2+ {K = ±1), and the curve with triangles to the DGDR2+ 
{K = ±2) . The solid curve is the sum of all components. The strength functions are calculated 
with the smearing parameter equal to 1 MeV. 



This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig2.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig3.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig4.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig5.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig6.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig7.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig8.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig9.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl0.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figll.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl2.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl3.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl4.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl5.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl6.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl7.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl8.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "figl9.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig20.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig21.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig22.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig23.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig24.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig25.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig26.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig27.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig28.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig29.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig30.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig31.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig32.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig33.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig34.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig35.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig36.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig37.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig38.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig39.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig40.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig41.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig42.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



This figure "fig43.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/9908003vl 



